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Abstract 

Boresight  Error  (BSE),  defined  as  the  angular  deviation  between  the  true  position 
and  the  apparent  position  of  a  target  as  indicated  by  a  radar,  is  a  very  important 
figure  of  merit  for  a  tracking  radar.  Ideally,  zero  BSE  is  desired  but  seldom  achieved. 
Hence,  a  capability  to  accurately  predict  BSE  in  the  design  phase  of  a  new  radar 
system  or  to  impact  modification!  of  an  existing  system  becomes  imperative.  Prior 
work  on  the  subject  matter  is  somewhat  sketchy  and  limited  in  scope.  Therefore,  this 
dissertation  undertakes  a  thorough  and  comprehensive  investigation  of  BSE  using  a 
systems  concept  so  that  the  final  product  is  applicable  to  a  variety  of  situations.  A 
monopulse  tracker  was  chosen  for  this  study  because  it  possesses  superior  angle 
tracking  capability  and  it  is  used  in  a  majority  of  modem  radar  systems. 

Although  theie  are  many  factors  intrinsic  and  extrinsic  to  a  radat  system  that  give 
rise  to  BSE,  the  most  significant  contributor  is  the  protective  radome.  An  incoming 
plane  wave  suffers  depolarization  and  phase  front  distortion  as  it  travels  through  the 
radome.  The  net  result  of  these  undesirable  changes  is  BSE,  resulting  in  tracking 
error  and  inaccurate  target  location  estimates.  Surface  integration  and  Geometric 
Opucs  (GO)  arc  two  methods  commonly  used  to  investigate  the  effects  of  a  tadome 
on  BSE.  Building  on  the  "consensus"  that  the  ray-trace  receive  GO  propagation 
technique  offers  the  "best"  compromise  between  accuracy  and  computational 
intensity,  this  research  effort  employed  a  GO  technique  which  greatly  expanded 
previous  ray-trace  receive  techniques  to  include:  1)  a  uniquely  Jefmed/developed 

vii; 


mathematical  description  for  each  surface  within  arbitrary  multi-layer  tapered 
radomes,  surface  descriptions  are  generated  using  a  tangent  ogive  reference  surface, 
2)  an  "ideal"  taper  function  concept  for  obtaining  optimum  BSE  prediction 
performance,  3)  a  generalized  technique  for  calculating  specular  reflection  points 
within  the  radome,  and  4)  the  total  refractive  cffec's  along  ray  propagation  paths. 
Fortran  computer  model  results  were  compared  with  limiting  case  data  (BSE  =  0*), 
published  experimental  data,  and  production  system  acceptance  test  data.  Limiting 
case  validation  was  accomplished  using  1)  single  and  multi-layer  tapered  "a:r" 
radomes  by  setting  the  relative  permittivity  and  permeability  of  each  layer  equal  to 
one,  2)  hemispheric  radomes  with  the  aper.ure  gimbal  point  located  at  the  sphere 
center,  and  3)  radome  boresight  axis  scanning  (through  the  tip).  For  these  cases, 
"system"  modeling  error  was  less  than  .06  mRad  for  all  scan  angles  and  polarizations 
of  interest.  "Excellent"  (BSE  within  ±1  mRaa)  results  were  obtained  using  a 
hemispheric  half-wave  wall  radome  witn  a  displaced  aperture  gimbal  point;  predicted 
SSE  values  were  within  ±  1  mRad  of  published  surface  integration  and  measured 
experimental  data.  Likewise,  modeled  BSE  pi  evictions  for  the  production  system 
were  within  i0.5  mRad  of  measured  data  over  a  30°  scan  range.  Validated  model 
results  were  then  used  to  determine  overall  ray  refractive  effects  on  BSE  prediction 
and  found  to  provide  only  marginal  benefit  at  the  expense  of  computational 
efficiency.  In  conclusion,  all  the  initial  goals  and  aims  of  this  research  were  not  only 
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met  but  exceeded  in  this  dissertation. 


RADOME  DEPOLARIZATION  EFFECTS  ON 
MONOPULSE  RECEIVER  TRACKING  PERFORMANCE 


I.  Introduction 

The  ability  to  accurately  locate  and  track  hostile  targets  either  actively  or 
passively  is  a  key  feature  in  the  survival  and  mission  effectiveness  of  a  vast  majority 
of  airborne  military  platforms.  Monopulse  processing  techniques  are  particularly  well 
suited  for  angle  tracking  Extensive  research  and  development  efforts  over  the  past 
several  decades  has  led  to  the  proliferation  cf  monopulse  radar  tracking  systems  on 
both  fixed  and  airborne  platforms  [1].  Requirements  and  techniques  for  integration 
of  a  complete  monopulse  radar  system,  i.e.,  aperture,  receiver,  and  processor,  onto 
an  airborne  platform  is  not  unl'ke  conventional  radar  systems  except  for  space 
limitations.  Environmental  protection  of  the  radar  system  is  essential  and  is  typically 
in  the  form  of  a  protective  cover  called  a  radome.  The  radome  is  generally 
composed  of  lew-ioss  dielectric  materials  and  streamlined  so  as  not  to  interfere  with 
the  aerodynamic  performance.  Although  designed  to  be  "transparent"  to  the 
operating  range  of  frequencies,  an  meaning  plane  wave  while  passing  through  the 
radome  is  subject  to  amplitude  attenuation,  phase  from  distortion,  depolarization,  etc. 
Hence,  the  monopulse  radar  system  generally  processes  a  distorted 
electromagnetic  (EM)  wavefront,  which  results  in  tracking  errors  and  degraded 
performance. 
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1.1  Motivation 


This  dissertation  provides  a  detailed  analysis  of  such  degrading  effects  on  the 
accuracy’  and  tracking  performance  of  a  monopulse  radar  system.  This  research 
effort  is  also  aimed  at  predicting  monopulse  tracking  performance  degradations  in  the 
event  a  particular  "system"  component  has  to  be  modified  to  meet  changing 
requirements.  In  this  context  the  term  "system"  applies  to  an  integrated  system 
consisting  of  radome,  aperture,  and  monopulse  receiver/processor  components. 
Although  numerous  analyses  have  been  performed  on  individual  system  components 
m  the  past,  modeling  capabilities  and  analysis  of  the  overall  system  is  limited  and  are 
typically  component  specific  and  seldom  flexible  enough  to  accommodate  any  future 
modifications.  This  research  effort  combines  previously  developed  modeling 
techniques  and  analysis  procedures  to  model  and  analyze  the  entire  "front-to-back” 
performance  of  the  integrated  "system",  i.e.,  incident  EM  field  to  boresight  tracking 
error  curve.  For  given  changes  to  key  component  parameters  the  synergistic  effect 
of  the  overall  system  can  be  predicted.  Results  will  provide  a  robust  analysis 
procedure  and  valuable  modeling  tool  by  which  system  modification  effects  can  be 
predicted  and  analyzed. 

1.2  Problem  Statement 

The  problem  addressed  under  the  research  effort  involves  accurate 
characterization  of  radome  depolarization  and  phase  front  distortion  effects  on 
monopulse  receiver  tracking  performance.  Degraded  tracking  performance  is  best 
characterized  through  analysis  of  system  boresight  error  (BSE)  under  varying 
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conditions  presented  by  "system"  components.  Angular  "system"  BSE  is  defined  as 
the  difference  between  the  angle  indicated  by  the  monopulse  tracking  system  (angle 
for  which  the  monopulse  system  generates  NO  tracking/pointing  correction  signal) 
and  the  true  target/source  location.  Source/target  location  data  is  typically  provided 
to  aircraft/missile  guidance  and  control  systems  for  use  in  establishing  platform 
responses.  Any  "system"  BSE  results  in  an  inappropriate  platform  response,  i.e., 
incorrf.ct  guidance/control  decision,  potentially  limiting  mission  effectiveness  or 
causing  total  system  failure.  "System"  BSE  must  be  accurately  predicted  and 
minimized  to  ensure  cost  effective  and  reliable  systems  are  developed.  Previous 
research  efforts  have  successfully  used  BSE  prediction  as  a  metric  for 
establishing/validating  monopulse  tracking  perfoimance  [2,  3,  4.  5,  6],  Predicted 
BSE  estimates  within  ±  1  mRAD  of  measured  values  are  generally  regarded  as 
"excellent"  [5]  and  are  typically  obtained  by  propagation  models  using  computationally 
intense  surface  integration  techniques  [2,  4],  Ray-tracing  propagation  techniques 
were  considered  as  a  means  to  reduce  computational  intensity.  Comparative  studies 
of  available  propagation  techniques,  i.e.,  ray-trace  receive,  ray-trace  transmit,  plane- 
wave  spectrum,  and  surface  integration,  established  a  "consensus"  that  the  ray-trace 
receive  (also  called  fast  receiving  and  backward  ray-trace)  propagation  technique 
offers  the  best  compromise  between  accuracy  and  computational  intensity  for 
moderate  to  large-sized  radar/radome  systems  [3,  5,  6).  Efforts  using  the  ray-trace 
receive  propagation  lechmque  have  successfully  characterized  BSE  "response"  to 
varying  "system"  parameters,  ‘.e.,  radome  layer  thickness  and  electrical  properties, 
E-Field  frequency  and  polarization,  operating  temperature,  aperture  illumination. 
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etc.  [4,  6],  However,  these  efforts  experience  limited  success  in  predicting  actual  BSE 
values  which  compare  favorably  with  measured  data. 


1.3  Approach 

The  ray-trace  receive  propagation  technique  is  used  for  analysis  and  modeling  to 
achieve  the  goal  of  comput  ral  efficiency.  A  survey  of  previous  research  efforts 
which  used  this  propagation  technique  revealed/identified  several  areas  requiring 
"improvement"  for  obtaining  accurate  analysis  and  modeling  results.  The  effects  of 
ray  refraction  (deflection/spreading)  upon  propagating  through  the  tadome  were  not 
accounted  for  in  previous  efforts.  Therefore,  refractive  effects  on  BSE  prediction  had 
not  been  established  and  needed  to  be  accurately  accounted  for  in  the  current 
analysis  and  modeling  development  [3,  6]. 

The  "simple"  radome  structure  models  of  previous  efforts  generally  provide 
limited  modeling  capability,  accurately  modeling  single  layer  constant  thickness 
radomes  but  providing  minimal  flexibility  for  multi-layer  tapered  radome  designs 
These  "simple"  modeling  techniques  are  unacceptable  for  BSE  refractive  effect 
characterization,  requiring  improvement  for  detailed  analysis  and  modeling  purposes. 
To  accurately  account  for  refractive  effects  m  arbitrarily  shaped  single  and  multi-layer 
radome  designs  (consisting  of  both  constant  and  tapered  dielectric  layers),  closed- 
form  radome  surface  equations  are  used  in  conjunction  with  Snell’s  Law  of  Refraction 
to  establish  ray  propagation  paths  through  the  radome.  Multi-layer  radome  surface 
equations  are  developed  using  a  tangent  ogive  reference  surface,  allowing  virtually 
all  circularly  symmetric  radome  shapes  to  be  accurately  analyzed  and  modeled. 
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Reflected  E-Fields  within  the  radome  can  account  for  a  significant  portion  of  the 
efror  between  predicted  and  measured  BSE  values  [5.  6).  Incident  E-Field  reflection 
points  are  typically  calculated  via  some  form  of  "hit-or-miss”  technique.  Given  an 
established  reference  plane  which  is  perpendicular  to  the  direction  of  propagation, 
rays  are  traced  from  uniformly  spaced  grid  points  (typically  spaced  at  one-quarter  or 
one-half  wavelength  intervals)  through  the  radome  to  inner :  urface  reflection  points. 
Rays  are  reflected  and  checked  to  see  "if1  they  intersect  the  aperture  plane  [4],  Rays 
wnich  intersect  the  aperture  are  included  in  reflected  E-Field  calculations  at  the 
nearest  aperture  sample  pomt/element  location,  non-interserting  rays  are  excluded. 
Tne  current  analysis  and  modeling  technique  improves  on  this  "hit-cr-miss''  technique 
by  calculating  ray  reflection  points.  Fermat’s  principle  and  variational  calculus 
techniques  are  applied  to  multi-layer  surface  equations  resulting  in  a  system  of 
equations  which  provide  reflection  point  solutions. 

1.4  Research  Contributions 

The  following  list  is  a  summary  of  research  contributions  achieved  in  addressing 
the  problem  of  accurately  predicting  "system"  BSE.  The  contributions  are  a  result 
of  applying  the  approach  stated  in  Section  1.4  and  provide  "improved"  BSE  prediction 
capability.  Research  contributions  include: 

(1)  A  uniquely  defined  analysis  and  modeling  procedure  for  "systems"  including 
a  multi-layer  tapered  radome,  providing  closed-form  equations  for  each  radome 
surface  relative  to  a  tangent  ogive  reference  sui  face. 
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(2)  Established  "Ideal"  taper  function  criteria  for  use  with  (1),  resulting  in  a 
procedure  which  produces  a  "near"  optimum  taper  function  for  BSE  prediction. 

(3)  A  generalized  analysis  procedure  and  solution  technique  for  calculating 
radome  reflection  points  on  an  arbitrarily  shaped  reflecting  surface.  The  generalized 
solution  is  reduced  to  a  system  of  two  non-linear  transcendental  equations  with  two 
unknowns  for  the  case  of  a  tangent  ogive  reflecting  surface. 

(4)  Propagated  E-Field  expressions  and  analysis  results  using  Geometric  Optics 
and  surface  equations  of  (1)  accounting  for  total  refractive  effects  along  ray 
propagation  paths. 

(5)  A  validated  "system"  model  based  on  (I)  thru  (4)  which  predicts  BSE 
performance  using  "front-to-back"  propagation  characteristics.  Validation  is  based 
on  empirical,  published  experimental,  and  production  "system"  acceptance  test  data. 

(6)  Refractive  effect  characterization  in  principle  and  diagonal  scan  planes  using 
the  validated  model  of  (5). 
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II.  Background 


2.7  Introduction 

This  chapter  outlines  and  highlights  background  information  existing  on  the 
overall  analysis  and  modeling  effort  for  radomes,  antennas,  and  monopulse 
techniques.  Each  section  presents  an  overview  of  the  key  component  parameters 
relevant  to  this  research.  Additionally,  applicable  assumptions  are  identified  and 
justified. 

2.2  Radomes 

A  radome  is  a  protective  covering  which  is  generally  composed  of  low-loss 
dielectrics  and  designed  to  be  transparent  to  Electromagnetic  (EM)  waves  of  interest. 
The  actual  shape  of  the  radome  is  application  specific  and  despite  best  efforts  to  the 
contrary,  affects  both  the  aerodynamics  of  the  aircraft  and  electrical  performance  of 
the  radar.  Generally,  highly  streamlined  shapes  have  less  drag  and  can  better 
withstand  precipitation  damage;  however,  electrical  performance  usually 
deteriorates  [7],  Generally,  the  more  streamlined  a  radome  shape  becomes  the 
greater  the  incidence  angle  a  propagating  EM  wave  experiences.  Higher  EM  wave 
incidence  angles  typically  increase  reflections  on  the  radome’s  surface,  resulting  in  less 
energy  being  transmitted  through  the  radome  “all.  Additionally,  wavefront 
"depolarization''  occurs  as  the  EM  energy  propagates  through  the  radome  wall. 

Depolarization  occurs  when  an  EM  wave  propagates  through  boundaries 
exhibiting  discontinuous  electrical  and  magnetic  properties  namely,  permittivity  e, 
permeability  a,  and  conductivity  a.  The  effects  of  the  discontinuities  on  the 
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propagating  wave  can  be  expressed  in  terms  of  a  transmission  coefficient  T  and  a 
reflection  coefficient  F,  both  of  which  are  generally  complex  quantities.  For  uniform 
plane  wave  illumination  {a  wave  possessing  both  equiphase  and  equiamplitude  planar 
surfaces),  at  an  oblique  incidence  angle  on  the  interface  of  layered  media,  both  T  and 
T  are  functions  of  1)  the  constitutive  parameters  on  either  side  of  a  boundary, 
2)  the  direction  of  wave  travel  (incidence  angle),  and  3)  the  orientation  of  the  electric 
and  magnetic  field  components  (wave  polarization)  [8,  9]. 

Standard  techniques  for  analyzing  plane  wave  propagation  across  electrically 
discontinuous  boundaries  generally  begin  by  establishing  a  plane  of  incidence.  The 
plane  of  incidence  is  commonly  defined  as  the  plane  containing  both  a  unit  vector 
normal  to  the  reflecting  interface  (boundary)  and  a  vector  in  the  direction  of 
incidence  (propagation).  To  analyze  the  reflection  and  transmission  properties  of  a 
boundary  for  waves  incident  at  oblique  angles  with  arbitrary  polarization,  for 
convenience,  the  electric  or  magnetic  fields  are  decomposed  into  components  parallel 
and  perpendicular  to  the  plane  of  incidence  .  When  the  electric  field  is  parallel  to 
the  plane  of  incidence  it  is  commonly  referred  to  as  the  Transverse  Magnetic  (TM) 
case.  Likewise,  when  the  electric  field  is  perpendicular  to  the  plane  of  incidence  it 
is  commonly  referred  to  as  the  Trans-.erse  Electric  (TE)  case  (8). 

An  EM  wave  is  propagated  across  the  boundary  by  applying  the  appropriate 
parallel  or  perpendicular  transinission/reflection  coefficient,  i.e.,  Tg,  T,,  Fj,  oi  Tx, 
to  the  corresponding  component  of  the  electric  or  magnetic  field.  Eas  (1)  thru  (3) 
illustrate  this  process  where  the  superscript  i  represents  the  incident  field  components, 
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the  superscript  r  represents  the  reflected  field  components,  and  the  superscript  i 
represents  transmitted  field  components.  From  Eq  (2)  it  is  apparent  why  an  EM 
wave  often  experiences  depolarization  upon  propagating  through  a  homogeneous 
radome  layer.  Conditions  which  cause  the  parallel  and  perpendicular  transmission 
coefficients  to  differ,  either  in  amplitude  or  phase,  can  result  in  the  total  transmitted 
field  exhibiting  a  polarization  which  differs  from  the  incident  field. 


e[  =  rj[  ;  E'  =  rX 


(1) 


E[=Tjl  ;  E[  =T,£,| 


(2) 


(3) 


Figures  1  and  2  illustrate  typical  variations  which  occur  in  perpendicular  and 
parallel  transmission/reflection  coefficients.  Comparison  of  these  figures  shows  how 
the  magnitude  of  the  coefficients  varies  as  a  function  of  constitutive  parameters  and 
incidence  angle.  Incidence  angles  which  reduce  the  reflection  coefficient  to  zero  are 
referred  to  as  Brewster  angles  [8].  The  depolarization  effect  resulting  from  the 
parallel  and  perpendicular  transmission  coefficients  having  unequal  magnitudes  is 
graphically  illustrated  in  Figure  3.  Clearly,  a  wave  "depolariza'ion”  between  the 


incident  and  transmitted  electric  field  has  occurred. 
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Figure  1.  TE  Transmission  and  Reflection  Coefficients 
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Figure  2.  TM  Transmission  and  Reflection  Coefficients 
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Figure  3.  Graphical  Representation  of  Radome  Depolarization  Elfects 

Magnitude  and  phase  differences  between  parallel  and  perpendicular  reflection 
coefficients  will  result  in  depolarization  upon  reflection.  Any  depolarization  occurring 
along  a  direct  or  reflected  propagation  path  generally  introduces  a  polarization 
mismatch  between  the  incident  EM  field  and  radar  aperture,  which  reduces  aperture 
efficiency  and  degrades  the  overall  performance.  Since  the  degree  of  radome 
depolarization  depends  upon  the  incidence  angle  and  polarization  of  the  incident 
wave,  and  the  electrical  properties  of  the  radome  itself,  radome  transmission  and 
reflection  coefficients  emerged  as  important  parameters  for  analysis  and  modeling 
puiposes  in  the  present  research  effort. 

Previous  efforts  have  effectively  utilized  reflection  and  transmission  coefficients 
to  characterize  overall  radome  electrical  performance  [2.  4],  Equations  have 
emerged  and  elaborate  code  has  evolved  for  calculating  transmission  and  reflection 
coefficients.  This  research  effort  capitalized  on  past  efforts  by  using  equations 
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derived  by  Richmond  for  calculating  the  coefficients  [10].  Modifications  were  made 
to  the  Richmond  derivation  to  make  the  calculated  output  coefficients  equivalent  to 
the  more  elaborate  Periodic  Moment  Method  (PMM)  code  developed  by  the  Ohio 
State  University  (OSU).  The  PMM  code  is  capable  of  analyzing  periodic  structures 
imbedded  in  an  arbitrary  number  of  dielectric  slabs  of  finite  thickness  [11]. 
Although  compatibihty  with  the  more  sophisticated  code  was  not  necessarily  a  firm 
requirement,  1.  was  addressed  to  1)  assure  flexibility  and  wider  use  of  the  analysis  and 
modeling  technique,  and  2)  provide  an  alternate  source  for  comparison  and  validation 
of  particular  components  of  the  model.  In  both  cases,  Richmond’s  development  and 
the  OSU  PMM  code  were  developed  for  planar  dielectric  slabs  of  finite  thickness  and 
infinite  extent  in  both  height  and  width.  Their  use  in  approximating  T  and  T  for 
"locally  planar"  analysis  and  modeling  cases  has  been  very  effective  and  widely 
accepted  among  the  technical  community.  Generally,  "locally  planar"  conditions  exist 
at  any  given  point  cn  a  specified  surface  where  the  radii  of  curvature  of  the  surface 
at  that  point  are  large  compared  to  a  wavelength.  As  explained  and  justified  in  detail 
later,  this  research  effort  concentrated  on  system  component  designs  and  test  conditions 
by  which  locally  planar  approximations  could  be  enforced. 

The  primary  reference  surface  considered  in  the  research  effort  was  the  tangent  ogive 
for  several  reasons.  First,  many  practical  radomes  possess  at  ieast  one  surface 
(boundary  between  two  adjacent  layers)  which  is  either  truly  a  tangent  ogive  surface 
or  some  slight  variation  thereof.  Second,  a  tangent  ogive  surface  can  be  expressed 
by  a  closed-form  analytic  expression  making  it  a  useful  reference  surface  for 
generating  closed-form  expressions  for  non-ogive  surfaces.  Third,  the  tangent  ogive 
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radome  shape  has  been  analyzed  extensively.  As  such,  a  wealth  of  empirical  and 
measured  data  exists  for  comparison  and  validation  purposes.  Lastly,  a  majority  of 
practical  tangent  ogive  iadomes  have  physical  dimensions  that  are  large  compared 
to  their  design  wavelength.  Therefore,  the  radius  of  curvature  at  any  point  on  the 
radome's  surface  through  which  the  radar  aperture  must  "iook"  is  larj  e  in  comparison 
to  one  wavelength.  Tuir  condition  makes  locally  planar  techniques  valid  fcr  analyzing 
overall  system  performance  and  at  the  same  time  satisfies  the  locally  planar 
requirement  fer  calculating  transmission  and  reflection  coefficients. 

There  are  two  widely  used  methods  for  analyzing  radoine  transmission  and 
reflection  properties.  The  surface  integration  method  calculates  the  field  distribution 
due  to  the  antenna  on  the  inner  radome  surface  at  a  number  of  ciscrete  points  and 
the  far-field  response  as  an  integration  over  the  outer  surface  of  the  radome, 
accounting  for  the  transmissivity  effects  (12).  This  method  is  the  most  accurate  but 
is  also  much  more  time  consuming  to  implement  (computationally  intense).  The 
second  method  is  ray  tracing  (geometrical  optics)  where  transmitted  rays  are  assumed 
to  pass  directly  through  the  radome  wall  and  reflected  rays  are  assumed  to  originate 
at  the  point  of  incidence.  The  tiansmitted  electric  field  at  a  given  point  is  found  by 
considering  an  incident  ray  which  would  pass  through  ,he  point  with  the  slab  removed 
and  then  weighting  the  field  of  this  ray  with  the  appropriate  insertion  transmission 
coefficient  [4],  Although  the  ray  tracing  method  does  not  accurately  predict  sidelobe 
performance  outside  those  nearest  to  the  main  beam,  it  accurately  characterizes  the 
properties  of  the  main  beam  itself. 


14 


Figure  4  illustrates  an  application  of  the  ray  tracing  method  to  an  arbitrarily 
shaped  radome  enclosing  a  phased  array  antenna.  Given  a  specific  reference  plane 
location,  a  ray  normal  to  the  reference  plane  is  traced  through  the  radome  to  a 
specified  element  location  on  the  aperture.  At  the  ray  intersection  point  on  the  outer 
surface  of  the  radome,  a  surface  normal  vector  is  calculated  and  the  ray’s  angle  of 
incidence  determined.  Invoking  the  locally  planar  approximations,  transmission  and 
reflections  coefficients  are  calculated  and  used  to  weight  appropriate  field 
components  of  the  incident  EM  wave. 


Figure  4.  Ray  Tracing:  Main  Beam  and  Image  Lobe  Rays 


The  main  beam  (direct)  rays  are  determined  directly  by  weighting  the  appropriate 
parallel  and  perpendicular  field  components  of  the  incident  wave.  The  image  lobe 
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(reflected)  rays  identified  in  Figure  4  are  quite  different.  A  typical  radome  creates 
an  image  lobe(s)  anytime  significant  reflections  ex;st  within  the  radome.  As  seen  in 
the  figure,  an  image  lobe  ray  actually  experiences  the  effects  of  both  reflection  and 
transmission  as  it  propagates.  For  ray  tracing  analysis  image  lobe  rays  must  be 
weighted  by  both  a  transmission  and  reflection  coefficient. 

2.3  Antennas 

Many  characteristics/parameters  are  used  to  describe  the  performance  of  an 
antenna,  including:  radiation  pattern  (field  or  power  pattern),  input  impedance, 
polarization,  gam,  directivity,  half-power  beamwidth  (HPBW),  first-null  beamwidth, 
side  lobe  level  (SLL),  etc.  Although  all  these  characteristics  are  considered  important 
for  various  analysis  purposes,  not  all  were  considered  to  be  key  to  the  success  of  this 
research  effort.  Two  of  these  characteristics  were  considered  of  primary  importance  in 
the  overall  system  analysis  for  this  research  effort,  namely,  the  radiation  pattern  and  the 
polarization. 

The  radiation  pattern  of  an  antenna  is  a  graphical  representation  of  its  radiation 
properties  as  a  function  of  spacial  coordinates.  A  radiation  pattern  typically  displays 
the  variation  m  the  "field  strength"  (field  pattern)  or  the  "power  density"  (power 
pattern)  as  a  function  of  angle.  Henceforth,  any  reference  in  this  dissertation  to 
antenna  radiation  pattern,  01  simply  antenna  pattern,  is  implicitly  referring  to  the 
antenna  "field  pattern"  and  will  be  designated  as  f(8,<J>).  The  two  variables  0  and  0 
allow  for  the  field  pattern  to  vary  in  two  dimensions  and  f(0,t}>)  is  considered  to  be 
the  field  pattern  as  observed  in  the  antenna’s  fur-field. 
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Tne  primary  antenna  type  considered  was  a  planar  phased  array  with  corporate  feed. 
A  corporate  feed  structure  is  an  equal  path-length  branching  network.  The  electrical 
path  length  from  the  centra!  terminal  to  all  phase  shifters  is  the  same  and  the  feed 
provides  a  broadside  wavefront  for  the  phase  shifters  to  operate  on  [14].  This  choice 
was  believed  reasonable  considering  that  a  vast  majority  of  modern  airborne 
platforms  utilize  this  specific  technology  to  meet  their  aperture  needs.  Additionally, 
a  relatively  "simple"  method  exists  for  obtaining  monopulse  operation  from  a  phased 
array  antenna.  Details  of  this  method  and  the  radiation  patterns  required  for 
accurate  monopulse  performance  are  provided  m  Section  2.3,  Monopulse  Processing. 

A  phased  array  aperture  is  an  antenna  whose  main  beam  maximum  direction  or 
pattern  shape  is  primarily  cont-clled  by  the  relative  phase  of  the  element  excitation 
currents  on  the  array  [33].  The  spaciai  orientation  of  the  basic  radiating  elements 
which  form  the  array  structure  is  somewhat  arbitrary,  depending  on  the  desired 
overall  array  performance.  However,  the  basic  planar  phased  array  antenna  typically 
consists  of  equally  spaced  elements  on  a  rectangular  grid.  As  with  linear  arrays, 
grating  lobes  are  one  of  the  major  concerns  during  array  design.  Grating  lobe .  are 
additional  maxima  which  appear  in  the  field  pattern  w  ith  intensities  nearly  equal  to 
the  intensity  of  the  mam  beam.  The  effects  of  radiative  grating  lobes  are 
minimized/eliminated  by  ensuring  the  distance  between  any  two  adjacent  elements 
within  the  array  is  less  than  one  wavelength  (A.).  Generaily,  the  spacing  between  any 
two  adjacent  array  elements  is  maintained  at  one-half  a  wavelength  ( ?./2). 

In  analyzing  the  field  pattern  of  array  antennas  the  principle  of  "pattern 
multiplication"  is  often  used.  As  shown  in  Eq  (4),  pjtietn  multiplication  allows  the 
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array  field  pattern  f(0,4>)  to  be  factored  into  the  product  of  the  element  pattern, 
fe(e,4>),  and  the  array  factor,  AF(6,<J)).  The  element  pattern,  ft(6,<j>),  is  the  field 
pattern  of  a  typical  element  within  the  array.  The  complex  C„  term  in  Eq  (4) 
represents  mutual  coupling  effects  for  the  nth  array  element.  Although  the  value 
of  C.  is  dependent  on  element  location,  mutual  coupling  effects  are  generally 
assumed  to  be  identical  far  all  elements  within  the  array.  This  is  particularly  true  for 
densely  populated  arrays  which  use  amplitude  tapers  for  pattern  control. 

N 

/(M)  =/t(e.4.)/lF(e,4>)  =  /£(0,<t>)]T  C^expC/pr.r,')  (4) 

/>•! 


The  array  factor  AF(0,<}))  is  the  field  pattern  of  the  array  when  the  elements  are 
replaced  by  point  sources  (isotropic  radiators)  excited  by  the  same  current  amplitudes 
and  phases  as  the  original  elements  [13],  The  vector  r  j  is  a  position  vector  from  the 
coordinate  origin  to  the  element  position  and  r  is  a  unit  spherical  radius  vector  for 
the  same  coordinate  system.  The  /,  term  in  Eq  (4)  represents  the  complex  excitation 
current  of  the  n,k  array  clement  and  may  be  expressed  as  given  in  Eq  (5),  where  A, 
and  a„  are  real  quantities  representing  the  amplitude  and  phase  of  the  excitation 
current,  respectively.  Assuming  the  array  of  elements  and  feed  network  form  a 
passive  reciprocal  structure,  the  phased  array  will  exhibit  identical  field  patterns  for 


bosh  transmission  and  reception  [14].  This  assumption  may  not  be  valid  if  the  feed 
network  contains  any  non-linear  or  directive  components  such  3s  ferrite  phase  shifters 
or  amplifiers.  For  the  research  effort  a  passive  reciprocal  array  structure  was  assumed 
allowing  for  reciprocity  to  be  invoked  whenever  necessary. 

The  polarization  of  the  antenna  array  is  the  second  characteristic  of  importance 
to  the  research  effort.  Mutual  coupling  between  array  elements,  seen  earlier  to  cause 
the  element  field  patterns  to  vary  from  element  to  element,  can  cause  a 
depolarization  effect,  i.e.,  the  phased  array  may  exhibit  a  polarization  unlike  the 
polarization  of  the  individual  elements  [13].  For  maximum  coupling  of  energy 
between  an  EM  wave  and  an  antenna  to  occur  the  polarization  of  the  wave  and  the 
antenna  must  be  matched. 

In  practice  a  true  polarization  match  between  an  incident  EM  wave  and  receiving 
antenna  is  seldom  achieved.  Design  and  manufacturing  limitations  both  contribute 
to  degraded  antenna  performance,  resulting  in  the  antenna’s  polarization  state  being 
less  than  "pure".  Also,  the  medium  through  which  the  EM  wave  must  propagate  to 
reach  the  antenna  aperture  will  likely  degrade  the  purity  of  the  incident  wave’s 
polarization.  This  latter  case  was  discussed  in  detail  in  Section  2.1 ,  Radomes,  where 
sigmficani  polarization  variations  were  seen  to  occur  as  a  result  of  radoir.e  design  and 
materials.  The  electric  field  transmitted  by  the  radome  was  seen  to  vary  significantly 
as  the  angle  of  incidence  of  the  incident  electric  field  was  varied. 

To  account  for  the  effects  of  polarization  mismatch  between  antennas  and 
incident  waves,  the  concept  of  a  Polarization  Loss  Factor  (PLF)  is  introduced.  The 
term  loss  is  used  since  a  polarization  mismatch  results  m  a  condition  where  the 
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amount  of  power  extracted  from  an  incoming  wave  by  an  antenna  will  not  be 
maximum  [13].  The  PLF  between  an  antenna  and  incident  electric  field  may  be 
defined  as  given  by  Eq  (6)  where  pw  is  a  complex  unit  vector  defining  the  polarization 
of  the  incident  wave,  p,  is  a  complex  unit  vector  defining  the  polarization  of  the 
antenna,  and  is  the  angle  between  the  two  unit  vectors  [15]. 

PLF  =  |p„*  pj2  =  Icos^l2  (6) 

The  PLF  can  vary  between  zero  and  one  with  a  PLF  =  1  representing  an 
ideally  polarized  (perfect  match)  condition  and  a  PLF  =  0  representing  an 
orthogonal  (total  mismatch)  condition.  These  two  extreme  cases  are  referred  to  as 
"co-polarized"  (PLF  =  1)  and  "cross-polarized"  (PLF  =  0)  conditions  [13].  Neither 
of  these  two  extremes  ate  likely  to  occur  in  practice.  However,  the  terms  co- 
polarized  and  cross-polarized  are  generally  used  co  describe  an  antenna’s  response, 
as  done  throughout  this  dissertation.  A  polarization  reference  direction  is  established 
to  coincide  with  the  polarization  direction  of  the  aperture  elements.  Fields  are  then 
resolved  into  components  which  are  perpendicular  (cross-polarized)  and  parallel  (co¬ 
polarized)  to  this  reference  direction. 

This  convention  was  used  in  Section  2.1  where  the  reference  polarization  direction 
was  established  as  the  antenna's  polarization  direction.  An  electric  field  radiated  by 
the  antenna  was  then  propagated  through  the  radome  resulting  in  a  transmitted 
electric  field  polarization  vector  which  had  "rotated".  Applying  reciprocity,  the  overall 
antenna/radcme  system  would  now  respond  to  EM  waves  with  polarization  states 
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different  from  the  aperture.  The  depolarization  effect  of  the  radome  alters  the 
response  of  the  original  system,  a  system  designed  by  assuming  ideal  transmission  and 
reflection  properties  for  the  radome. 

2.4  Monopulse  Processing 

To  actually  define/describe  what  constitutes  a  monopulse  receiver/processor 
system  is  difficult  and  varies  from  author  to  author.  However,  the  concept  that 
"monopulse"  implies  determining  the  relative  bearing  of  a  radiation  source  by 
analyzing  the  received  characteristics  of  a  “single  pulse"  of  energy  from  that  source 
is  consistently  seen  throughout  literature.  This  "single  pulse"  tracking/locating  is 
typically  accomplished  by  one  of  three  basic  monopulse  techniques.  These  three 
techniques  are  characterized  by  the  manner  in  which  angle  information  is  extracted 
from  the  received  signal  and  include  amplitude-comparison,  phase-comparison,  and 
amphtude-phase-comparison  (a  combination  of  the  two  previous  techniques)  [16). 
System  processing  requirements  for  implementing  each  technique  can  vary 
significantly  as  will  be  seen  throughout  the  course  of  this  discussion. 

Amplitude  Comparison  Monopulse  (ACMP)  processing  is  perhaps  the  simplest 
form  of  "single  pulse"  processing.  A  basic  ACMP  system  typically  consists  of  an 
antenna  system,  receiver/processor,  and  a  display,  and  is  generally  used  as  a  baseline 
reference  for  analysis  and  discussion  purposes  In  the  baseline  system,  the  antenna 
system  employed  typically  consists  of  a  pair  of  radiating  elements  having  identical 
field  patterns,  say  f(6).  These  radiating  elements  are  physically  located  in  such  a 
manner  that  their  field  pattern  maximums  occur  at  angles  of  ±  9^  with  respect  to  the 
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antenna  system’s  boresight  axis  (0  =  0‘),  0„  is  referred  to  as  the  "squint"  angle.  This 
same  field  pattern  relation  may  also  be  achieved  by  a  phased  array  antenna  via 
electronic  beam  steering  techniques.  Henceforth,  these  two  field  patterns  will  be 
referred  to  as  f,(0)  for  the  0  -  0^  pattern  and  f,(0)  for  the  0  +  0^  pattern.  The 
voltages  corresponding  to  each  antenna  output  will  be  referred  to  as  v,  and  v2  for  the 
fi(0)  and  f2(0)  field  patterns,  respectively.  These  voltages  may  be  expressed  as: 


=  Ale,e'f(Q  -0^)  =^,ejP'/1(0) 

(7) 

=  A2eJ*--f((!  +  e^)  =A2e,hf2(0) 

(8) 

where  Ah  and  A j,  /3j,  are  proportional  to  the  amplitude  and  phases  associated 
with  an  impinging  EM  wavefront  on  antenna  #1  and  antenna  #2,  respectively.  These 
equations  indicate  that  the  antenna  voltages  developed  are  in  general  complex 
quantities.  At  this  point  no  consideration  has  been  given  to  the  EM  wave 
polarization  match/mismatch  conditions.  Rather,  it  has  been  assumed  here  that  the 
polarization  of  the  impinging  wave  identically  matches  the  polarization  of  the 
antennas  being  used,  a  condition  not  existent  in  a  practical  situation.  This  matched 
polarization  assumption  is  carried  throughout  the  remainder  of  the  monopulse 
discussion.  Polarization  effects  identified  earlier  in  Section  2.2  as  both  cross-polarized 
and  co-polarized  antenna  field  pattern  effects  will  be  incorporated  into  the  equations 
later  since  the  study  and  analysis  of  these  effects  was  the  primary  concern  of  the 


research  effort. 


A  basic  ACMP  system  typically  accomplishes  angle  tracking  by  forming  a  ratio 
consisting  of  the  "difference"  (del  or  A)  between  the  antenna  voltages  to  the 
"sum"  (r)  of  the  antenna  voltages.  Since  phased  array  antenna  systems  can 
independently  form  any  desired  set  of  radiation  patterns,  within  practical  limits,  the 
"sum"  and  "difference"  patterns  are  usually  formed  directly  rather  than  by  addition 
and  subtraction  of  outputs  of  the  individual  field  patterns  f,(6)  and  f2(0)  [1],  As  such, 
outputs  of  the  individual  beams  may  not  be  available  for  analysis  anywhere  in  the 
system.  The  ACMP  ratio  formed  is  as  expressed  in  Eq  (9). 

A  _  del  _  v,  -  ^ 

E  sum  v,  +  v2 

Antenna  vohages  v,  and  v,  are  generally  complex  quantities  resulting  m  the 
monopulse  processing  ratio  of  Eq  (9)  also  being  complex.  However,  a  couple  of 
simplifying  assumptions  are  often  made  for  analyzing  the  basic  ACMP  system.  These 
assumptions  include  the  antennas  having  collocated  phase  centers  and  the  source  of 
radiation  being  located  in  the  far-field  (plane  wave  incidence)  These  assumptions 
result  in  a  condition  where  p,  equals  02  causing  the  monopulse  processing  ratio  to  be 
purely  real  (in  this  case  the  relative  phase  difference  between  the  sum  and  difference 
is  zero  degrees). 

Because  of  its  role  m  indicating  error,  i.e.,  a  source  located  left/right  of  boresight 
results  in  a  ratio  with  relationships  indicating  both  magnitude  (how  far  off  boresight) 


1  - 


(9) 


1  + 
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and  phase  (which  direction  off  ooresight),  the  ratio  signal  is  often  referred  to  as  the 
"error  signal",  as  will  be  done  throughout  the  remainder  of  the  dissertation.  A  steep 
slope  on  the  error  signal  curve  represents  a  high  degree  of  sensitivity  to  changes  in 
source  location  whereas  a  shallow  slope  represents  a  low  degree  of  sensitivity.  The 
error  signal  generated  by  the  ACMP  receiver/processor  is  generally  used  for  either 
indicating  the  location  of  a  radiating  source  via  some  form  of  visual  display  or  to 
correct  the  pointing  direction  of  the  boresight  axis.  The  actual  use  of  the  error  signal 
was  not  the  focus  of  the  research  effort;  the  error  signal  itself  is  analyzed  under  varying 
conditions  presented  by  the  antenna  and  radome  components. 

Phase  Comparison  Monopulse  (PCMP)  is  similar  to  ACMP  in  that  both  require 
the  same  basic  components  (antenna  system,  receiver/processor,  display).  However, 
the  design  and  functions  of  these  components,  especially  the  antenna  and 
recewer/processor,  may  vary  significantly.  The  basic  PCMP  system  utilizes  a  pair  of 
antenna  elements  with  identical  field  patterns  f(9),  similar  to  the  ACMP  system.  A 
major  difference  in  the  two  monopulse  systems  is  that  the  antenna  elements  in  the 
PCMP  system  are  not  squinted  but  physically  displaced  by  distance  d.  Thus 
separation  effectively  separates  the  antenna  phase  centers.  The  antenna  phase  center 
is  a  point,  located  on  the  antenna  or  in  space,  whereby  fields  radiated  by  the  antenna 
and  referenced  to  this  point  are  spherical  waves  with  ideal  spherical  wavefronts  or 
equiphase  surfaces  [15],  Under  far-fteld  conditions  parallel  ray  approximations  are 
used  to  estimate  antenna  output  voltages.  The  expressions  given  by  Eqs  (10) 
and  (11)  represent  typical  source  to  observer  approximations  where  R  is  the  distance 
from  the  antenna  coordinate  axis  and  r,/r;  are  the  distances  from  the  antenna 
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element;,  to  the  source.  These  approximations  for  the  amplitude  and  phase  remain 
valid  provided  the  condition  that  d  <  <  R  is  maintained. 


AMP  :  r,  =R  =  r2  (10) 

PHASE:  r,  =  R  +  ||jsin0  ;  r2  =  fl  -  jsin6  (H) 

Given  these  approximations  and  considering  only  single  source  illumination,  the 
voltages  present  at  each  of  the  antennas’  terminals  may  be  expressed  as  given  by 
Eq  (12)  and  (13)  where  A  and  f3  represent  the  amplitude  and  phase  of  the  source, 
respectively.  The  approximate  forms  of  v,  and  v2  in  Eqs  (i2)  and  (13)  clearly  show 
that  while  the  magnitudes  of  v,  and  v,  vary  with  source  location,  they  are  equivalent 
for  all  values  of  0.  Therefore,  no  angle  discriminating  information  is  contained  in  the 
voltage  magnitudes.  However,  the  equations  do  indicate  a  difference  between  the 
phase  terms  in  the  approximate  voltage  forms  of  v,  and  v2.  Each  phase  term  is  seen 
to  vary  with  the  antenna  separation  distance  d,  the  wave  number  K  =  2 n/>.,  and  the 
source  location  0.  Because  angle  discriminating  information  is  contained  only  in  the 
phase  difference  between  the  two  antenna  voltages,  the  term  phase-comparison  is 
applied.  Using  Eqs  (12)  and  (13)  the  phase  difference  A0  between  voltages  v,  and 
v2  may  be  expressed  as  given  by  Eq  (14). 

.  -/(Kr,  •  P)  ..Kit 

v,  =  /s (0,)  — - =/(0)A€-'lw>“‘8*P»  - -  (li) 

r\  R 
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v2  =  /2(02)  — - =  f(6)A  e'/umte-H  £ - 

r,  f? 


(13) 


A4>  =  4>,  -  4>2  =  JCdsind  =  sinO  (14) 

The  A <f>  expression  of  Eq  (14)  characteristically  represents  an  interferometric 
variation  of  Aa>  as  0  is  varied  and  exhibu  avorable  characteristics  for  processing 
source  tracking/locating  data,  specifically,  a  linear  response  with  a  positive  slope 
passing  through  zero  at  0  =  0'  the  source  location.  This  condition  is  consistent 
with  results  shown  previously  for  the  ACMP  processing  system  with  one  major 
exception.  The  PCMP  system  exhibits  several  positive  slope  lines  passing  through 
zero  m  the  presence  of  single  source  illumination  (usually  referred  to  as  ambiguities). 
Since  the  location  of  these  ambiguities  is  dependent  on  the  antenna  separation 
distance  d  and  wavelength  X  of  the  arriving  signal,  care  must  be  taken  in  the  PCMP 
processing  to  ensure  erroneous  angle  data  is  not  used  Similar  processing  techniques 
are  often  used  for  both  the  ACMP  and  PCMP  type  systems.  In  reality  even  practical 
ACMP  systems  exhibit  these  interferometric  effects.  It  is  practically  impossible  to 
perfectly  align/generate  squinted  beam  patterns  with  coincident  phase  centers,  a  fact 
neglected  in  the  previous  discussion  of  the  ACMP  system. 

As  mentioned  in  Section  2.2,  Antennas,  the  primary  antenna  configuration 
considered  for  analysis  and  modeling  throughout  the  research  is  the  phased  array 
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aperture  with  a  corporate  feed.  This  choice  was  deemed  reasonable  considering  that 
a  vast  majority  of  modern  airborne  platforms  use  phased  array  technology  to  meet 
their  aperture  requirements.  The  research  effort  concentrated  on  a  "simple"  method 
for  obtaining  monopulse  operation.  Figure  5  depicts  how  an  array  antenna  can  be 
divided  into  symmetrical  quadrants/sections.  The  outputs  of  all  elements  within  each 
quadrant/section  are  summed  to  produce  sum  signals/beams,  one  per 
quadrant/section,  which  are  then  in  turn  combined  to  form  appropriate  antenna 
system  sum  and  difference  voltages  [!}. 


Figure  5.  Phased  Array  Monopulse  Processing 

The  component  blocks  identified  by  the  large  4>  label  in  Figure  5  represent  the 
complex  weightings  which  drive  the  individual  array  elements.  The  "corporate"  feed 
structure  feeding  the  weights  is  designed  such  that  the  path  length  from  the  central 
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terminal  (one  per  section)  to  the  weighting  functions  is  the  same  and  thus  the  feed 
provides  a  broadside  wavefront  for  the  weighting  functions  to  operate  on  [14]. 
Complex  weighting  functions  allow  for  both  amplitude  and  phase  variations  to  be 
applied  across  the  array  aperture,  i.e.,  from  element  to  element.  By  providing 
amplitude  variation  across  the  aperture  it  is  possible  to  obtain  optimum  radiation 
patterns.  For  the  research  effort,  uniform,  cosine  (radial),  and  Taylor  amplitude 
distributions  were  considered.  Uniform  and  cosine  amplitude  tapers  were  considered 
primarily  for  comparison  with  existing  measuied  data  and  the  Taylor  amplitude  taper 
considered  since  it  generally  provides  the  highest  antenna  gain  for  a  given  sidelobe 
level.  It  has  also  been  shown  that  for  monopulse  processing  systems  the  Taylor 
weighting  provides  the  greatest  monopulse  sensitivity  slope  factor  and  efficiency  for 
a  given  sidelobe  ratio  [17],  Phase  control  is  typically  used  in  a  phased  array  to 
perform  beam  steering.  However,  electronic  beam  steering  via  element  phase  control 
may  introduce  some  asymmetrical  antenna  properties.  For  completeness  and  overall 
analysis  and  modeling  robustness,  electronic  scanning  of  the  phased  array  aperture  was 
considered  within  the  scope  of  the  research  effort.  Throughout  the  analysis  and 
modeling  validation  phase,  comparisons  between  mechanically  scanned  and 
electronically  scanned  results  are  made  when  appropriate. 
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III.  Multi-Layer  Radome 


3.1  Introduction 

This  section  provides  the  motivation  and  development  of  the  multi-layer  radome 
equations  used  in  the  analysis  and  modeling  of  the  overall  radome,  amenna,  and 
monopulse  receiver  system.  The  approach  taken  in  regard  to  the  multi-layer  radome 
orientation  and  structure  was  driven  by  several  factors.  First,  a  method  was  required 
to  analyze  and  model  tapered  radomes  since  such  structures  minimize  internal 
reflections  and  variation  of  insertion  phase  delay  as  a  function  of  antenna  look  angle. 
Second,  because  Geometric  Optics  (GO)  is  employed,  a  surface  normal  vector  is 
required  for  a  given  point  within  the  radome.  The  task  of  determining  a  surface 
normal  vector  is  simplified  if  the  surface  is  mathematically  describable  by  a  function 
which  is  differentiable  with  respect  to  all  variables  of  the  specified  coordinate  system. 
In  this  case,  the  gradient  of  the  function  provides  the  surface  normal  vector.  Lastly, 
the  prohte.ation  of  "tangent  ogive"  shaped  radomes  in  missile  and  aircraft  applications 
make  this  a  realistic  shape  for  analysis  and  modeling  purposes.  The  quotation  marks 
are  added  to  emphasize  that  a  radome  may  exhibit  a  shape  whose  overall  appearance 
is  consistent  with  a  tangent  ogive’s  appearance  without  actually  containing  an}  "true" 
tangent  ogive  surfaces.  A  tangent  ogive  surface  is  considered  since  it  will  be  used  as 
a  "reference"  surface  for  the  remainder  of  the  development.  The  analysis  proceeds 
by  considering  a  single  layer  radome  consisting  of  a  tangent  ogive  reference  surface 
separated  from  a  second  surface  by  a  specified  amount.  Expressions  derived  for  the 
single  layer  radome  are  then  extended  to  the  multi-layer  radome. 
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The  tangenl  ogive  radome  is  perhaps  the  most  commonly  used  and  modeled 
radome  shape.  It  can  be  thought  of  as  a  compromise  between  the  desirable 
aerodynamic  performance  of  a  cone  and  the  desirable  electrical  performance  of  a 
hemisphere  [6).  The  tangent  ogive,  illustrated  in  Figure  6,  is  a  surface  of  revolution 
generated  by  rotating  a  circular  arc  about  a  cord  (every  plane  tangent  to  its  base  is 
perpendicular  to  ihe  base  plane).  The  surface  of  a  tangent  ogive  is  completely 
described  by  its  generating  radius  R  and  base  diameter  W.  For  a  tangent  ogive  which 
has  its  base  plane  in  the  x-y  plane  such  that  its  tip  lies  on  the  z-axis,  as  illustrated  in 
Figure  6  for  a  right-hand  coordinate  system,  any  point  which  lies  on  the  ogive  surface 
satisfies  Eq  (15)  where  0  <  z  <  L.  The  ogive  length  L  along  the  z-axis  is  given 
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by  Eq  (16).  The  "fineness  ratio"  of  a  radome  is  defined  as  the  ratio  LAV.  For  typical 
radomes  the  fineness  ratio  is  greater  than  one  and  Eq  (16)  can  be  used  to  show  that 
R  >  W. 


g(x,y .z)  =  JR2  -  z2  -  JxT*  y1  +  5K/2  -  R  =  0 


(15) 


L  =  Jr?  -  (R  -  wj 2f 

In  its  rectangular  form,  the  right-hand  side  of  Eq  (15)  can  be  seen  to  represent 
the  radii  of  curvature  for  a  family  of  circles  which  lie  in  planes  parallel  to  the  s-y 
plant  and  which  vary  in  length  along  the  z-axis.  Therefore,  to  generate  a  second 
surface  located  on  either  the  inside  or  outside  of  the  reference  tangent  cgive  surface 
an  additional  radius  factor  can  be  subtracted  or  added  to  the  right-hand  side  of  the 
equation,  subti  action  generates  an  inner  surface  whereas  addition  generates  an  outer 
surface.  It  is  this  additional  radius  factor  which  will  be  established  in  the 
development  of  the  single  layer  radome  equations. 

3.2  Single  Layer  Development 

For  a  single  layer  radome  the  geometry  of  Figure  7  is  used  for  developing  the 
surface  equations.  It  is  assumed  that  the  desired  taper  function  t(0)  for  the  radome 
layer  is  known,  i.e,  t(0)  can  be  specified  based  on  the  desired  ladome  electrical 
performance  as  a  function  of  antenna  lock  angle  0.  Tr.c  reference  ogive  surface  is 


specified  by  the  parameters  R,  and  W,  with  R,  >  W,  yielding  fineness  ratios  greater 
than  one  and  satisfying  the  spnerica’,  coordinate  form  of  the  tangent  ogive  equation 
as  given  by  Eq  (17).  Since  the  surface  ts  symmetric  m  4!  the  analyst:  is  simplified  by 
considering  the  4*  =  0'  plane  (x-r.  plane)  for  equation  development,  the  Jesuits  of 
which  are  readily  extendable  for  any  arbitrary  4>"piar.e. 


>;.y 


Figure  7.  Single  Layer  Radome  Geometry 


As  shown  in  Figure  7  tz(0)  represents  the  radius  parameter  required  to  be 
added/subiracted  from  the  reference  surface  in  order  to  generate  the  outer/inner 
surface,  respectively.  The  goal  is  to  determine  the  function  tz(0)  required  to  modify 
g,t((r,0)  such  that  the  desired  second  surface  is  generated.  For  a  specified  value  of  0 
the  law  of  cosines  may  be  applied  to  the  geometry  of  Figure  7  to  solve  for  R„(0).  By 
computing  the  roots  of  a  quadratic  equation  while  enforcing  the  conditions  R,  >  W„ 
and  0  <  9  <  jr/2,  Ro(0)  can  be  derived  as  summarized  in  Eqs  (18)  thru  (22). 


Applying  Law  of  Cosines 
R)  »  «>)  *  £  -2KrRa(6)  cos(tt/2  -  0) 


(18) 
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With  R„(8)  determined  the  process  of  determining  tz(0)  continues  by  calculating 
the  coordinates  identified  as  x'(0)  and  z'(0)  in  Figure  7.  Once  calculated  per 
Eq  (23),  these  coordinates  are  substituted  into  the  expression  given  by  Eq  (24)  to 
obtain  the  final  form  expression  for  tz(0).  The  final  tz(0)  expression  is  then  added 
to  the  reference  surface  as  shown  in  Eq  (25)  to  obtain  g(r,0).  Eq  (25)  represents  the 
spherical  coordinate  expression  for  the  radome’s  second  surface  and  possesses  the 
characteristic  of  being  differentiable  with  respect  to  all  coordinate  variables  over  the 
desired  range  of  0  values.  It  is  evident  from  Eq  (25)  that  g(r,0)  identically  equals 
zero  when  r  =  A(0)  for  the  desired  range  of  r  and  0  values.  Therefore,  the 
equation  of  the  second  surface  can  be  simplified  as  shown  in  Eqs  (26)  thru  (28). 

x'(6)  =  Jr2  -  zW  *  WI2  -  R, 

(23) 

z'(8)  »  [*,,<6)  +  r(0)]cos(8) 


rz(6)  =  /{R„(0)  ♦  i(0) ]2  -  z'(0)2  -  x'(Q) 

rz(0)  =  i4(e)sin(0)  -  Jr2  -  47(0)cos2(0)  -  V/J2  +  R, 
where:  4(0)  =  Ro(0)  +  r(0) 


«(r,0)  =  gttf(.rfi)  +  M(0)  =  0 


g(rfi)  = 


|  Jr)  -  rW(0)  -  //?,2  -  A2(6)cos5(6) 
I  +  [i4(0)  -  r]sin(6) 


(25) 


=  0 


g(r,8)  -  /t(8)  -  r  =  R0( 6)  *  f(6)  -  r  =  0 


S(r,0)  =  ^  -  A?cosJ(6)  -  ATrsin(0)  +  r(0)  -  r  =  0 
IV 

for  Kr-R ,  _ 


£(r,0)  =  -  (R,  -_ivr/2)W(0)  -  (/?,  -  IK,/2jsm(0)  -  ((0)  -  r  =  0 

for  [w,!2  +  f(rc/2)j  s  r  s  |£.r  +  r(0)J  and  0  s0  s  n/2 


(26) 


(27) 


(28) 
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3.3  Multi-Layer  Development 

Extending  the  single  layer  derivation  process  to  the  multi-layer  radome  is  the  next 

task.  The  radome  geometry  and  coordinate  orientation  for  an  m-layer  radome  are 

shown  in  Figure  8.  For  the  multi-layer  development,  let: 

m  —  Total  number  cf  radome  layers,  including  any  resulting  from  the 
introduction  of  an  artificial  reference  surface. 

k  =  Total  number  of  radome  layers  inside  the  reference 
surface,  0  <  k  <  m. 

i  =  Surface  number  under  consideration,  1  <  /  <  m+1. 

t,(0)  =  Taper  of  the  i®  layer  relative  to  the  i®  surface,  specified  in 
the  f  direction. 

ts,(8)  =  Taper  "sum"  (equivalent  taper)  of  the  i®  layer  relative  to  the 
reference  surface. 


x.y 


Figure  8.  Multi-Layer  Radome  Geom-try 
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For  the  m-layer  radome  there  are  a  total  of  m  +  1  surfaces  for  which  equations 
must  be  generated.  As  with  the  single  layer  radome  development,  the  surface 
labeled  r  in  Figure  8  is  a  '  true''  tangent  ogive  surface  used  as  a  reference.  The 
reference  surface  r  need  NOT  be  an  actual  surface  of  a  radome  layer.  An  "artificiar 
reference  surface  may  be  introduced  anywhere  between  the  radome’s  inner  and  outer 
surfaces  provided  the  taper  functions  ts,(0)  and  t,(6)  can  be  specified  for  each 
radome  surface. 

Tne  following  derivation  and  analysis  process  applies  to  any  arbitrarily  shaped 
radome  which  when  oriented  as  shown  in  Figure  8,  1)  has  all  layer  surfaces 
independent  of  (j>,  i.e.,  symmetric  about  the  z-axis  in  any  given  (J>  plane,  2)  has  at  least 
one  surface  which  satisfies  the  tangent  ogive  equation  g,eI(r,0)  or  at  least  a  cross- 
sectional  area  such  that  a  tangent  ogive  curve  lies  entirely  between  the  radome’s 
inner  and  outer  surface,  thereby  allowing  for  an  artificial  reference  surface  to  be 
introduced,  and  3)  is  constructed  such  that  a  taper  function  t,(0)  can  be 
deterxnined/descnbed  for  each  layer,  including  any  artificial  layers  introduced  by  an 
artificial  reference  surface. 

Considering  the  surface  of  the  multi-layer  radome,  the  first  task  is  to  derive  an 
expression  for  the  total  taper  sum  ts,(0).  As  previously  defined,  ts,(0)  represents  the 
"sum"  or  equivalent  taper  of  the  i“  radome  surface  relative  to  the  reference  surface 
Once  calculated  ts,(0)  can  be  substituted  into  the  A(0)  expression  of  Eq  (26)  to 
obtain  g,(r,0).  G,(r,0;  represents  the  surtace  equation  describing  the  shape  of  the «“ 
surface.  For  a  given  0  value  the  taper  function  for  the  inner  most  radome  surface 
j  ts(f0)  can  be  expressed  as  shown  in  Eq  (29). 

i 

l 


37 


M,(0)  =  -E'/0) 


(29) 


The  taper  value  of  the  inner  most  surface  for  a  specified  0  value  is  simply  the 
negative  of  the  sum  of  all  t,(6)  values  prior  to  the  reference  surface.  Therefore,  the 
taper  value  for  the  i*  surface  is  obtained  by  adding  the  sum  of  all  t,(0)  values  prior 
to  the  i*  surface  to  ts,(0)  as  shown  in  Eq  (30).  Given  Eq  (30)  the  final  mathematical 
expression  for  ts(0)  is  obtained  by  the  development  process  provided  in  Eqs  (31) 
thru  (35),  resulting  in  the  final  expression  as  given  by  Eq  (36). 


0-1)  0-1)  i 

»,(0)  =  E  f/0> +  «i(0)  =  E  r>(°)  -  E*i(0) 

y-i  j-i  i-i 

i»i 

for  1  s  j  s  m*  1 


Expanding  The  Finite  Sum  Expression : 

-  [«,(0)  ♦  0(0)  ♦  ...  ♦  t,.,(0>] 

«,(0>  * 

-  [t,(0)  *  i2(6)  ♦  ...  *  6,(6)] 

For  k  >  (i-1) 

»,(«)  =  -{',0)  +  t.„(e) "  •  •  -  0(9)} 


(30) 


(31) 


(32) 
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For  k  <  (i- 1) 


(33) 


W((0)  =  +{W0>  -  <t.2(0)  +  ...  +  r,_,(0)} 


Fori  =  (i-1) 

»,(0)  =  0 


(34) 


Define  Signum(i-k-\)  as  Follows- 


sgn(i-k- 1) 


-1,  i  <  (i+1) 
0,  i  =  (i+1) 
+  1,  j  >  (i+l) 


(35) 


Letung  q  -  max(i-l,i)  and  s  =  rain  (i-1,  A)  +  1  then, 

»,( 0)  =  sgn(i-k-l)  Y,  f  (0)  (30) 

/or  1  <  i  s  m*I 

The  ts,(9)  function  of  Eq  (36)  represents  the  total  taper  value  to  be  added  to  the 
reference  surface  for  any  given  value  of  6.  It  can  be  substituted  directly  into  the 
single  layer  radome  Afu)  exDression  of  Eq  (26)  to  provide  the  final  form  expression 
for  g,(r,6).  The  final  expressions  obtained  are  as  given  by  Eqs  (37)  and  (38).  As  with 
B.-tifo0)  ind  gtr.6),  g,(r,0)  is  independent  of  (J)  and  is  differentiable  with  respect  to 
both  r  and  0  over  the  specified  variable  ranges. 
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8t(rfi)  •  i4,(0)  -  r  =  /?„( 6)  +  »,(©)  -  r  -  0 


(37) 


g,(r,6)  =  '  (K  -  ^/2)W(6)  -  (*,  -  H^/2)sin(0)  +  ts,(0)  -  r  =  0 

/or  1  £  i  £  w  + 1  ,  0  i0  sn/2 
and  [WJ2  ♦  ts,M 2)}  s  r  s  [l,  +  ts,(0)j 


(38) 


Since  ts,(0)  can  be  explicitly  solved  for  in  Eq  (38),  this  equation  proves  to  be  very 
useful  if  the  .ntroduction  of  an " artificiat '  reference  surface  is  required.  Given  the  r-0 
profile  for  one  of  the  radome’s  actual  surfaces,  i.e.,  a  set  of  (r,0)  points  which  lie  on 
the  actual  radome  surface,  a  corresponding  ts(0)  profile  is  obtainable  via  direct 
substitution  of  the  (r,0)  values  into  Eq  (38).  Numerical  interpolation  and 
approximation  techniques  such  as  Lagrangian,  Hermite,  cubic  splines,  etc.,  may  then 
be  employed  to  develop  the  required  ts(0)  function  [18]  The  only  restriction  on 
the  best  fit  approximation  is  that  the  ts(0)  function  derived  be  at  least  once 
differentiable  with  respect  to  8. 

The  GO  ray  tracing  method  employed  requires  a  surface  normal  to  be  obtained 
for  any  arbitrary  point  on  the  i*  surtace.  Since  Eq  (38)  mathematically  describes  the 
i*  surface  and  varies  only  as  a  function  spherical  coordinates  r  and  8,  a  spbencal 
coordinate  gradient  may  be  utilized  to  calculate  surface  norma!  vectors.  The 
sphencil  gradient  vector  components  G,(r,8)  and  GA(r,8)  are  next  calculated.  From 
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the  definition  for  the  spherical  gradient  of  a  function  varying  only  in  r  and  6,  Eq  (39), 
the  normal  vector  components  are  obtained. 


Vg,(r,0)  <=  “[g,(r,6)]r  +  I  A[g|(r,e)]  8 
-  Gr(r, 6)  t  +  G^(r,Q)  6 


(39) 


G,(r,  8)  =  !:[<?.  (r,6)]  =  -! 


(40) 


G“-(r’6)  =  7^igl(r>0)] 


\j>_\ 

rae 


[  y'/t;  -  kjco£(Q)  -  K,  sin(0)  ♦  K,(0)  -  r 


(41) 


GAr,  8)  = 


(l/2)(-/f?)(2cos9)(-sin9)  Kr cos0  \  d 


t  (r;  -  k]  cosJ(0) 


(42) 


G,hyfi)  = 


K,  cos(9)  (  Krsin(0) 


-  1 


I  /fi2  -  JC,W(9) 
wAere  Kr  =  Rr-  WJ2 


(43) 
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The  expressions  given  by  Eqs  (40)  and  (43)  for  G,(r,0)  and  G,h(r,6)  are  the  result 
of  taking  partial  derivatives  of  g,(r,0)  with  respect  to  r  and  6,  respectively.  These 
results  are  substituted  into  Eq  (39)  to  obtain  the  surface  normal  vector. 

3.4  Taper  Functions 

Two  specific  radome  taper  functions  are  considered,  namely,  1)  a  constant  taper 
thickness  as  measured  along  the  radome  surface  normal  direction  and  2)  an  "ideat' 
taper  function  based  on  constant  electrical  thickness  as  a  function  of  look  angle  0. 
The  constant  thickness  taper  development  is  included  here  for  completeness  and 
validation  purposes.  Radomes  with  constant  thickness  tapers  are  perhaps  the  most 
extensively  analyzed  in  literature.  Therefore,  a  wealth  of  empirical  and  measured 
range  data  exists  for  verification  of  analysis  and  modeling  results.  The  concept  of  an 
"ideat'  taper  function  is  introduced  as  a  means  by  which  an  optimum  radome  taper 
may  be  developed.  Similar  to  the  multi-layer  surface  equation  development,  the 
"ideat'  taper  function  development  is  unique  to  this  research  effort. 

Under  the  assumption  of  locally-planar  surfaces  at  ray-radome  intersection  po.nts, 
transmission  characteristics  of  Electromagnetic  (EM)  waves  can  be  derived  using 
dielectric  slab  propagation  methods.  As  an  EM  wave  propagates  through  a  dielectric 
slab  it  will  experience  a  phase  change,  i.e.,  a  phase  shift  as  it  propagates  from  the 
incident  surface  to  the  transmit  surface.  For  both  the  constant  and  "ideat'  taper 
function  developments  it  is  necessary  to  introduce  the  concept  of  "electrical 
thickness".  Letting  <}>,  represent  the  amount  of  phase  shift  as  measured  at  points 
along  the  surface  normal  direction,  the  "electrical  thickness"  d c  of  a  lossless  dielectric 
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slab  can  be  expressed  as  shown  in  Eq  (44).  In  this  equation  dh  represents  the  slab 
thickness  as  measured  along  the  surface  norma!  direction,  A  is  the  free-space 
wavelength  of  the  propagating  wave,  e,  is  the  lelative  dielectric  constant  of  the  slab, 
and  9,  is  the  incidence/design  angle  between  the  surface  normal  vector  and  wave 
propagation  direction  vector.  Letting  4>,  -  Nir  in  Eq  (44)  results  in  what  is  called 
an  "Nth-order  wall"  satisfying  the  expression  given  by  Eq  (45)  [6,  7], 


-  sin2  6,  (Wavelengths) 


(44) 


For  an  Nth-Order  Dielectric  Slab  : 

dn  =  — - -  for  N  =  1,2,3,... 

2^ -  sin2  6^ 


(45) 


Two  commonly  considered  cases  are  when  N  =  1  and  N  =  2,  corresponding  to 
"half-wave"  wall  and  "full-wave"  wall  designs  at  the  specified  design  angle  of  9,.  All 
Nth-order  wall  designs  exhibit  characteristics  such  that  reflections  from  the  transmit 
surface  cancel  reflections  from  the  incident  surface,  thereby,  resulting  in  maximum 
transmission  of  the  incident  wave.  Also,  at  the  design  angle  both  parallel  and 
perpendicular  polarization  components  experience  zero  reflection,  maximum  and 
equal  transmittance,  and  equal  insertion  phase  delays.  These  characteristics  and 
Eq  (45)  strictly  hold  for  lossless  materials  but  aie  good  approximations  for  the  low- 
loss  materials  normally  used  in  radome  structures  [7j. 
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Figure  9.  Geometry  For  Constant  Taper  Function 


3.4.1  Constant  Taper  Function.  For  constant  thickness  radomes  an  "optimum" 
combination  of  dielectric  constant  and  design  angle  8,  is  selected  by  theoretical  or 
empirical  methods.  With  N,  er  and  8,  specified,  Eq  (45)  is  used  to  calculate  the 
appropriate  normal  thickness  d0.  Relative  to  a  tangent  ogive  reference  surface,  a 
constant  taper  generates  an  additional  ogive  surface  located  either  inside  or  outside 
of  the  reference  surface.  The  new  surface  generated  must  satisfy  Eq  (46)  where  the 
subscript  "c"  is  used  to  denote  "constant"  surface  quantities  as  identified  in  Figure  9. 

gc(rc, 8)  =  '  r] cosJ0  -  r  sin8  +  W J2  -  Rc  =  0 

(46) 

-  rc  =  Kcsm6  +  J  r]  -  it]  cos2  6  for  Kc  =  WJ2  -  Rt 
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Given  the  tangent  ogive  reference  surface  parameters  R,  and  V/„  spherical 
coordinate  r,  may  be  expressed  as  given  by  Eq  (47).  From  Figure  9  it  is  apparent 
that  tsc(0)  =  rt  -  r,  for  all  0  <  0  <  ir/2.  Using  this  relationship  and  Eqs  (46) 
and  (47)  the  taper  sum  function  of  Eq  (48)  is  obtained  where  the  additional 
relationship  that  Kc  =  K,  has  been  utilized.  This  taper  sum  expression  is  substhuted 
into  the  g,(r,0)  expression  of  Eq  (38)  to  obtain  the  final  surface  equation  for  a 
constant  thickness  taper.  The  derivation  process  is  completed  by  differentiating  tsc(e) 
with  respect  to  0.  Results  are  shown  in  Eq  (49)  which  is  the  required  derivative  tenn 
for  calculating  surface  gradient  components  of  Eq  (43).  li  can  be  shown  that  Eq  (49) 
is  valid  over  the  requiied  0  range. 


! 

3.4.2  "Ideal'  Taper  Function,  if  for  every  possible  combination  of  antenna 
element  location  and  propagation  direction  the  radome  normal  thickness  d,  can  be 
modified,  i.e.,  varied  m  accordance  with  Eq  (45)  such  that  each  propagating  ray 
|  expenences  equivalent  amplitude  attenuation  and  phase  delay,  the  radome  would 

possess  "ideal'  propagation  characteristics.  The  term  "ideal'  implies  the  radome  is 
’tmnspareni"  to  an  incident  uniform  piane  wave,  i.e.,  introduces  minimal  or  no 
variation  on  either  the  amplitude  or  phase;  this  ''ideal'  situation  is  impossible  to 
achieve  under  realistic  conditions.  However,  given  a  specific  element  location  within 
the  array  and  a  desired  scan  plane,  in  conjunction  with  specified  radome  geometry, 
it  ts  possible  to  utilize  Eq  (45)  to  develop/estabiish  what  will  be  referred  to  as  an 
'ideal'  taper  function. 

Considering  the  geometry  established  m  Figure  10  where  g,(r,9)  is  a  given  surface 
equation  for  the  l-th  radome  surface,  the  goal  is  to  determine  taper  function  t,(0) 
such  that  g,*,(r,6)  can  be  represented  as  the  sum  of  g.(r,9)  and  t,(6)  for  0  <  r  <rB1J 
and  0  <  6  <rr/2.  Given  S'  is  varied  in  the  arbitrary  4>'  =  0  plane,  0  <  8'  <  rr/2, 
it  is  possible  to  calculate  an  " ideal  i  (6)  utilizing  Eq  (45).  Tne  taper  function 
developed  is  ideal'  only  for  the  element  location  under  consideration,  i.e ,  an  clement 
located  at  the  aperture  center.  This  is  a  reasonable  choice  considering  how 
amplitude  weights  are  typically  maximum  at  this  point.  Given  any  other  element 
location,  the  radome  thickness  d,  or  t,(0)  value  seen  by  an  emanating  ray  will  be  less 
than  optimum,  i.e.,  will  generally  not  satisfy  Eq  (45)  for  the  given  0  scan  direction 
and  calculated  incidence  angle  6,.  Eqs  (50)  thru  (55)  summarize  the  process  fot 
calculating  an  " ideal 1  taper  function  t.(0). 
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Figure  JO.  Geometry  For  Establishing  "Ideal"  Taper  Functions 


Given  8' ,  g.O'.O)  =  0,  and  4>'  -  0  : 
r,'  -  r  3  «,(r,0)  >•  0 
?'  -  r('  s=in0'  i  +  r’  cos9'  j 
f' 

t’  =  =  sinB'  2  *  cos0'  z 

i  •'(  i 


(50) 


For 


at  P‘  :  h: 


i  V«,(r/,e')  j 


0,  =  cos’1  { -  r,'  •  }  (51) 


Given  8(  and  sin2 Ot  -  1  -  cos2 9,  »  1  -  (('  •  j2 


t/ien  d_  -  ■ 


NX 


(52) 


2 \/e,  +  }2  '  1 


On  77ie  i+l  Surface,  £„,(r,9)  =  0, 

.  F,'  -  d,  «/' 

J>.W.‘-4A  .  ^I-TTT 


(53) 


8  =  cos'-‘(P14l  *z|  =  cos" 


F!  ~dnAi  \ 


r;  •*- W  •<)! 


I'/  I  j 


(54) 


Given  0  ,  g,(r,0)  =  0  anti  r,  =  r  3  £,(r>6)  =  0  , 

*(0)  =  l^.i  !-'•,  =  i'T  "  <*.»/!  "  r( 


(55) 


For  typical  radome  geometries  including  a  tangent  ogive  the  above  process 
seldom  leads  to  a  closed-form  solution  for  t.(0).  Numerical  techniques  are  employed 
to  calculate  0  and  t,(0)  by  varying  0 '  over  its  full  range  of  values.  Curve  fitting 
techniques  xo  are  used  to  arrive  at  a  functional  form  approximation  for  :,(0).  it  is 
evident  from  Figure  10  that  the  surface  normal  vector  fails  to  exist  at  the  radome  tip. 
Numerical  solution  techniques  must  account  for  tms  in  approximating  t,(0)  values  for 
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small  values  of  9'.  For  radomes  with  high  fineness  ratios  (greater  than  two)  this 
presents  little  difficulty  since  in  the  radonie  tip  region  the  value  of  6  rapidly 
approaches  the  value  of  8'. 

The  numerical  solution  procedure  is  applied  to  the  tangent  ogive  radome  case 
previously  cons.dered  by  letting  g,(r.9)  represent  a  tangent  ogive  surface.  A  Fortran 
computer  program  was  written  to  accept  arbitrary  ogive  surface  parameters, 
Radius  (R)  and  Width  (W),  and  arbitrary  dielectric  constants.  The  program  varies 
the  value  of  8'  over  its  range  and  generates  a  corresponding  table  of  0  and  t,(9) 
values.  For  a  test  case  the  program  was  run  with  ogive  surface  parameters  of 
R  -  234.472  X,  W  =  36.1358  X,  and  with  a  dielectric  constant  given  by  e,  =  4.8. 
The  values  of  8  and  corresponding  t,(6)  values  generated  are  plotted  in  F’gure  11  as 
indicated  by  the  symbols.  These  data  points  were  input  into  a  computer  with 
curve-fitting  capability  which  m  turn  provided  a  "best  fit"  equation;  Eq  (56)  is  the 
Gaussian  equation  returned  from  the  curve-fitting  routine.  Tne  dashed  line  in 
Figure  11  is  a  plot  of  Eq  (56)  showing  a  "near"  perfect  fit  to  the  original  input  data 
points. 


Gaussian  Function  For  "IDEAL"  Taper. 
r,(6)  -  0.2267  *  0.607 1  exp  j - .5  j  '{  | 


(56) 
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IV.  Radome  Reflection  Points 


4.1  Introduction 

This  chapter  provides  the  motivation  and  methodology  foi  calculating  wavefront 
reflection  points  on  a  radome’s  inner  surface.  A  Geometric  Optics  (GO)  approach 
is  used  to  analyze  and  predict  overall  reflection  performance.  As  noted  previously, 
a  significant  portion  of  the  discrepancy  between  measured  and  numerically  predicted 
boresigbt  errors  is  caused  by  reflections  within  the  radome  interior  [6].  Hence,  to 
obtain  accurate  and  reliable  boresight  error  estimates  such  interna!  reflections  must 
be  taken  into  account  in  the  development  of  sound  radome  models.  Primary 
reflection  points  within  a  radome  account  for  a  major  portion  of  the  total  reflected 
energy  received/transmitted  by  the  aperture.  The  following  development  concentrates 
cn  calculating  pr.mary  reflected  rays  only  as  a  first  order  approximation  to  overall 
reflection  effects.  A  primary  reflected  r?y  is  considered  to  be  any  ray  which 
encounters  at  most  a  single  transmission  point  and  single  reflection  point  while 
propagating  from  outside  the  iadoine  to  a  point  within  the  radome  on  the  aperture 
plane. 

The  geometry  and  conventions  established  m  Figure  12  are  used  throughout  the 
discussion.  Tile  point  identified  as  (xry„z,)  represents  the  "roiated"  antenna  element 
position  for  a  particular  element  in  the  apertme  plane  with  its  associated  "critical'  or 
reflection  pomt(s)  given  by  (x(,y.,zc).  Ail  critical  points  he  on  an  arbitrary  radome 
surface  specified  by  g(x,y,z)  =  0  as  indicated  in  Figure  12.  The  vector  k  represents 
the  desired  direction  of  propagation  with  components  <k*V},k,>.  For  this 


geometrical  representation,  the.  problem/objective  can  be  stated  as  follows:  "Giver 
the  point  (xnynz,)  and  the  desired  propagation  dnection  k,  determine  the 
coordinate(s)  of  the  poiin(s)  (x(,yc,z()  satisfying  g(x,y,z)  =  0  such  that  a  ray 
emanating  from  (x„ypz,)  intersects  the  radome  surface  at  (xc,);.ze)  and  leflects  in  the 
k  direction. 
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4.2  Generalized  Variational  Solution 


The  method  of  approach  to  .solve  the  problem  is  to  invoke  Fermat’s  principle  and 
then  apply  variational  calculus  to  obtain  the  critical  points.  Firs'.,  a  planar  reference 
surface  f(x,y,z)  =  0  is  established  as  in  Figure  12  with  normal  vector  k  and  located 
at  a  distance  D  from  the  aperture  plane.  As  shown,  the  point  (x,.yi,z,)  is  the  "Jina1" 

A 

intersection  point  of  a  ray  which  emanates  irorr  (:;„yc,zc),  propagates  in  the  k 
direction,  and  strikes  the  f(x,y,z)  =  0  normally.  The  problem  statement  may  now 
be  modified  by  application  of  Fermat’s  principle:  “Given  the  points  (x„y„z,)  and 
(x„y„7<),  subject  to  the  constraints  that  (x„y„zr)  be  in  the  aperture  plane  and  (x^z,) 
lie  on  the  reference  plane,  find  the  path  which  requires  minimum  propagation  time 
from  (Xpy„z,)  to  (x,,y(,z()  such  that  the  point  (xc,yc,z.)  is  on  the  path  and 
g(x,,yc,zc)  =  0."  ns  now  stated,  the  problem  is  a  classical  variational  calculus  problem. 

The  following  development  is  based  on  classical  variational  calculus  techniques. 
Conventions  and  notations  used  here  closely  parallel  those  used  by  Gelfand  ar.d 
Fomm  and  were  adopted  primarily  for  convenience  and  consistency  purposes.  As 
given  by  Gelfand  and  Fomin  [19],  the  functional  J[x,y,z]  for  a  variational  calculus 
problem  involving  Fermat’s  principle  m  three  dimensional  space  :s  as  given 
by  Eq  (57). 


*  »  /  F[x,y,z,x,y,z]  * 

3  v[x(r),y(0,z(0]  3 


(57) 
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In  Eq  (57)  i,  represents  the  fixed  "miliar  time  at  starting  point  (x„y„z,)  and 

H 

variable  t,  is  the  “final'  time  at  the  ending  point  (x(,y„z,).  The  functions  x(t),  v(t)  and 

z(t;  a-  e  parametric  representations  of  corresponding  rectangular  coordinates  x,  y,  and 

Bf 

z.  The  function  v(x.y,z)  represents  the  positional  velocity  of  the  medium  through 

1 

which  the  wave/particle  must  traverse  in  going  from  the  initial  to  fina1  points.  Since 

I 

the  problem  at  hand  requires  that  the  ray  intersect  the  radome  wall  at  (x0yc,zt)  along 

the  path,  the  functional  of  Eq  (57)  is  rewritten  as  the  sum  of  two  functionals  as  given 

| 

by  Eq  (58),  with  variable  time  *.  being  introduced  as  the  "critical"  time  at  which  the 

| 

wave/particle  reaches  the  critical  point.  The  subscripts  "1"  and  "2"  are  used  to 

I 

■ 

distinguish  between  the  two  different  regions  of  integration  and  an  additional 

1 

0$ 

constraint/boundary  condition  is  added  as  a  result.  Since  the  overall  solution  which 

I 

minimizes  J[x,y,z]  must  be  continuous,  i.e.,  a  continuous  path  from  the  initial  to  final 

points,  boundary  conditions  must  be  enforced  at  (xc,yt,zt)  such  that  x,(tc)  =  x,(tc), 

I 

yiCQ  =  and  z,(t{)  =  z2(tc). 

I 

I 

1 

BEj 

,  ,  ^ 

i 

5 

1 

B 

For  a  functional  of  the  form  given  by  Eq  (58),  the  variation  may  be  expressed  as 

fif 

given  by  Eqs  (59)  thru  (61)  [19).  A  necessary  condition  for  the  functions  x,(t),  y,(t), 

p§ 

and  z,(t)  to  minimize  J,[x,y,zj,  and  x_('),  y;(t),  and  Z;(t)  to  minimize  J,[x,y,z],  is  that 

1 

the  variation  expressions  given  by  Eqs  (60)  and  (61)  vanish  for  all  admissible 

1 

I 
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I 

1 

functions,  i.e,  6  J,[x,y,z]  =  oJ2[x,y,z]  =  0.  A  necessary  condition  for  each  variation 
to  vanish  is  that  the  system  of  Euler’s  Equations  given  by  Eqs  (62)  and  (63)  be 
satisfied  for  each  variation  6J,[x,y,z]  and  6J2[x,y,z]. 


6J[*,y,z]  =  + 


(59) 


dt 


(60) 


+  +Fu,6zi)  r  +  (F1  "Vu,  ~yAh -ifui 


tj  1 


+  [F2*£x7'Fiyt‘yi*Fuf‘h)'i " (v^rVV^a,) 


(6?) 


F->.  -  iN  =  0  •’  F>*  -  iN  =  0  ;  F'a  iN  -  ° 


(62) 


^  -  iN  - 0  ■’  -  iN  - 0  F^  -  iN  = 0  (63) 
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Assuming  a  homogenous  medium  exists  over  both  regions  of  integration,  the 
velocity  of  propagation  within  the  medium  is  independent  of  position,  resulting  in 
vi(xi-yi'zi)  =  f,  =  constam  and  v./x^z,)  =  v,  =  constant.  For  this  case,  the 
Euler’s  Equations  gi  ’cr,  cy  Eqs  (62)  and  (63)  reduce  to  expressions  given  in  Eqs  (64) 
and  (65)  after  the  appropriate  partial  derivatives  of  y„z,]  and 

are  taken.  Each  of  these  equations  represents  a  system  cf 
first-order  differential  equations.  As  expecteo  for  a  homogeneous  medium,  the 
general  solution  to  ihese  two  systems  of  equations,  and  hence  the  general  solution  to 
the  Euler’s  Equations,  is  a  set  of  linear  equations  of  :he  form  given  by  Eq  (66). 
Therefore,  any  set  of  x(t),  y(t).  and  z(t)  functions  satisfying  Eq  (66)  cause  tK  integral 
terms  of  both  <S  JjJx,^,^,)  and  *J;[xj,y  ,zJ)  to  vanish  and  are  candidate  solutions  for 
minimizing  the  functional  of  Eq  (58). 
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(66) 
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Given  that  candidate  functions  are  chosen  which  satisfy  the  system  of  Euler’s 

Equations,  the  variation  expressions  given  by  Eqs  (60)  and  (61)  can  be  recast  in  the 

‘  -r  *  • 

form  given  by  Eq  [(/)  for  each  of  the  regions  being  considered.  These  expressions 

are  further  simplified  by  noting  that  the  second  term  of  each  equation  identically 

•  •; 

*  j-  ■ 

s  - 

equals  zero  for  all  admissible  functions  x(t).  y(t).  and  z(t),  when  the  required  partial 

■  /;&?> 

' 

derivatives  are  taken  as  illustrated  in  Eq  (b8).  Ey  applying  the  results  of  Eq  (68)  to 

if  ••  •. 

Eq  (67),  the  variation  expression  given  by  Eq  (69;  is  obtained  where  the  evaluation 

i 

process  has  been  applied  at  the  upper  and  lower  limits  and  the  terms  combined  into 

; 

a  single  variation  expression. 
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"  (^i  "*1^11,  -*tFit,) 

(F^Xz+F^y^F^z^" 

*(v*v*vv4.; 


(67) 
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v^FF2  \l *!*?.*£ 

For  n  -  1,2 


=  0 


(68) 


(69) 


W(*,>v:]  =  (fuSx,  +  F,v  6y,  ♦  f^oj,  -  (r,4  6*,  -  F„  6>’,  -  F„82j  ) 


-  (  ^*2  +  f2^2  "  *2^)  '  (fz>:6xz  +  F:-V;6yZ  *  F8^( 
'  K/ 


Considering  the  initial  point  is  specified  at  t  —  t,,  i.e.,  (x„y„z.)  is  known,  and  that 
continuity  conditions  at  t  =  l£  must  be  enforced  as  expressed  in  Eq  (70).  a  final 
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simplified  variation  expression  is  obtained  as  presented  in  Eq  (71).  The  expression 
of  Eq  (71)  must  vanish  for  all  admissible  variations  about  critical  point  (x^y^x,)  and 
final  point  (x^y,,^).  As  such,  two  separate  cases  are  next  considered,  Case  !  is  for  no 
variation  about  the  critical  point  and  Case  li  L  for  no  variation  about  the  final  point. 


For  (X',yt-zs)  Ghen  ■  =  6yi  =  6r,  =  0 


By  Continuity  at  t  -  t£ :  (bx,  =>  SXjj  ;  (6y,  =  iy2j  ;  (fiz,  =  SzJ 


w&wi 


(70) 


*  *  *r2£- ^Z/J 


(71) 


For  Case  I  all  variation  about  (x„y,.z,)  is  assumed  to  be  zero  such  that 
Sx,  =  Sy,  =  4z,  =  0.  This  assumption  leads  to  the  expression  given  by  Eq  (72)  which 
when  rewritten  using  the  vector  convention  established  in  Eq  (73),  results  in  the  final 
dot  product  expression  given  by  Eq  (74).  Since  the  dot  product  of  the  twe  sectors 
identically  equals  zero,  the  velocity  and  incremental  vectors  must  be  orthogonal  to 
each  other  at  the  critical  point. 


For  6xf  *  6>y  ~  ~  0 

wiwi  -  + ('V'vk  ♦ 


(?-) 
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Using :  Fi  =  ;  F  =  ;  F,  = 

1  FU)  31  F(»)  1  F(») 


where  F(*)  =  F[x,y,z,x,y,i]  =  '/i2(r)  +  y2(r)  +  z2(f) 


Letting:  r(t)  =  [i(t),y(t).z(?)] 

z(f)  =  [iW,y(t),z(f)l 

!F(f)  |  =  t/  x2(:)  -  y2(t) ^  z2(f) 


(73) 


JW.Ml 6l  L 0m. 

f,(*)  F2(.)J  '  Ft(*)  F2(*)  /c  (f|(.)  F,(«)J  c\  (74> 


fftc)  -  WJ}- {I*.  tye,K)  =  ° 


In  <•  gruDhical  sense  incremental  vector  (<5x(,Ay(,<5zc)  repiesems  the  change  in 
position  of  the  critical  poim  on  the  inner  radomc  surface.  Therefore,  allowing  for 
variation  of  the  critical  point  requires  that  conditions  of  the  original  problem 
statement  be  enforced,  ie.,  (x(  -t-  i  xt,yc  +  s  yc,z_ + <5  zt )  must  satisfy  g(x,y.z)  --  0.  The 
results  of  applying  this  constraint  are  summarized  in  Eqs  (75)  thru  (78).  In  arriving 
at  the  finai  expression  given  by  Eq  (78)  an  incremental  form  of  a  Tay’or  .cries 
expansion  is  employed.  Enforcing  the  conditions  that  1)  the  point  (xt,yt,zc)  satisfy 
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1. 

|t§ 

-l,v. 

:  >C;V; 

g(x,y.z)  -  0  art.1  2)  the  incremental  variations  Sxc,  Sy„  and  Szr  are  very  small,  results 

•i-‘  •  v>* 

in  the  simplified  expression  for  the  dot  product  of  Vg(x,,)e,zc)  and  incremental  vector 

{ixc,5y.,£z<)  as  given  by  Eq  (78) 

.7 

Viz<)  r  0  (75) 

: 

>lv. 

VX  - 

Using  The  Incremental  Farm  of  Taylor  Series 

-1*  'it# 

i  (’6) 

f(x*h,y+k,z*m)  -f(x,y,z)  <  [hf.rkfy*aftj 

t 

:■  M 
Mi 

.  *  '.fs* 

■  fiighe,  Order  Terms  ( H.O.T) 

■"■■f.*  ;  * 
•  ■'*  \* 

i 

s{xe  " bx.  >y< *  6>V .  =  «i-c.v.,z,)*  [sM  *  gyS yc  +  g,bzc) 

*  H.o.r  -  o 

•'*?/&# 

’"i-’-C  ^  ’ 

For  yUc,yc,z{j  -  0  and  bi. < < !  ,  5yc<  - 1  ,  oz(<<]  —  H.O.T  ~ 0 

:  '’**  ’  -'  / 

f78) 

j 

The  do;  product  >>i  the  gradient  jpd  incre.t.tntal  sectors  identical!)  equals  zero. 

C^Vv.  ' 

indicating  the  two  vectors  are  Ojthogonal  to  ocv r  other.  Comp,. risen  ot  the  velocity 

6! 

•  -• 

’•A*-  ' 

_ _  _ _ -  _ _  _ _ _  _ _ 

expression  given  bv  Eq  (74)  with  the  final  gradient  expression  of  Eq  (78),  clearly 
shows  two  different  vectors  orthogonal  to  a  common  incremental  vector,  indicating 
the  gradient  and  velocity  vectors  are  coplanar  to  each  other.  This  condition  is  true 
for  all  admissible  incremental  vectors.  If  the  incremental  vanabons  S\.  <5y„  and  62, 
are  once  again  assumed  to  be  very  small  and  the  reflecting  surface  is  assumed  to 
have  a  large  radius  of  curvature  in  relation  to  the  variations,  a  locally  planar 
approximation  of  the  surface  at  the  critical  point  may  be  made.  As  such,  all 
admissible  incremental  vectors  satisfying  Eqs  (74)  and  (78)  lie  on  a  locally  planar 
surface.  This  assumption  further  restricts  the  orientation  of  the  velocity  and  gradient 
vectors  dictating  that  they  be  collinear,  the  condition  expressed  by  Eqs  (79)  and  (80) 
where  >.  and  are  non-zero  real  numbers  As  defined,  vector  nc  represents  the  unit 
normal  vector  at  the  critical  point. 


K<fe)  -  K^t>\  *  xvs{x'>y<,ie) 


(79) 


letung-  n.  -  and  Xc  =  XjVs( Wc,;gj 

|  vS(x‘’>’c-Z<l\  (80) 


?fi<)  -  fyO  =  xc  ", 


The  geometry  illustrated  in  Figure  13  is  next  considered  to  establish  the 
significance  of  th.s  relationship.  Using  the  conventions  indicated,  Eqs  (81)  thru  (84) 


62 


'S-v 


777 
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were  derived  to  show  that  the  condition  specified  by  Eq  (79)  is  equivalent  to 
satisfying  Snell’s  Law  of  Reflection  at  the  critical  point,  a  necessary  and  expected 
condition  when  applying  GO  assumptions.  Relationships  established  in  Eq  (81)  are 
obtained  by  taking  the  dot  product  and  cross  product  Eq  (80)  with  nt  as  indicated. 


Figure  13.  Geometry  for  Verifying  Snell’s  Law 


For  Case  II  a!!  variation  about  the  critical  point  is  assumed  to  be  zero,  i.e, 
SxK  =  6 yc  =  Szc  =  0,  and  a  derivation  procedure  similar  to  the  one  applied  to  the 
first  case  is  carried  cut.  A  simplified  variation  expression  for  Case  II  is  obtained 
from  Eq  (71)  and  may  be  expressed  as  given  by  Eq  (85).  As  before,  the  vector 
notation/convention  established  in  Eq  (73)  is  applied  to  obtain  the  final  dot  product 
expression  shown  in  Eq  (86). 


For  bxc  -  byc  -  6zc  -  0 


3i(x,y,z] 


(85) 


■{bXf.tyM,)  *  0  (“) 


Similar  to  the  incremental  vector  about  the  critical  pom:,  the  incremental  vector 
(fx„(5y(,6Z,)  graphically  represents  tne  change  in  position  of  the  final  point  on  ihe 
reference  surface.  Any  variation  of  the  final  point  must  contorm  to  the  original 
pioblem  statement  whicn  requires  that  the  end  point  lie  on  the  plane  specified  bj 
f(x,y,7)  =  0.  Therefore,  for  all  admissible  incremental  vectors  about  the  final  pom:, 
the  point  (x,+  £xr.y(+  6y,.z;-e  6z,)  must  satisfy  f(x.y,z)  =  0  Eqs  (87)  and  (88) 
summarize  the  results  of  enforcing  this  endpoint  boundary  condition 
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(87) 


For  f[xt  +  hx{,yl*hyi, z;+  bZ/)  '  0 
f(xf^6xryf*6y/,z/*bzJ)  =/(*/, >,.*,)  * 

b 

+  H.o.r  =  o 


Fer  f{x/,y/tz^  ~  0 , 6xf<<  1 , 8jy-»<  1  ,  6z/<<  1  -  i/.O.T.  =  0 
'fen  (,'  -0  ”  V“\X/,yr?./)'{bxf.&y/,Sz/)  *  0 


m 


Evaluation  of  f(*  y,z)  at  a  variabit  end  pom!  ,o  accomplished  by  once  again  using 
an  mere,  -ental  form  of  Taylor’s  ses  iss  expand -r.,  resulting  in  the  final  simplified  det 
product  e/piession  given  in  Eq  (88).  This  expression  cleat  !>  showi  an  orthogonal 
re.rtionship  between  the  increment,.!  vector  and  g.adient  vector  Vf(x,y,z)  at  the  final 
time.  Since  'he  incien.ental  ve,.,os  (-.x„^y,.iv,)  is  orthogonal  to  both  the  ’  elociiv  and 
gw, die  :it  vee’ors  at  the  final  time,  evident  by  comparing  Eq  (86)  with  Eq  (881.  end 
all  fums-ib!"  ir>«.remema!  vectors  must  he  cn  the  planar  refeietice  surface,  the 
ire-vircn!?!  and  grad.eit  ve--  iors  themselves  must  be  collinear  This  relationship  <s 
expressed  in  Eq  (89;  Aiiere  two  ..Jditional  far  ois  have  beer,  taken  into  consideration 
First,  the  form  of  candidate  solutions  is  a  set  of  lira  «n  f  jtiClions  given  by  Eq  (66) 
which  result  from  ctms.denng  a  homogeneous  medium  As  such,  the  velocity  vector 
remains  con, tan!  over  time.  Second  tre  reieience  plane  f(x,v,r.)  =  0  is 

A 

•.  harantenzed  by  the  problem  statement  as  t'svir.g  a  surface  normal  of  k .  \  <■.. 


✓V 

Vf(~.,y,z)  =  k  for  all  (x,y,z)  on  the  surface.  Applying  these  conditions,  the 
relationship  established  in  Eq  (89)  is  s’gnificant  in  that  it  ciearly  satisfies  the 
requirements  of  the  original  problem  statement. 


For  f^t)  -  constant  Vre[rc , tf]  and  -  £ 

then  4(f/)  =  W(*/. >/.£/)  -  tyt)  =  £  »  (kx  ky,k.) 


(89) 


At  this  point  die  necessary  conditions  and  reldtionsl  ps  have  been  established  for 
minimizing  J|xv  i'  an d  are  summarized  as  fellows,  i)  t).s  assumption  of  a 
homogeneous  medium  in  hoth  regions  uf  integration  establishes  that  only  functions 
of  ‘he  torm  expressed  by  Eq  i()6)  arc  admissible,  2)  requ.nng  continuous  solutions 
at  t  =  i,  ten  all  admissible  variation  about  the  critical  point,  white  satisfying  the 
surface  constraint  that  g(\,y„zj  =  0,  results  in  a  relationship  satisfying  Snell’s  Law 
of  Reflection  and  is  expressed  as  given  by  Eq  1,79;.  ')  a  homogeneous  medium  over 
the  second  region  of  integration  and  a  constraint  that  the  una1  point  he  on  a  planar 
reference  surface  specified  by  f(x,y,z)  =  0  with  surface  corma1  vector  k ,  results  m 
the  constant  velocity  relationship  shown  in  Eq  (89),  satisfying  the  velocity 
condition/reflectioii  direction  requirement  at  t  =  t„  and  4)  specifying  the  initirp  point 
(x„y,.z,)  results  in  a  fixed  end  point  condition  at  t  -  t,  eliminating  three  unknowns 
frem  candidate  function  solutions.  Beginning  wi,n  candidate  ..options  of  ’he  torm 
given  by  Eq  (66),  each  of  Wit  above  mentioned  conditions  and  relationsh.ps  are 
applied  to  ob'ain  a  specilic  solution  for  i^e  problem  Using  candidate  functions  of 
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the  form  given  by  Eq  ( 66)  the  fixed  end  point  boundary  condition  at  t  =  t,  is  first 
applied  to  readily  eliminate  three  unknowns.  Eq  (90)  summarizes  the  results  of 
applying  the  fixed  end  point  boundary  condition  where  the  initial  time  t,  is  arbitrarily 
set  equal  to  zero  for  convenience  purposes. 


*i(*,)  =  a*,  *  nc 
>,(',)  =  bt,  *  be 

=  Cti  +  Co  l','° 


x,(0)  --X'  =  a0 

>'i(0)  =  »  =  b0 

Z.(0)  =  z,  =  c„ 


x, (0  =  at  *xr 

y, (r)  =  bt  +  yr  (90) 
2,(0  =  Ct  +  zr 


The  vector  notation  established  in  Eq  (73)  is  next  applied  to  the  candidate 
functions,  resulting  the  unit  vector  velocity  equations  of  Eq  (95).  The  quantity  v„  is 
introduced  ;  s  the  positional  velocity  in  the  "medium"  and  is  a  result  of  assuming  that 
v,  «=  v,,  i.e.,  the  positional  velocity  over  both  regions  of  integration  are  identical. 
These  results  are  applied  to  Eqs  (79)  and  (89)  where  thtee  unknowns  are  determined 
from  Eq  (89)  and  substit  ited  into  Eq  (79),  resulting  in  the  system  of  equations  given 
in  Eq  (93). 


fit)  =  — for  rejr  ,] 
•fa1  +  b2  «■  c2  v» 

h(‘)  -  -  <A’ii9  f0r  r-k.J 

v„  *'  /J 

where :  -  ya2  +  b1  +  c2  *  -j A1  +  B‘  *  C2 


(91) 
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For  r2(tc )  =  * 


4(fc)  -  — =  (k, ,ky,kt)  ~  A  =  vmkx  ,  B  =  ,  C  =  vmk. 


For  i \(tj  -  f2(t)  =  l^g(xc,yctzc) 
7-  -  **  =  k8A*c>ye'Zc) 

m 

Tn~k>‘  Xs>(x'-y'’z') 

-7  -  K =  x^M£.i£) 


When  combined  with  an  arbitrary  reflection  surface  equation,  the  v„  expression 
0  "q  (91)  and  velocity  equations  of  Eq  (93)  represent  a  system  of  five  equations  and 
five  unknowns,  namely,  a,  b,  c,  a,  and  tt.  They  are  valid  under  the  two  previously 
stated  assumptions,  homogeneous  and  identical  mediums  over  both  rtgons  of 
integration,  and  apply  to  any  arbitrary  reflecting  surface  satisfying  g(x,y,z)  =  0  which 
is  differentiable  with  respect  to  all  variables. 


4.3  Tangent  Ogive  Solution 

A  tangem  ogive  surface  satisfies  the  stated  requirements  and  is  introduced  as  the 
refecting  surface  of  the  radome.  The  expression  given  m  Eq  (94)  represents  a 
tangent  ogive  which  has  its  base  plane  in  the  x-y  plane  and  :ts  tip  on  the  positive 
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From  Equation  90: 


*,  . 
—  -  at  *  x, 
R 


y, 


bt  +  y. 


R  C‘  * 


(96) 


From  Equation  87 : 


=  \Ja“  r  b2  ■*  c2 


(97) 


Normalizing  the  Surface  Equation : 


g(*.y.z)  =  0 
R 


/f- Z:  -  Jx2  +  f  *  K  =  0 


wft««  K  -  —  -  1 
2R 


(98) 


c>om  Equation  93 . 

(s  -  *,)  =  i  #,(*,, yt,f/l 
(6  *  >-g;[xc.yc,zc) 

(c  -*.)  =  Ag{( X',yc,it) 


(99) 


l 
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The  observation  is  now  made  the  by  introducing  a  spherical  coordm.-te  change 
of  variables  as  given  by  Eq  (100),  Eq  (97)  ts  satisfied  for  all  a  and  /9  over  the 
specified  ranges.  The  original  system  of  equations  is  now  reduced  to  four  equations, 
shown  in  Eqs  (98)  and  (99),  3nd  four  unknowns,  namely,  a,  p.  X,  and  Tc. 


a-sinacosp  ,  i>=smc!sinp  ,  c  =  cose 
where  0  <  p  <  2it  and  0  <  a  <  — 


(100) 


Further  simplifications  are  possible  by  considering  a  couple  of  additional 
constraints  on  the  reflected  rays  First,  only  reflected  rays  with  positive  h, 
components  are  of  interest.  This  results  from  restricting  the  antenna  scan/observation 
region  to  angles  which  look  "through"  the  radome.  Second,  only  reflected  rays  with 
reflection  angles  greater  than  zero  or  less  thar  nil  Rads,  as  mcasuied  from  the 
surface  normal  at  the  critical  point,  are  considered.  As  a  result  of  this  restriction 
Eq  (83)  can  be  used  to  verify  that  a  must  be  less  than  zero.  Third,  as  zc  approaches 
L,  both  x,  and  y,  approach  zero  As  a  result,  the  parual  derivatives  of  g(x,\,z)  with 
respect  to  x  and  y  fail  to  exist  ai  (O.O.L)  end  two  expressions  of  Eq  (99)  become 
invalid  and  an  insufficient  number  of  equations  rem—r.  to  solve  the  problem. 
Therefore,  z{  cannot  equal  L  and  care  must  be  taken  ns  zt  approaches  L.  Given 
k,  >(!,/.<  G,  and  0  <  zt  <  L,  the  last  expression  of  Eq  (99;  can  be  used  to  verity 
that  c  must  be  greater  than  zero,  resuming  in  the  rcstr.ctcJ  range  of  a  value.-  shown 
tn  Eq  (iQO). 

■*> 


Since  the  scaled  surface  expression  of  Eq  (98)  is  independent  of  X,  X  may  be 
eliminated  using  Eq  (99)  provided  care  is  taken  to  ensure  that  the  partial  derivatives 
exist  and  are  non-zero.  In  eliminating  X  three  special  cases  are  cons'dered.  Each 
case  is  distinguishable  from  the  others  by  the  relationship  between  the  values  of  xc 
and  yt.  For  each  case  considered  the  goal  *s  to  begin  with  £qs  (98)  and  (99)  and 
eliminate  as  many  unknowns  as  possible,  preferably  obtaining  closed-form  solutions 
for  all  variables.  Each  of  the  ihree  special  cases  is  further  divided  in  to  subcases 
based  on  raceme  geometry  and  aperture  scan  angle  restrictions.  A  total  of  eleven 
special  cases  and/or  subcases  exist.  The  following  development  is  provided  for  the 
case  which  yields  a  majority  of  critical  point  solutions.  Complete  solution  details  for 
additional  cases  and  subcases  thereof  are  provided  in  Appendix  A. 

For  the  first  case,  Case  1,  both  x,  and  yc  are  assumed  to  be  non-zero.  As  a 
result,  the  partial  derivative  expressions  of  Eq  (98)  are  substituted  into  Eq  (99)  to 
eliminate  X  as  shown  in  Eq  (101).  By  equating  the  first  two  expressions  of  Eq  (101) 
and  applying  the  change  of  variables  presented  in  Eq  (100),  the  process  of  deriving 
an  expression  for  T,  is  accomplished  as  shown  in  Eq  (102).  At  this  point  two  sub¬ 
cases  are  required.  Case  i.A  is  considered  here  and  Case  l.B  derivation  procedures 
are  provided  in  Appendix  A.  Case  LA.  assumes  that  7t  can  be  explicitly  solved  for 
from  Eq  (102)  resulting  ir.  the  expression  as  shown  :n  Eq  (103).  This  Tc  expression 
and  the  change  of  variables  g.ven  by  Eq  (100)  are  substituted  into  Eq  (96)  to  obtain 
the  critical  point  coordinate  expressions  given  by  Eqs  (104)  thru  (106)  With  both  X 
and  T,  ehmmaied,  a  system  if  two  equations  with  two  unknowns  remains,  namely  a 
and  f 
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r ' 


V-  V  -  V  ...v  v;  - .  . 


■  •  >  1  '*  71 


0  • 


; 


■*  m 


<  -  s 


dt*.i 


■. ~ -v-f  ■-  V V: 


For  Lose  l,  xc*Q  and  yc*0 


..  (£  - K)'j*\  * y2c  _  - *y)\ fsT*?c  _  (£ - K)fi~-sj 


a  -  k,  b  -  k  .  ,  , .  . 

-  {*-K)yc  =  (b-ky)ic 

C  * C 


where  xltc)  =  afc  +  ir  and  y(tc)  =  ir~  ■*  j>r 
r.sina^cosp  -^sinp}  =  xr(sinacosp  -ky j  -  yr(sinacosp  -tx) 


for  Case  IA :  sinct^cosp  -  kxSai$)  *  0 


.  x,(sina  sinp  -  jfcy)  -  yr(sina  casp  - 
sin  a  (i.ycos  P  -  ^sin  p ) 


cosp[i,(sinKsinP  -  ky)  -  y,(sina  ccsp  -  k,)\  +  _ 


(iycosp  -  Jcxsin  p ) 


-inp [i  (sine sin p  -  k  )  -  ylsuxa  cosp  -  k  )\ 

y  = - 1 — -\ - - - - - - —  ...  I  +  y 

(<:ycosP  -  P ) 


cot  a  [i,(sina  smP  -  fcy)  -  yr(sina  cosp  -  £y)| 
(it}cosP  -  ijSinRj 


m 


r  I 


M 
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One  equation  used  to  solve  for  a  and  p  is  formed  by  substituting  the  coordinate 
expressions  of  Eqs  (304)  thru  (106)  into  the  scaled  surface  expression  of  Eq  (98). 
The  second  equation  is  formed  by  equating  the  right-hand  expression  of  Eq  (101) 
with  either  of  the  other  two  expressions,  as  shown  in  Eqs  (107)  and  (108),  with  x„ 
y.,  and  zc  replaced  by  Eqs  (104)  thru  (106),  respectively,  in  choosing  which  of  the 
two  expressions  to  use  consideration  should  be  given  to  the  limiting  cases  as  either  xt 
or  vc  approach  zero.  Note  that  the  final  expressions  to  be  solved  are  non-linear 
transcendental  type  equations  with  no  explicit  solution  possible.  Therefore,  numerical 
solution  techniques  must  be  implemented  to  solve  for  remaining  unknowns. 


[g  -  *  9l  (e-kgr-"Z 


(107) 


(g  -  kjjg  >  f:  _  (c  -  k){T-7c 

y< 


(108) 
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V.  Propagation  Technique 


5 1  Introduction 

When  ar.  Electromagnetic  (EM)  wave  travels  from  a  source  and  arnves  at  a 
distant  observation  point,  the  medium  through  which  it  propagates  will  general!) 
introduce  amplitude,  phase,  .md  polarization  variations.  As  stated  in  the  Background 
section,  all  three  variations  are  expected  to  occur  as  an  EM  v.ave  propagates  through 
a  radome.  Surface  Iniegution"  and  "Ray- Tracing"  techniques  can  be  used  to  account 
for  the  variations.  This  rcseatch  effor.  capitalized  on  existing  propagation  analysis 

| 

!  and  modeling  techniques  which  were  modified  and  extended  to  achieve  specific 

I 

I 

!  research  goals,  ip.  selecting  a  wav*  propagation  technique  two  key  faciois  had  to  be 

i 

|  taken  into  consideration.  Firm  the  propaganda  technique  chosen  for  analysis  and 

i  modeling  had  to  provide  flexibility,  i.c.,  the  technique  had  to  allow  for  key  system 

I  parameter  variations  while  providing  a  response  which  was  both  realistic  and 

: 

techmiaily  sound.  Second,  responses  to  mch  variations  had  to  occur  in  :•  timely 
manner,  i.e.,  consideration  hud  to  be  given  to  minimizing  compuiational  intensity 
whue  providing  accurate  results.  In  conjunction  with  the  dtsue  to  limit  piopagatioit 
teennique  development,  consideration  of  these  kej  factors  led  to  selection  of  tne 
"Ray-Tracing  (Geometric  Opti-s)  propagation  technique  for  analysis  and  modeling 
Purposes. 

A  parametric  investigation  of  radome  analysis  methods  [.'j  developed  a  general 
theory  for  iadon>e  analysis.  The  investigation  determined  relative  at  curacies  and 
ranges  of  validity  for  three  computer  uict  o  analysis  methods  -  surn.ee  mtegtation 


ray-trace  receive,  and  ray-trace  transmit.  The  study  concluded  that  the  ray-trace 
receive  method  was  tSie  'most  attractive"  because  of  fast  computation  time  and 
reasonably  accurate  results.  Following  this  effort,  Kozakoff  [5]  conducted  an 
improved  radome  analysis  methodology  study  which  incoiporated  reflected  E-Field 
contributions  in  the  ray-trace  receive  analysis,  achieving  'reasonably  good" 
computational  accuracy  for  most  applications  where  smer.na/radome  dimensions 
exceeded  5  wavelengths.  The  ray-trace  receive  technique  was  later  pursued  by 
Klemer  [6]  who  characterized  BSE  sensitivity  ■'behavior  using  multiple  radome  designs 
(various  wall  thicknesses  and  dielectric  properties)  and  aperture  characteristics  (size 
element  type,  amplitude  taper,  eic.),  establishing  a  ’consensus  that  the  rav-trace 
receive  formulation  offers  the  best  comprom.se  between  accuracy  and  computational 
effort  for  moderate  to  large-sized  radomes,  i.c,  radoir.es  wnh  dimensions  greater  than 
ten  wavelengths,"  typically  the  case  for  most  aircraft  radomes.  Also,  radome?  of  this 
type  are  generally  characterized  by  “fineness"  maos,  a  intio  of  radome  length  to 
radome  base  widih,  in  the  range  of  two  to  three.  These  two  limitations  ensure  that 
“iocally-plu.  -ar"  conditions  are  met  at  typical  radar  frequencies 

Considering  an  end  goal  of  the  research  effort  was  to  analyze  and  model 
mcnopulse  radar  tracking  performance,  the  ray-’race  receive  propagation  technique 
proved  to  be  very  useful.  Monopuisc  tracking  performance  is  typically  characterized 
through  BSE  measurement  procedures  conducted  on  radar  test  ranges  ^  typical  test 
procedure  for  characterizing  tracking  performance  begins  by  scanning  the  radar 
aperture,  either  electronically  or  aiechan.  tally,  in  a  specified  directum  This  direction 
:s  Designated  as  the  system  "iracking/pomting"  direction  The  location  of  a  far  field 
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source  is  varied  about  the  fixed  "tracking"  direction  m  a  predetermined  scan  pattern 
until  the  monopuise  error  signal  goes  to  zero,  indicating  the  energy  received  by  the 
aperture  is  arriving  from  a  direction  coinciding  with  the  systc  m  "tracking”  direction. 
The  final  source  location  causing  this  condition  is  designated  as  the  "actual"  source 
location.  Optical  instrumentation  test  data,  obtained  from  sources/sensors  located  on 
the  apertu'e  and  source,  is  used  to  calculate  the  angular  difference  between  the  radar 
"tracking"  direction  and  "actual"  source  location  This  difference  represents  the 
"system"  BSE  in  the  specified  '  tracking"  direction  This  process  is  easily  duplicated 
for  analysis  and  modeling  purposes  using  the  ray-trace  receive  propagation  technique 
The  ray-trace  receive  propagation  technique  is  applied  to  both  direcl "  and 
"pnrvarv  refiec-ea  rays,  ^s  preciously  defined,  pnmars  reflected'  rays  are  rays  which 
experience  at  most  one  reflect.on  priu.  ...  inteuepung  the  aperture  plane  Both 
direct  and  reflected  rays  experience  tne  effects  of  propagation  through  the  radome 
wait,  i.e.,  ampiuude.  phase,  ana  pclanzatiro  distortion  upon  tra.nsmissie  All 
'pnmar  "-fleeted"  rays  experience  .laditionnl  distortion  upar  reflecting  ‘  on, ,  .  ,nner 
radome  suna.e  At  each  ray-r.idome  imersecuon  point,  resrardless  wr.ich  ivpe  4  ray 
is  being  considered,  a  "hzclh-pianai  upproxima*'oi  is  nude  n.'  plant  ,  f 
incidence '  is  established  The  '■■..dh-pic.'-.jf  .oreiti. .  is  recurc.1  or  ,  i!-.-...u,r,g 
complex  :rt  emission  ana  refect, or  s_“  cient  .ip:  -ox.-u-.. :  -  >•  .  "iut,  f,  r  wuh 
We  es'annshct.  mai  o>r::acnti  :v.  .r:ni  mat:  "’v  .  .:r’~  a.  _n  .ruience  the 
■meets  he  r-uone  ■‘..-•.re  :s  ,  a.--  no  -tie  -.  ">v 


'  A" 


,  ....  .  '  rJ?  I,  T?  -A  .  .-  '■  ,  '  <  -i  ~  ^ 


■m 


-A 


5.2  Local  Plane  of  Incidence 

The  local " plane  of  incidence’'  is  defined  in  the  following  manner.  Given  a  point 
on  a  surface  with  normal  vector  n  and  an  incident  E-Fieid  at  that  point  propagating 
m  direction  k ,  the  "plane  of  incidence"  is  the  plane  containing  both  n  and  k .  Tne 
specific  unit  vector  pair  chosen  for  E-Field  decomposition  is  dependent  or.  particular 
analysis  and  modeling  objectives,  varying  according  to  tlte  propagatior  echn:_ue 
implemented 

The  definitions  stated  in  Eqs  (!0d|  tnru  Jill)  and  depicted  in  Figure  14  are  the 
orthogonal  unit  vector  pairs  chosen  to  define  the  "plane  of  incidence"  and  thr 
corresponding  transmission  and  reflection  unit  vectors  These  particular  vector  par 
combinations  were  primarily  chosen  for  compatibility  purposes.  As  developed  by 
Munk  [20]  and  later  used  in  development  ot  the  PMM  code  [11],  these  definitions 
were  arbitrarily  chosen  such  that  the  reflection  coefficient  lor  the  E-Field  equate 
minus  one  for  the  case  when  the  dielectric  boundary  is  a  perfect  conductor  These 
definitions  and  conventions  are  used  thrcechout  tne  remainder  of  the  dissertation 


For  the  mcid'i 


For  the  'tansmitted  ray, 


5  3  Direct  Ray  Propagation  Path 

Tht  geometry  established  in  Figure  15  is  used  m  analyzing  the  propagation  effects 
along  u  direct  ray  path  For  the  geometry  shown,  the  following  definition.;  apply: 

EJ  =  incident  E-Field  on  the  element  due  to  a  direct  path 
E,  =  Reference  E-Field  representing  a  uniform  plane  wave. 
er  =  Complex  unit  polarization  vector  of  reierence  E-Field. 

E,  —  Component  of  b,  responsible  fo:  pro:  ..gation  in  k  direction. 


E,  =  Incident  E-Field  op  outer  tadome  •>  rfr,. 

E,  =  Transmit  E-Field  on  .mi-  dome  surface 
INpi  =  ii:cident/„uter  surface  ,  •<  <r  pt> r 
TXpt  =  Tra"smu'mner  surface  ■--'w.  ;on.. 

k,  =  Incident  E -Fif'd  propagation  direction 
k,  =  Transmit  E-Field  prcpags'ion  u  -• 

d,  =  Propagation  path  length  fro—  ~*e  er.  -r  <1  pi  x  urpi  sr 
the  outer  rsdome  surface. 

do  =  Propagation  path  length  frotr  ~  i  ntv:  .a.  me  s  rface  to 
arbitrary  aperture  elemen'  -oordin.  'es 

Given  an  arbitrary  aperture  element  locaie:  ‘.angular  -dm-it-  ...  ,.z. 

<nd  a  reference  E-Fieid  located  a  distance  D  tr.~  **•:  ..perm--  _  a.  f.o.n.  „„  tne 
incident  E-Field  E®  incident  on  the  elemen  due  .  -opagattor  dong  a  direct  path 
needs  to  be  determined.  The  fina  expressiu'  de-.e  -.'ilfyE'T  s  .ccnu  '  .  ■  wo 
important  analysis  and  modeling  cna-acte-istics  Fits.  -  ausr  'r,-_-e  w  -um 
is  known  to  adversely  affect  monop.ise  tra^...rg  per  —  ..nee,  ail  p.iase  ue...-  non., 
the  direct  propagation  path  must  be  n-curai- »  ..cc.  aucs  t.  Seconu.  ..m.  ur.alys.. 
involving  arbiti urily  polarized  refe--r.ee  E-T  e.u.  s  desired, ....  fii.a.  cvp.c*».<..r.s  must 
be  gen  ral  in  nature  reflecting  a  d  oendance  or.  'eference  lie  <j  polar  rum  tue'- 

The  following  equations  summ.i  :ze  the  re  ...its  !  using  ."r  ray-tt.>'  -;er-t- 

prnpagation  technique  to  arrive  ai  the  E:,  e-rretsion  The  'error  ..orr.r.  it-' 
of  £,  responsible  lor  propagation  in  direr,  r  *  .  fir-t  determined  .;  jccore_.net 
with  Eq  (112)  and  propagateo  to  the  outer  rcaorre  .urmce  :.  -r n  the  inciornt 


h  i  ic-d  E,.  .  v  stu  wn  in  Eqs  (113)  and  (114)  a  free-space  phase  delay  corresponding 
jistaiKe  c,  i  ten  added  to  account  for  the  propagation  phase  delay.  The 
jjontn  ..  nt  fret-  pace  phase  constant  and  is  equivalent  to  2n  divided  by  the 
f'»“-spas wave  .  -^gth  X,,  (Kc  -  2tr/X0). 

G'ven  E  =  Eter  then  =  Jsr{er  -  0^2) 

l  =  V***  =  E,‘'JKAitr  -  [K'tyi)  iU3) 


Defining  Et  -  Ere  ,K’i'  and  et  -  "if,) £, 


(114) 


Decomposing  E,  miu  parallel  (/?)  and  perpendicular  (j )  components  as  established 
by  the  ,-c.ii  plane  of  incidence. 


El  -  E,(e,  •u‘,,in‘u  =  £,*£« 


K  -■  £,(v"lX  r  K 


(115) 


Propagating  parallel  and  perpendicular  components  through  the  radorne  and 


weigniing  the  components  bv  planar  transmissivity  approximations  results  tn  the 
expression  given  by  Eq  (116). 


£t  *£,''  + £,x  -  +  TiE,iu' 

*  £.{T«(*i  •«»)“« + 


(116) 


Calculating  the  portion  of  E,  responsible  for  propagation  in  direction  k,  and 
utilizing  the  definitions  for  the  nansmit  ray  plane  of  incidence  unit  vectors  the 
following  relationships  hold.  The  propagating  E-Field  component  is  seen  to  be 
equivalent  to  E,  as  shown  in  Eq  (119). 


-(*<•*.)*<}  (117) 


“‘x'K  »  *  0  ■  =  (118) 

\hXKA 


+  T.(«. -“I)"!}**.  = 0  -  K  = E.  (119) 


The  transmitted  E-Field  E,  is  next  propagated  to  the  aperture  element  location 
as  summarized  in  Eqs  (120)  thru  (123)  wheie  the  final  expression  for  E”  is  as  given 
by  Eq  (123).  The  Final  expressions  for  both  Eqs  (121)  and  (122)  are  obtained  by 
noting  th.  the  dot  product  of  A. ,  with  both  incident  plane  unit  vectors,  parallel  and 
perpendicular,  identically  equals  to  zero. 
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Ed  =  E.e'lK,i> 


-  E,ciK‘^-d^  {Ta(?, . «;,)«;  ♦  Ti(?i .  <)«: } 


«i*“«  =  {‘r  -  («,  *  ^i}  *  >2»  =  «,  •  “5  -  («",  *  ){K  ‘  «i) 


(120) 


(121) 


S,  *  fii  =  K  -  («,  ‘  £iK}  (<\  *  *,  )(*,  *  «1) 


(122) 


=  £t  where 


Ed  -  e 

«?=T.(«,*“i)^+TJL  (V“1K 


(123) 


The  final  E”  expression  in  Eq  (123)  reflects  the  two  important  characteristics 
previously  discussed,  polarizatio  i  dependance  and  accurate  phase  tracking.  An 
additional  dtsm8,  factor  not  appearing  in  Figure  15  is  included  in  the  total  phase 
delay  term  of  the  complex  expression.  Although  this  additional  phase  delay  factor 
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is  not  evident  in  previous  analysis  and  modeling  efforts,  its  inclusion  in  this 
development  is  believed  reasonable  and  justifiably  correct  since  noil-trivial  tapers  in 
radome  thickness  are  considered.  When  considering  various  techniques  for 
approximating  complex  planar  transmission  coefficients,  i.e.,  PMM,  Richmond’s 
equations,  etc.,  care  must  be  taken  to  ensure  that  the  phase  of  the  calculated 
transmission  coefficients  is  accurately  interpreted.  In  most  cases  the  calculated  phase 
represents  the  total  phase  delay,  often  referred  to  as  the  "Insertion  Phase 
Delay"  (IPD),  of  the  planar  phase  front  as  measured  along  the  surface  normal 
direction.  Referring  to  Figure  16,  this  delay  corresponds  to  the  difference  between 
the  phase  of  the  normal  phase  front  phase  <)>,  and  the  incident  phase  front  4>,. 

Given  a  specific  set  of  incident  (INpt)  and  transmit  (TXpt)  points  as  indicated 
in  Figure  16  the  actual  IPD  corresponds  to  the  difference  between  the  transmit  phase 
front  4>,  and  incident  phase  front  <j),.  As  shown  in  Figure  16,  cj),  can  be  obtained  by 
adding  an  additional  free-space  phase  delay  corresponding  to  the  distance  x. 
Eqs  (124)  and  (125)  summarize  the  process  for  calculating  <J>,  where  the  following 
definitions  are  appiied  to  the  geometry  established  in  Figure  16  : 

d„  =  Distance  as  measured  along  the  surface  normal  direction. 
dp  =  Distance  as  measured  along  the  propagation  path, 
d,  =  Distance  trotn  norma!  transmit  point  to  TXpt. 

6,  =  Ray  incidence  angle. 


Figure  16.  Additional  Phase  Delay  Geometry 
For  dt  =  ]Jdp  -d*  And  x  =  dlsmOl  (H4) 

♦,-♦.-*,<**1118,  (125) 

5.4  Reflected  Ray  Propagation  Path 

For  analyzing  propagation  effects  along  reflected  ray  paths  the  geometry 
established  in  Figure  17  is  used  As  with  the  direct  ray  propagation  case,  the 
E-Field  Ej  incident  on  the  element  due  to  a  reflected  ray  path  must  be  determined. 
Accurate  phase  tracking  and  reference  E-Field  polarization  dependance  are  once 
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again  important  characteristics  to  be  contained  in  the  final  expressions.  Vector 
definitions/conventions  for  all  components  along  the  path  from  the  reference  E-Field 
plane  to  the  radome  reflection  point  RFpt  are  unchanged  from  those  used  in  the 
direct  ray  propagation  analysis. 


Figure  17.  Reflected  Ray  Incident  E-Field  Geometry 


For  components  not  previously  defined  the  following  definitions  are  made  as 
illustrated  in  Figure  17: 

E*  =  Incident  E-Field  at  the  element  due  to  reflected  path 
E“  =  Incident  E-Field  at  the  reflection  point  due  to  direct  path. 
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E,  =  Reflected  E-Field  on  inner  radome  surface. 

RFpt  =  Reflected/inner  surface  reflection  point. 

r*. 

k,  =  Reflected  E-Field  propagation  direction. 

n,  =  Surface  normal  vector  at  RFpt. 

Results  of  the  direct  path  analysis  are  applied  directly  to  Figure  17  for  obtaining 
the  incident  E-Fteld  Er  at  RFpt.  Since  the  aperture  element  location  was  arbitrarily 
chosen  for  the  direct  path  development,  the  direct  path  results  are  applicable  for  any 
point  located  within  the  radome  structure,  including  points  not  located  on  the 
aperture  plane.  Reflection  point  RFpt  satisfies  this  condition  and  Eq  (123)  may  be 
used  to  express  E[?  as  shown  in  Eq  (126) 


=  Eg  e°  where 


Eg  ^  E 

4  =  T„(ef.«;)i;+Ti(cr.«:)«i 


(126) 


At  the  reflection  point  a  " locally-planar "  assumption  is  made  and  a  "plane  of 
incidence"  established  m  accordance  with  Figure  14  definitions.  Consistent  with 
definitions  given  by  Eqs  (109)  thru  (111),  Eqs  (127)  and  (128)  define  the  orthogonal 
unit  vectors  used  for  decomposing  E“  at  RFpt.  The  unit  vector  designations  have 
been  changed  from  u  to  v  to  avoid  confusion  between  the  transmission  point  and 
reflection  point  plane  of  incidence  vectors. 
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”r  Xk, 

\nrxic,\ 


(127) 


nrX^,  .r  £ 

v*  =  — — 7-  ;  v«  =  -v*xK 

\nrxi,\ 


(128) 


Using  these  definitions  the  incident  E-Field  Er  is  decomposed  in  to  parallel  and 
perpendicular  components  at  the  inner  radome  surface.  After  decomposition  each 
component  is  weighted  by  its  respective  reflectivity  coefficient  resulting  in  the 
reflected  E-Field  E,  expression  as  given  by  Eq  (129). 


£•  _  f»3  pi  p  «r  p  pi/i 

S  ‘  =1i£*  va  +  *  x  £jt 


(129) 


The  reflected  E-Field  E.  is  then  propagated  to  the  aperture  element  accounting 
for  a  path  length  phase  delay  corresponding  to  distance  d3.  resulting  in  the  final 
expression  for  E*  as  given  Eqs  (130)  thru  (132).  The  final  expression  for  e*  given 
by  Eq  (132)  is  obtained  after  the  dot  product  operations  indicated  in  Eq  (131)  are 
carried  out. 
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%  -  -  ^{T„(^-v')v;  ♦  TA(aMK}  (130) 


£”*/?  p  R  ~*R  i 

t  -  Et  t,  where 


ER  -  £  e->KM  *4  *4i “9i) 

?f  =  r,(?".vi)v;+rx(^.vi)v: 


(131) 


<'  "  r»  V"  {  T»  (*r  *  ‘  «i)  *  (*,  ‘  <i)K  *  ♦.*)} 

+  r*  *1  { T,  («,  *  «a)(u£  •  «;)  +  Tx  (e,  •«!)(«{•  *!)} 


(132) 


5.5  Refractive  Effects 

Typical  ray-tracmg  propagation  techniques  neglect  the  effects  of  refraction  in 
propagating  a  wave  through  the  radome  structure.  Such  effects  are  generally 
assumed  to  be  negligible,  thereby  simplifying  the  analysis  and  modeling  development. 
The  simplified  analysis  and  modeling  procedures  are  unable  to  accurately  characterize 
the  second-order  effects  of  net  ray  deflection  on  overall  boresight  error  [6],  The 
following  development  summarizes  the  process  by  which  refractive  effects  are 
analyzed,  modeled  and  characterized  under  the  research  effort.  Refractive  effects  are 
mcorpoiated  into  the  propagation  model  at  two  different  levels  For  the  first  level 
of  development  refractive  effects  are  included  only  within  the  radome  structure.  At 
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this  level,  the  transmit  point  (TXpt)  is  calculated  from  a  given  incident  point  (INpt) 

A 

and  incident  propagation  direction  k,  accounting  for  total  refractive  effects  between 
radome  layers.  Locally-planar  approximations  are  enforced  such  that  the  output 

A 

propagation  direction  k ,  remains  unchanged  from  the  input  propagation  direction. 

The  second  level  of  development  accounts  for  "overall"  refractive  effects  along  ray 

propagation  paths,  including  a  transmit  propagation  direction  k ,  which  differs  from 
/\ 

k,  and  is  dependant  on  radome  electrical  properties  and  geometry. 

Two  key  factors  are  required  to  accurately  account  for  refractive  effects  within 
the  radome  structure.  First,  refractive  effects  across  boundaries  with  discontinuous 
electrical  and  magnetic  properties  must  be  accounted  for.  These  effects  are  well 
established  and  are  obtained  through  Snell’s  Law  of  Refraction.  For  ncn-conductive 
materials  Snell’s  Law  of  Refraction  can  be  expressed  as  given  by  Eq  (133)  [SJ. 


Hirsute,  =  y' -*,e,  sin0, 


M, 

n|  m. 


sin6. 


(133) 


In  Eq  (133)  both  the  incident  angle  0,  and  transmit  angle  6,  are  measured  relative 
to  the  surface  normal  vector  at  the  incident  point.  The  complex  permittivity  er  and 
complex  permeability  jx,  characterize  the  electrical  and  magnetic  properties  of  the  two 
media.  Generally  0,  is  a  complex  transmission  angle  which  becomes  purely  real  when 
lossless  media  are  involved.  Foi  typical  radome  material',  dielectrics  with  very  low 
loss-tangents  (5,  <  <  1),  the  magnitudes  of  c,  and  ii,  are  used  in  Eq  (133)  for 
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approximating  the  transmission  angle.  Other  cases  involving  materials  with 
appreciable  loss  may  be  handled  by  computing  an  "equivalent"  real  transmission 
angle  [8], 

The  second  key  factor  is  that  each  surface  of  the  radome  structure  be  accurately 
defined.  Surface  equations  and  gradient  components  from  the  multi-layer  radome 
development  of  Chapter  III  satisfy  this  requirement.  Each  layer  of  the  radome  being 
analyzed  is  uniquely  defined,  including  a  surface  normal  vector  at  any  arbitrary  point 
on  a  specified  radome  surface.  With  each  surface  specified  and  both  INpt  and  k, 
given,  Eq  (133)  is  repeatedly  applied  to  the  radome  structure  until  TXpt  is 
determined.  Under  locally-planar  assumptions  surface  normal  vector  directions  at 
each  surface  intercepted  between  INpt  and  TXpt  equal  the  normal  vector  direction 
at  INpt.  Under  these  conditions  the  ray  propagation  direction  remains  unchanged 

A  A 

in  passing  through  the  radome.  i.e..  k,  =  k,.  Therefore,  the  first  level  of 
development  accounts  only  for  refraction  within  the  radome  structure  and  not  along 
the  entire  propagation  path. 

The  second  level  of  development  accounts  for  "overall"  refractive  effects  along 
the  propagation  path  Locally-planar  approximations  are  employed  only  as  a  means 
for  approximating  transmission  and  reflection  properties  of  the  radome  They  are  not 
used  to  approximate  surface  normal  vector  components  along  ray  propagation  paths 
Rather,  normal  vector  components  are  calculated  at  distinct  ray/surface  intercept 
points  using  the  multi-layer  surface  equations  developed  in  Chapter  III.  Therefore, 
repeated  application  of  Snell's  Law  of  Refraction  to  a  ray  incident  on  the  radome’s 
outer  surface  results  in  1)  a  transmit  point  (TXpt)  on  the  inner  surface  which  is 
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displaced  from  the  transmit  point  resulting  from  a  ray  passing  straight  through  the 

A 

radome,  and  2)  a  transmit  propagation  direction  Hr,  which  differs  from  the  incident 

A  A  A 

propagation  direction  kt,  k,  and  k,  may  in  fact  lie  in  separate  planes. 


The  geometry  of  Figure  18  shows  how  ray  refractive  effects  are  incorporated  in 
ray-trace  receive  analysis.  The  reference  E-Field  Er  represents  a  Uniform  Plane 
Wave  (UPW)  incident  from  a  distant  source.  A  ray  is  traced  from  the  reference 
E-F;eld  location  in  direction  k,  to  the  incident  point  INpt  on  the  radome’s  outer 
surface.  Snell’s  Law  of  Refraction  is  repeatedly  applied  to  the  radome  structure  until 
the  transmit  point  TXpt  is  established  on  the  inner  surface.  At  TXpt  a  final 
apphcation  of  Snell’s  Law  determines  the  transmit  propagation  direction  k ,  The  ray 
is  propagated  from  TXpt  to  the  point  where  it  intersects  the  aperture  plane.  This 

94 


1 


process  is  continued  for  each  ray  emanating  from  the  reference  E-Field  plane. 
Because  refractive  effects  have  been  included,  a  uniformly  spaced  bundle  of  rays 
emanating  from  the  reference  plane  intercept  the  aperture  plane  at  non-uniformly 
spaced  intercept  points/element  locations.  This  requires  aperture  integration  to  be 
performed  on  a  non-uniform  grid  of  sample  points. 


Since  aperture  integration  is  required  over  uniformly  spaced  or  known  element 
locations,  consistent  with  the  overall  system  analysis  and  modeling  approach,  an 
alternate  propagation  technique  must  be  established  Given  E,  is  a  Uniform  Plane 
Wave  (UPWj,  all  rays  emanating  from  the  reference  plane  are  initially  weighted  by 
equivalent  amplitude  and  phase  terms  regardless  of  their  initial  location  on  the  plane 
Figure  19  represents  an  equivalent  method  for  incorporating  refractive  effects  into 
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the  overall  propagation  model.  A  non-uniformly  spaced  bundle  of  rays  is  propagated 
from  the  reference  plane  which  in  turn  intercept  the  aperture  plane  at  uniformly 
spaced  points.  By  invoking  reciprocity,  these  propagation  paths  are  easily  calculated 
from  known  aperture  element  coordinates.  The  multi-layer  surface  equations  of 
Chapter  III  in  conjunction  with  Snell’s  Law  of  Refraction  facilitate  the  use  of  an 
iterative  Newton’s  method  for  calculating  the  ray  paths  via  ray  transmission 
techniques. 

For  calculating  INpt  and  TXpt  by  a  transmit  technique  a  specific  element  location 

A  A 

and  desired  propagation  direction  k ,  are  given.  An  estimate  of  -k,  is  determined 
and  Newton’s  method  is  used  in  conjunction  with  the  m  ter  radome  surface  equation 
to  calculate  a  candidate  TXpt.  The  ray  is  then  propagated  from  TXpt  to  INpt  by 
repeated  application  of  Snell’s  Law  resulting  in  an  estimate  of  -k,.  This  estimate  is 
compared  to  k,  by  a  simple  dot  product  operation  to  determine  its  relative 

A 

"closeness ".  If  not  within  established  limits  the  -  k ,  estimate  is  updated  and  the 
process  is  repeated.  This  iterative  procedure  is  continued  until  the  calculated 
estimate  of  -  k,  satisfies  established  equality  limits.  Once  equality  limits  are  satisfied 
the  corresponding  INpt  and  TXpt  points  are  u^ed  for  ray-trace  receive  calculations. 
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VI.  Reference  E-Field  Development 


6.1  Introduction 

The  Reference  E-Field  E,  is  a  key  element  used  in  analyzing  and  modeling  the 
overall  radome,  radar,  and  monopulse  processing  "system".  Accurate  tracking  of 
wavefront  polarization  along  propagation  paths  is  essential  to  analyzing  and 
understanding  depolarization  effects  of  various  system  components.  Specifically, 
radome  depolarization  effects  must  be  isolated  from  antenna  element  polanzation 
sensitivity  effects.  Development  of  the  incident  E-Field  expressions  in  Chapter  V  is 
indicative  of  (he  level  of  effort  required  to  maintain  polarization  "purity",  i.e.,  no  loss 
of  generality  in  regard  to  E-Field  polarization  states.  Reference  E-Field  analysis  and 
modeling  must  account  for  two  important  considerations.  First,  location  flexibility 
must  exist,  i.e.,  E,  must  be  arbitrarily  located  anywhere  within  the  system’s  field  of 
view.  Variable  E,  locations  are  required  for  aperture  far  field  pattern  generation  and 
monopulse  tracking  error  charactei  ization.  Second,  E,  must  possess  arbitrary 
polarization  states.  Tins  requirement  allows  for  radome  depolarization  and  system 
polarization  sensitivity  effects  to  be  fully  characterized. 

6.2  Reference  E-Field  Location 

Two  separate  coordinate  systems  are  considered  ir.  developing  expressions  for  an 
arbitrarily  located  reference  E-Field.  Aperture  far-field  radiation  patterns  and 
monopulse  tracking  error  characteristics  are  generally  established, 'characterized  in  a 
coordinate  system  relative  to  the  aperture  scan  direction.  Reference  E-Field 
locations  used  for  analysis  and  modeling  must  vary  in  accordance  with  this  same 
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relative  coordinate  system.  Given  that  an  overall  "system"  is  analyzed  and  modeled 
an  absolute  coordinate  system  is  used.  Therefore,  final  expressions  for  reference 
E-Field  location  are  developed  such  that  1)  they  reflect  a  dependance  on  the 
aperture  scan  direction  angles  and  2)  are  referenced  to  the  absolute  coordinate 
system  being  used. 


2 


Figure  20.  Primed  and  Unprimed  Reference  E-Field  Coordinate  Systems 

The  absolute  coordinate  system  used  for  developing  E-Field  location  expressions 
is  unchanged  from  previous  developments  and  is  represented  as  the  unprimed 
rectangular  coordinate  system  shown  in  Figure  20.  The  radome  orientation  remains 
fixed  with  its  base  located  in  the  x-y  plane  and  its  tip  on  the  z-axis.  Prior  to 
gimballing  and  scanning  the  radar  aperture  is  in  the  x-y  plane.  Given  aperture  scan 
angles  8.  and  4>,  and  gimbal  point  coordinates  (0,0, zt),  the  primed/relative  coordinate 
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system  shown  in  Figure  20  is  established.  The  location  of  E,  is  varied  relative  to  the 
scanned  apertuie  with  location  angles  specified  by  0j  and  <j>)  in  the  primed  coordinate 
system.  Location  vector  r '  for  E,  may  be  expressed  as  shown  in  Eq  (134)  where  unit 
vectors  x',  y',  and  z'  define  the  primed  coordinate  axis.  These  unit  vectors  are 
converted  to  unprimed/absolute  coordinate  system  components  by  application  of  the 
rectangular  coordinate  rotation  equations  of  Eq  (B.l),  resulting  in  the  equations 
summarized  in  Eq  (135). 


t  -  x  x  i-  y  y  *  z  z 

=  sin(0r')cos(4>/)x  +  sin (8/) sin (4>/))’’  +  cos(0r')z' 


(134) 


Converting  to  Unprtmed  Coordinates 
£’  =  (cos2 ^ cos 9,  +  sm24>,)i  -  cos<f>jsin4>J(l  -  cos05)y  -  cas^sinQ^^^-^ 
y  =  -  cos  4>,  sin^U  -  cos9a)i  +  ( sin2  ^  cos  0s  *  cos2 4^)  y  -  sm4>Jsin0J£ 
z  -  -  cos  4>s  sin  0,  i  -  sim})Jsm0Jy  +  cos95z 


Given  aperture  scan  angles  8y<f)t  and  relative  E-Field  location  angles  0 '/<{)'„  vector 
expressions  of  Eq  (135)  may  be  substituted  directly  into  Eq  (134)  to  obtain  the  final 
location  of  E,  in  unprimed  rectangular  coordinates.  The  final  expression  obtained 
reflects  the  desired  dependance  on  primed  reference  field  location  angles  and 
unpnmed  aperture  scan  angles. 
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6.3  Reference  E-Field  Polarization 

For  defining  an  arbitrary  polarization  state  of  E,  the  concept  of  a  complex  unit 
polarization  vector  is  introduced.  Utilizing  a  complex  polarization  vector  allows  for 
an  unlimited  number  of  polarizations  states  to  be  uniquely  specified,  including  linear, 
elliptical  and  circular.  A  '"'rizat.on  ellipse  is  used  to  establish  the  various 
polarization  states.  For  the  polarization  ellipse  shown  in  Figure  21  the  Axial 
Ratio  (AR)  and  tilt  angle  t  are  used  to  completely  describe  the  polarization  states 
of  E,  [13].  For  Figure  21  the  following  definitions  apply: 

OA  =  Semi-major  ellipse  axis 

OB  =  Semi-minor  ellipse  axis 

r  =  Ellipse  tilt  angle,  0  <  r  <  it 

AR  =  Ellipse  axial  ratio,  1  <  AR  =  OA/’OB  <  00 


Figure  21.  Reference  E-Field  Polarization  Ellipse 
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For  the  following  development  the  wave  propagation  direction  is  defined  as  r 

A  A  A 

where  r  =  0  x  4>  in  Figure  21.  The  orientation  and  rotational  behavior  of  Er 
defines  the  polarization  state.  Given  that  AR  and  r  are  specified,  Eqs  (136) 
and  (137)  completely  define  the  complex  unit  polarization  vector  e.  The  phase 
term  S  in  the  final  e  expression  determines  the  complex  nature  of  all  possible 
polarization  states.  For  all  elliptical  polarization  states,  including  the  special  case  of 
circular  polarization,  the  rotation  direction  of  E,  is  determined  by  sign  of  AR.  As 
viewed  in  the  propagation  direction,  an  AR  >  0  produces  left-hand  (counter¬ 
clockwise)  rotation  and  an  AR  <  0  produces  right-hand  (clockwise)  rotation. 

e  =  cosy  6  +  el1  siny  4>  (13f>) 

e  =  cot'l(,4R)  =  tan~’(— ) 

y  =  -  cos'1  (cos2e  cos2t)  (137) 

2 

6  =t3n-‘|i“2i) 

{  sin2t  ) 

Linear  polarization  states  are  achieved  by  letting  AR  approach  infinity  in  the  limit 
with  f  being  used  to  control  the  relative  orientation  of  Er.  For  an  AR  approaching 
infinity  the  definitions  m  Eq  (137)  can  be  used  to  show  that  6  approaches  zero 
while  y  approaches  the  value  of  r.  Suostitutmg  these  conditions  into  Eq  (136)  results 
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in  the  e  expression  of  Eq  (138).  This  final  expression  represents  a  time  independent 
position  vector  m  the  6 -t}>  plane  with  vector  component  magnitudes  determined  by 
the  value  of  r. 


e  =  cost  0  +  sinT  <f>  (138) 

ror  the  ray-trace  receive  formulation  being  implemented  two  reference  E-Field 
polarization  variations  are  considered.  For  the  first  case  the  reference  E-Field 
polarization  vector  e,  remains  constant  with  respect  to  the  absolute  coordinate  system 
regardless  of  the  reference  E-Field  location.  This  case  accounts  for  realistic  flight 
conditions  where  the  relative  polarization  of  a  fixed  source  remains  unchanged  for 
changes  is  aircraft  elevation  and  azimuth.  The  second  case  allows  the  polarization 
vector  e,  to  vaiy  as  a  function  of  the  reference  E-Field  location.  This  case  duplicates 
typical  test  range  conditions  where  source  and  observer  locations  are  varied  in  a 
circular  arc  with  respect  to  each  other. 

For  both  polarization  variations  the  reference  E-Field  polarization  vector  e,  is 
initially  defined  in  the  x-y  plane  with  a  -z  propagation  direction.  Under  these 
conditions,  the  complex  polarization  vector  definition  of  Eq  (136)  may  be  applied  to 
obtain  the  initial  polarization  vector  expression  for  e,  as  given  by  Eq  (139)  where  x 
and  y  are  unit  vectors  in  the  defined  absolute  coordinate  system.  For  the  constant 
polarization  case  Eq  (139)  represents  the  final  e,  expression  used  for  analysis  and 
modeling  validation.  For  the  varying  polarization  case,  unprimed  reference  E-Field 
coordinates  and  location  angles  0/(1),  (calculated  in  accordance  with  Section  6.2)  are 
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used  in  conjunction  with  the  rectangular  coordinate  rotation  expressions  of  Eq  (B.l) 
to  spatially  rotate  polarization  vector  components  e\  and  e'y  to  obtain  the  final  e, 
reference  E-Field  polarization  vector. 


e,  =  e^sinyx  +  casyy 

r  *  r  « 

=  e,x  +  eyy 


(139) 
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VII.  Monopulse  Processing 


7.1  Introduction 

This  chapter  presents  a  method  for  analyzing  monopulse  boresight  errors.  The 
development  begins  by  considering  the  E-Fields  incident  on  each  element  of  the 
phased  array  and  using  the  expressions  developed  in  Chapter  V  to  determine 
Co-Polarized  (CP)  and  Cross-Polarized  (XP)  incident  field  components  relative  to 
both  the  reference  E-Field  and  aperture  element  polarizations.  Each  polarization 
component  is  independently  processed  and  analyzed  such  that  radome  depolarization 
effects  are  distinguishable  from  antenna  element  polarization  sensitivity  effects. 
Boresight  errot  estimates  in  two  orthogonal  scan  planes  are  generated  from  complex 
mcnopulse  ratios.  Dual  plane  boresight  error  analysis  is  accomplished  in  both 
principal  and  diagonal  scan  planes. 

7.2  Co/Cross-Polarized  E-Field  Components 

Co/Cross-Polarized  field  components  3re  determined  relative  to  a  reference 
polarization  direction.  For  analysis  purposes  two  specific  reference  polarization 
directions  are  considered,  the  reference  E-Field  polarization  direction  e,  and  the 
aperture  element  polarization  direction  eA.  Depolarization  effects  due  solely  to 
propagation  through  the  radome  are  determined  by  ?t.  The  degree  of  radome 
depolarization  is  obtained  by  comparing  e,  with  the  polarization  of  the  E-Field 
incident  on  an  array  element. 

The  total  incident  E-Field  ET(m)  on  the  m"  array  element  is  given  by  Eq  (140) 
below  where  E°(m)  and  E”(m)  are  the  incident  E-Fields  calculated  in  Chapter  V  for 
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" Direct "  and  "Reflected"  ray  paths,  respectively.  Summation  of  the  n  reflected  E-Field 
terms  accounts  for  multiple  primary  reflected  ray  paths  for  the  m'k  element. 


ET(m)  =  E'(m)  *  £ 
i- 1 

=  Er(m)  eT(m) 


(140) 


£ T(m)  =  E  T(m)  T(m)  *  er  +  e  T(m)  •  dr  j 
n  1 

=  '  £f(»)  *«,  +  £  *«,[«, 

r-i  j 

*  |f®(m)  •  k,  +  £  £*(m j)  •  h,  |  A, 
=  E~P(m)  et  *  E^,(m)ftr 


(141) 


A  A  A 

Given  e,  and  propagation  direction  k  „  an  orthogonal  polarization  vector  /i(may 

A  /»  A  -A  A 

be  defined  as  h ,  =  k ,  x  er  Unit  vectors  e.  and  h ,  are  used  to  decompose  the  total 
incident  E-Field  of  Eq  (140)  into  CP  and  XP  components  by  a  simple  dot  product 
operation.  The  incident  E-Field  is  decomposed  into  orthogonal  components  as  shown 
in  Eq  (141).  Radome  depolarization  will  manifest  itself  by  introducing  an  XP 
response  as  indicated. 

The  primary  goal  is  to  characterize  radome  depolarization  effects  on  the  overall 
system,  specifically  the  effects  on  buresight  error.  For  this  development  the  element 
polarization  direction  eA  is  the  primary  vector  used  tor  CP/XP  field  decomposition 
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where  it  is  assumed  that  eA  represents  the  polarization  of  a  "typical"  array  element, 
neglecting  any  polarization  differences  due  to  the  element’s  location  within  the  array. 
Typical  planar  array  elements  include  dipoles,  open  ended  waveguides,  slotted 
waveguides,  and  horns  [1].  For  planar  arrays  constructed  with  these  elements  the 

A  A  A  „  A 

orthogonal  polarization  vector  h  A  is  defined  as  hA  =  ;iA  x  eA  where  nA  is  the 
aperture  piane  (i.e.,  plane  containing  array  elements)  surface  normal  vector. 
Aperture  polarization  vectors  eA  and  h  A  are  used  to  decompose  the  incident  E-Field 
expression  of  Eq  (140)  into  CP  and  XP  components  as  shown  below. 


£r(m)  =  ET(m)  T(m)  •  eA  *  eT(m)  •  A,) 

=  +  E  -aU 

h\  J 

+  %  (m)  •  ft,  ♦  £  £?(« ./) *  V  ^ 
>•1 

~  ECF{m)eA  +  Exp(m) 


(142) 


7.3  Element  Voltages 

By  decomposing  the  incident  E-Field  into  the  CP  and  XP  components  identified 
in  Eq  (142)  element  polarization  effects  may  be  accounted  for  m  determining 
element  voltage  responses.  Figure  22  represents  the  basic  element  model  used  for 
deriving  voltage  expressions  for  the  array  element. 
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Figure  22.  Aperture  Element  Model 

The  Aln  and  <j)m  terms  in  Figure  22  represent  the  amplitude  and  phase  weights 
used  for  controlling  the  shape  and  main  beam  pointing  direction.  Assuming  mutual 
coupling  effects  are  identical  for  all  elements,  regardless  of  their  location  within  the 
array,  the  CP  element  pattern  f^p(6,(j>),  XP  element  pattern  f*r(0,<J>),  and  element 
polarization  directions  ?A//iA  aie  identical.  This  is  a  reasonable  assumption  for 
analyzing  densely  populated  arrays  which  utilized  amplitude  tapers  for  pattern 
control.  Array  elements  near  the  aperture  edges  exhibit  atypical  mutual  coupling 
effects  which  are  offset  by  amplitude  weights  which  are  significantly  less  than 
elements  located  near  the  center.  Independent  CP  and  XP  element  patterns  allow 
for  varying  polarization  responses  to  be  analyzed  depending  on  the  specific  element 
being  used.  For  a  typical  element  as  shown  in  Figure  22  the  total  voltage  response 
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V„  can  be  expressed  as  shown  in  Eq  (143)  where  EjP(m)  and  EjP(m)  are  as  given 
in  Eq  (142). 

ym  =  K  A  £c7?(m)  */"(8,< »  El(m) } 

*  V„(m)  *  V„(jn) 

7.4  Monopulse  Error  Signal 

Monopulse  sum  and  difference  voltages  are  derived  by  applying  the  principle  of 
superposition.  The  aperture’s  far-field  radiation  patterns  are  obtained  by  coherently 
combining  the  individual  element  responses  as  given  by  Eq  (143).  As  identified  in 
Section  2.3,  a  "simple"  method  for  obtaining  monopulse  operation  is  implemented  as 
shown  in  Figure  5  [1].  The  system  depicted  in  Figure  5  is  for  single  plane  monopulse 
processing,  i.e.,  target/source  location  estimated  with  regard  to  a  single  scan  plane. 
For  this  system  the  monopulse  sum  voltage  VSUM  may  be  expressed  as  given  in 
Eq  (144)  where  the  sum  voltage  is  calculated  by  directly  summing  the  individual 
element  responses  from  all  elements  within  the  array.  The  monopulse  difference 
voltage  VDEL  is  obtained  by  dividing  the  aperture  into  two  symmetrical  half-planes  for 
which  the  element  responses  within  each  half-plane  are  summed.  The  half-plane  sum 
voltages  are  then  phase  weighted  and  summed  such  that  a  180’  phase  change  occurs 
in  VDEL  upon  crossing  from  one  half-plane  to  the  other.  Half-plane  locations  are 
arbitrary  provided  elements  are  symmetrically  located  on  either  side  of  the  dividing 
line. 
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^ SUM  ~  E  (k„(«)  ♦  v„m) 

« ■! 

/ 

=  E  +  //P(0 , 4>)  }  ] 


A  dual  plane/four  quadrant  monopulse  system  is  used  for  analysis  and  modeling 
purposes.  Typical  systems  employ  monopulse  operation  in  two  orthogonal  planes, 
i.e.,  azimuth  and  elevation  or  equivalently,  yaw  and  pitch.  For  these  systems  the 
aperture  is  divided  into  symmetrical  quadrants  and  a  monopulse  difference  voltage 
VDEL  is  generated  for  each  "plane"  of  monopulse  operation.  Element  voltages  within 
each  aperture  quadrant  are  summed  to  produce  quadrant  sum  voltages.  Quadrant 
sums  are  then  combined  to  produce  a  single  monopulse  sum  voltage  and  two 
independent  difference  voltages,  one  for  each  monopulse  plane  of  operation.  Sum 
and  difference  voltage  expressions  for  analysis  and  modeling  are  derived  by  assuming 
the  aperture  is  initially  located  in  the  x-y  plane.  As  shown  in  Figure  23  the  aperture 
is  divided  into  four  symmetrical  quadrants  A,  B,  C,  and  D.  Element  voltages  of  the 
form  given  by  Eq  (143)  are  summed  within  each  quadrant  to  produce  element 
quadrant  voltages  VA,  VB,  Vc,  and  VD.  For  element  coordinates  given  by 
[xt(m),y,(m),0]  the  resulting  quadrant  voltages  may  be  expressed  as  given  by  Eq  (145) 
From  the  quadrant  voltages  of  Eq  (145)  monopulse  sum  and  difference  voltages 
may  be  expressed  as  given  by  Eq  (146)  where  the  x-z  and  y-z  planes  correspond  to 
the  azimuth  (AZ)  and  elevation  (EL)  planes,  respectively.  The  180'  difference 
pattern  phase  shift  required  across  the  aperture  is  determined  by  the  phase 
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Figure  23.  Four  Quadrant  Monopulse  Geometry 


=  T,  {FCJ,(m)  +  for  xfm)  <0  and  ys(m)  >  0 

W 

VB  =  Y,  {VcAm)  *  F„,(m)}  for  xfm)  >  0  and  yt(m)  >  0 

m 

Vc  =  £  {Fcf.(m)  *  ^(m)]  for  xfm)  <  0  and  y,(m)  <  0 

M 


F0  =  £  (Kc/>(m)  +  for  xt(m)  >  0  and  >’c(m)  <  0 


Vsuu  -  +  VB  *  Vc  *  VD 

v£l  =  K  *  Vc)**'  +  +  Vo)*-'*  (146) 

V&L' -  K  +  Vm)*"  +  (*t  +  Vo)e',( 

terais  a  and  /3.  The  values  of  a  and  0  selected  determine  the  "nominal"  value  of  the 
relative  phase  difference  between  VSUM  and  VDEL.  The  "noininal”  phase  difference 
value  represents  an  ideal  system  and  is  typically  either  O’  or  90’.  The  nominal 
relative  phase  difference  is  achieved  under  ideal  conditions  with  plane  wave 
illumination  and  no  noise,  clutter,  or  other  multi-path  effects  present.  Since  it  is 
generally  impossible  to  predict  specific  deviations  of  the  relative  phase  from  its 
nominal  value,  radars  are  typically  designed  for  one  of  the  two  nominal  relative 
values  [1],  For  this  development  a  nominal  relative  phase  value  of  O'  is  established 
by  setting  a  =  0  =  jr/2. 

For  generating  the  mooopulse  boresight  error  signal  Er?(y,v)  from  the  complex 
monopulse  voltage  ratio  VDEL/V sl!M  an  "exact '  monopulse  processor  is  implemented. 
An  "exact"  monopulse  processor  is  arbitrarily  defined  as  a  processor  capable  of 
producing  the  real  part  of  the  complex  VDEL/VSUM  ratio  perfectly  for  each  angle 
coordinate.  For  a  nominal  relative  phase  value  of  O’  the  usual  practice  is  to  process 
only  the  real  part  ot  the  complex  ratio  'The  rationale  is  that  the  target  contributes 
only  to  the  real  part  while  noise,  interference,  and  clutter  contribute  equally  to  the 
real  and  imaginary  parts"  [!].  The  "exact'  monopulse  processor  boresight 
error  E„p(y,v)  may  be  expressed  as  given  by  Eq  (147)  where  the  eiror  is  defined  over 
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the  angular  range  of  y  at  a  frequency  of  v.  The  monopulse  sensitivity  K  is 
determined  by  the  slope  of  the  normalized  difference  pattern  and  has  units  of 
(v/v)/Rad  [2j. 


T^(Y,u)  *  Re 


^DEl.1  V SUM 


=  fC'1 


where  4>0  =  Phase  of  VDEL 


cos(({)D  -  <j>5)  (Radians) 

|  "suu\ 

and  <*>,  =  Phase  of  Vsl,u 


(147) 


The  expression  in  Eq  (147)  is  referred  to  as  an  "error  signal”  because  it 
represents  the  angular  off-set  of  a  source/target  relative  to  the  monopulse  boresight 
axis.  Under  ideal  conditions  a  far-field  source  located  on  the  monopulse  boresight 
axis  produces  a  difference  pattern  voltage  of  zero  and  a  maximum  sum  pattern 
voltage,  resulting  in  Eop(y,v)  going  to  zero.  A  source/target  not  located  on  boresight 
produces  a  bipolar  error  signal  with  magnitude  indicating  the  amount  and  sign 
indicating  the  direction  of  error. 

7.5  System  Boresight  Error 

The  monopulsc  boresight  error  E„p(y,v)  is  used  to  establish  and  characterize  the 
boresight  error  for  the  overall  radcme,  radar,  and  monopulse  processing  "system". 
In  this  context,  "system"  boresight  error  (BSE)  is  defined  as  the  angle  indicated  by 
Eq  (147)  when  the  aperture  has  no  pointing  error,  t.e.,  the  aperture  scan  direction 
equals  the  source  location  This  is  equivalent  to  fixing  the  aperture  scan  direction 
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while  repositioning  the  source  until  Eq  (147)  equals  zero.  The  angular  difference 
between  the  scan  direction  and  source  location  is  the  system  BSE.  As  defined,  the 
system  BSE  is  the  key  metric  used  m  characterizing  monopulse  tracking  errors  ir.  the 
presence  of  radome  depolarization  and  phase  front  distortion.  For  a  given  aperture 
scan  direction  the  system  BSE  may  be  expressed  as  in  Eq  (148)  where  it,  is  a  unit 
vector  in  the  aperture  scan  direction  and  unit  vector  n.p  represents  the  direction  of 
the  source  location  such  that  E„p(y,v)  =  0.  Hence,  one  has 

es se  =  cos'‘(v'v)  (148) 
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VIII.  Analysis  and  Modeling  Validation 

8.1  Introduction 

This  chapter  presents  the  procedures  for  validating  results  of  Chapters  HI 
thru  VII.  An  overview  of  the  Fortran  computer  model  is  presented  to  summarize 
approximately  10,000  lines  of  computer  code  written  to  support  the  validation 
process.  A  model  flow  diagram,  input  parameters,  and  output  data  format  are 
presented. 

Modeled  results  are  validated  using  limiting  case  data,  published  empirical  and 
experimental  data,  and  production  "system"  boresight  error  (BSE)  acceptance  test 
data.  Limiting  cases  (zero  BSE)  considered  include:  air  radomes,  hemispheric 
radomes  (over  limited  scan  ranges),  and  radome  boresight  scanning.  Published 
empirical  and  experimental  data  for  a  hemispheric  radome  with  a  displaced  aperture 
is  considered  for  initial  model  validation  with  measured  BSE  data.  Measured  BSE 
data  from  multipie  production  "system"  acceptance  tests  is  then  used  for  final  model 
validation.  The  validated  production  "system”  model  is  used  to  characterize  overall 
refractive  effects  on  BSE  prediction.  Effects  of  refraction  are  characterized  in  both 
principal  and  diagonal  scan  planes. 

8.2  Model  Development 

This  section  provides  an  overview  of  the  model  developed  for  validating  analysis 
results  obtained  in  Chapters  III  thru  VII.  Detailed  coding  procedures  are  not 
presented,  rather  basic  model  control  and  functional  flow  are  considered.  The 
developed  model  predicts  the  entire  "front-to-back"  performance  of  the  overall 
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"system".  Program  control  and  flow  is  dictated  by  the  key  "system"  component 
parameters  identified  in  previous  chapters.  The  model  was  developed  on  the 
VAX  8650  processing  system  and  is  comprised  of  one  mam  program  and 
approximately  40  subroutines,  all  written  in  a  standard  Fortran  computer  language. 
With  the  exception  of  two  1MSL  math/library  subroutines  used  for  solving  reflection 
point  equations  of  Chapter  IV,  the  model  is  self-supporting. 

Figure  24  is  a  functional  flow  diagram  for  the  overall  model.  Model  input  data 
is  provided  in  three  separate  data  files,  one  data  file  per  "system"  component.  The 
RADOME.1N  data  file  provides  ogive  reference  surface  parameters  and  radome 
layer  electrical  parameters.  In  conjunction  with  a  taper  function  subroutine,  this  data 
file  completely  describes  the  radome  structure.  The  RADAR.IN  data  file  provides 
control  variables  required  to  completely  establish  aperture  and  monopulse  processing 
characteristics.  Parameters  controlled  by  the  RADAR.IN  data  file 
include:  mechanical/electrical  scanning,  element  field  pattern  and  polarization, 
element  location,  aperture  amplitude  taper,  gimbal  arm  length  and  pivot  point 
location,  azimuth/elevation  monopulse  plane,  and  co/cross-polarized  system  Boresight 
Error  (BSE)  calculation.  Operation  of  the  overall  "system"  is  controlled  by  variables 
contained  in  the  SYSTEM.IN  data  file.  This  data  file  controls  propagation 
characteristics  and  aperture  scan  angles.  Propagation  characteristics  are  controlled 
such  that  the  model  generates  output  data  with  or  without  a  radome  present,  with  or 
without  reflection  points  included,  and  with  or  without  overall  refractive  effects. 
Aperture  scan  angle  control  variables  allow  for  three  aperture  scan  patterns  to  be 
implemented,  including:  1)  a  "diagonal"  scan  through  the  origin  in  an  arbitrary 


Figure  24.  Model  Flow  Diagram 
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4>,  =  constant  plane,  2)  an  "azimuth"  scan  where  6.  and  <j>,  are  varied  to  produce  a 
constant  elevation  angle,  and  3)  an  "elevation"  scan  where  0,  and  <J>,  are  varied  to 
produce  a  constant  azimuth  angle.  The  reference  E-Field  is  also  established  and 
characterized  by  variables  and  parameters  contained  in  the  SYSTEM.IN  data  file. 
Reference  E-Field  location,  axial  ratio,  tilt,  magnitude,  phase,  and  frequency  are  all 
determined  in  accordance  with  specified  input  data. 

The  flow  diagram  of  Figure  24  indicates  the  final  output  of  the  model  is  BSE 
data.  As  defined  in  Chapter  VII ,  system  BSE  is  the  key  metric  for  characterizing  the 
monopulse  tracking  error  of  the  "system".  An  output  BSE  estimate/prediction  is 
provided  for  each  combination  of  6,  and  <J),  generated  by  the  model.  Aperture  scan 
angles  0,  and  cf>,  are  varied  within  the  model  according  to  the  flow  diagram  Upon 
completing  all  0^4>i  combinations,  as  dictated  by  the  input  data,  a  table  of  scan  angles 
and  corresponding  BSE  estimates  is  output. 

User  selectable  intermediate  outputs  are  available  from  the  model  at  various 
points  along  the  processing  path.  These  outputs  allow  for  "system"  component 
responses  to  be  independently  analyzed  and  are  identified  as  the  A*,  B*,  and  C* 
points  in  Figure  24.  At  A*  co/cross-polanzed  incident  E-Field  components  are 
available.  This  output  contains  phase  front  distortion  and  depolarization  effects  due 
solely  to  propagation  through  the  radome.  The  level  of  radome  depolarization  may 
be  ascertained  by  comparing  the  incident  E-Fields  at  this  point  with  an  established 
polarization  reference.  Voit3ge  outputs  at  B*  are  used  primarily  for  characterizing 
the  "system’s"  far-field  radiation  patterns  and  monopulse  tracking  sensitivity.  Radar 
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performance  (half-power  beamwidth,  first  sideiobe  level,  etc.)  is  characterized  at  this 
point  by  generating  model  results  with  no  radome  present.  The  final  intermediate 
output  point  C*  allows  for  the  monopulse  error  signal  E„p(y,v)  to  be  analyzed  under 
varying  conditions  presented  by  the  radome  and  radar.  This  output  corresponds  to 
the  monopulse  "S-curve"  used  for  pointing  error  corrections  in  monopulse  tracking 
systems. 

A  series  of  comparative  analysis  data  is  presented  to  characterize  and  validate  the 
model’s  performance.  Limiting  cases  for  which  BSE  =  0*  are  first  considered  as 
a  means  for  establishing  bounds  on  numerical  processing  and  modeling  errors. 
Empirically  derived  and  published  expei  imental  results  are  used  for  limiting  case  data 
validations.  Finally,  measured  BSE  data  from  a  radome  production  facility  is  used 
to  complete  the  validation  process.  Refraction  effects  for  specified  radome  structures 
are  also  investigated. 

8.3  Limiting  Case  Validation 

The  "system"  analysis  results  and  computer  model  predictions  are  first  validated 
by  comparison  with  limiting  case  data,  i  e ,  special  test  cases  for  which  system  BSE 
is  known  to  be  zero.  Three  specific  such  cases  were  considered  under  the  iesearch 
effort.  Fust,  an  arbitrarily  thick  single  layer  radome  with  both  the  relative 
permittivity  e,  and  relative  permeability  n,  set  to  unity  will  result  in  a  system  BSE  of 
zero.  This  condition  effectively  places  an  "air”  radome  into  the  propagation  pa.h 
while  introducing  no  depolarization  or  phase  front  distortion  effects.  The  second 
limiting  case  considered  is  a  hemispheric  radome.  Regardless  of  the  number  of 
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radome  layers,  layer  thicknesses,  and  layer  electrical  properties,  a  "hemispheric 
radome  will  have  zero  boresight  error  if  the  receiving  aperture  is  located  at  the 
sphere  center"  [6],  Regardless  of  where  the  aperture  is  scanned  within  the 
hemisphere,  the  symmetrically  distorted  phase  front  incident  on  the  aperture 
produces  zero  BSE.  This  is  also  true  for  arbitrarily  shaped  radomes  when  the 
aperture  is  "looking"/scanned  through  the  radome’s  boresight  axis  (tip)  if  the  radome 
is  circularly  symmetric  which  is  the  case  for  all  cases  considered  m  this  effort.  All 
limiting  case  BSE  data  is  generated  using  a  mechanically  scanned  32-element 
5-wavelength  diameter  aperture  with  dimensions  as  shown  below. 
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Figure  25.  32-Element  Quadrant  Symmetric  Aperture 


119 


The  aperture  in  Figure  25  is  quadrant  symmetric  with  uniform  amplitude  weights 
applied  to  the  elements.  This  particular  aperture  configuration  is  of  interest  because 
1)  limiting  case  BSE  conditions  are  independent  of  aperture  and  monopulse 
processing  characteristics  and  2)  experimental  data  for  this  configuration  is  available. 
Radar  parameters  were  input  into  the  model  and  the  model  run  for  the  case  with  no 
radome  present.  The  normalized  sum  and  difference  patterns  of  Figures  26  and  27 
were  returned  from  the  model.  Data  plotted  in  Figure  26  reveals  a  Half-Power 
Beamwidth  (HPBW)  of  approximately  .2027  Rads  and  a  First  Sidelobe-Level  (FSLL) 
of  approximately  -15.0  dB  for  both  the  co-polarized  and  cross-polarized  field  pattern 
responses.  The  cross-polarized  response  is  obtained  by  setting  the  reference  E-Field 
polarization  orthogonal  to  the  element  polarization  without  changing  the  reference 
E-Field  scan  plane.  The  decreased  cross-polarized  response  at  higher  scan  angles 
results  from  polarization  mismatch  conditions.  The  HPBW  and  FSLLs  generated  by 
the  model  compare  very  well  with  experimental  and  empirical  data.  From  Figure  25 
it  is  observed  that  the  32-element  aperture  is  neither  truly  rectangular  nor  truly 
circular.  A  5-wave!ength  rectangular  aperture  with  uniform  amplitude  illumination 
will  exhibit  a  HPBW  and  FSLL  on  the  order  of  .1766  Rads  and  -13.26  dB, 
respectively.  An  equivalent  circular  aperture  has  a  HPBW  and  FSLL  on  the  order 
of  .2038  Rads  and  -17.6  dB,  respectively  [15],  The  .2027  Rad  HPBW  and  -15  dB 
FSLL  obtained  from  the  model  fall  well  within  these  limits  and  are  acceptable  for 
validating  the  model  for  cases  with  a  radome  present. 
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Azimuth  Angle  (Rads) 

Figure  26.  Co/Cross-Polarized  Normalized  Sum  Patterns 


The  normalized  difference  pattern  plotted  in  Figure  27  characterizes  the 
sensitivity  of  the  monopulse  system  and  represents  both  the  co-polarized  and  cross- 
polarized  responses  which  are  identical  for  the  case  considered  The  monopulse 
sensitivity  K  of  the  system  is  defined  as  the  slope  of  the  normalized  difference  pattern 
near  boresight,  i.e.,  at  an  azimuth  angle  of  zero.  From  the  data  used  to  plot 
Figure  27  a  value  of  K  -  .1158  (v/v)/Deg  was  obtained.  This  compares  very 
favorably  with  a  published  sensitivity  value  of  K  =  .11  (v/v)/Deg  [2],  With  radar  far- 
field  patterns  and  monopulse  sensitivity  accurately  characterized,  model  validation 
using  limiting  case  BSE  conditions  is  next  considered. 
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Azimuth  Angle  (Rads) 

Figure  27.  32-Element  Aperture  Normalized  Difference  Pattern 


8.3.1  Air  Radomes.  Various  combinations  of  air  radome  parameters  were  input 
into  the  model  using  the  32-element  radar  and  monopulse  processing  previously 
validated.  Figures  28  and  29  represent  typical  BSE  data  returned  from  the  model. 
Figure  28  represents  data  obtained  for  the  'TE"  polarization  case  using  a  constant 
thickness  radome  with  three  different  fineness  ratios.  Figure  29  is  equivalent  data  for 
the  "TM"  polarization  case.  The  'TE”  and  ’TM"  polarization  designations  aie 
determined  by  the  orientation  of  the  element  polarization  vector  2C  and  reference 
E-Field  polarization  vector  2,  relative  to  the  aperture  scan  plane.  The  designation 
'TE'  is  used  when  both  c,  and  2,  are  perpendicular  to  the  aperture  scan  plane. 
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Likewise,  the  'TM"  designation  is  used  when  both  2,  and  2,  are  parallel  to  the 
aperture  scan  plane  (this  convention  is  used  throughout  the  remainder  of  the 
dissertation). 


For  the  six  cases  shown  in  Figures  28  and  29  the  predicted  BSE  is  less 
than  .014  mRads  over  the  30  degree  scan  range  when  compared  to  a  theoretical 
value  of  zero  degrees.  Approximately  twelve  additional  test  cases  were  run  for 
varying  radome  structures,  including  multi-layer  tapered  designs  with  as  many  as  four 
layers.  For  ail  cases  considered  the  predicted  BSE  never  exceeded  the  .014  mRad. 
The  .014  mRad  r.umerical/modehng  error  is  approximately  two  orders  of  magnitude 
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below  expected  measured  BSE  values.  Typical  radome  BSE  maximum  values  fall  in 
the  range  of  2  to  8  mRads.  Hence,  the  .014  mRad  modeling  error  is  within 
acceptable  limits. 


Scan  Angle  (Degs) 

Figure  29.  Limning  Case  BSE,  Air  Radome  TM  Polarization 


S.3  2  Hemispheric  Radomes.  A  hemispheric  radome  induces  symmetrical  phase 
front  distortions  on  an  incident  plane  wave  and  as  noted  previously,  a  receiving 
aperture  located  at  the  sphere  center  will  indicate  zero  BSE.  Figures  30  and  31 
represent  BSE  data  satisfying  this  condition.  The  data  presented  is  for  a  specific 
radome  design  of  a  half-wave  constant  thickness  hemispheric  radome  smoothly 
terminated  in  a  cylindrical  section. 
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Figure  30.  Limiting  Case,  Hemispheric  Radome  TE  Polarization 


The  data  was  generated  using  ogive  reference  parameters  of  R  =  3.4X, 
W  =  6.8X,  and  a  dielectric  constant  of  t.  =  5.3.  For  both  the  TE  and  TM 
polarization  cases,  a  curve  showing  overall  refractive  effects  is  also  included.  As 
shown,  ref.active  effects  do  influence  system  BSE  but  do  not  significantly  impact  BSE 
magnitude.  Equivalent  BSE  results  were  obtained  for  both  single  layer  and  multi¬ 
layer  hemispheric  radomes  with  dielectric  constants  varying  between  4.8  and  9  3.  For 
all  hemispheric  test  cases  considered,  the  limiting  case  BSE  magnitude  was  less  than 
.06  mRads  for  ali  scan  angles.  The  erratic  behavior/fluctuation  in  BSE  estimates 
beyond  a  10'  scan  angle  is  attributable  to  an  interaction  between  radar  sidelobes  and 
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the  cylinder  section  termination.  Although  these  effects  tend  to  increase  the  average 
limiting  case  BSE,  the  .06  mRad  upper  limit  BSE  estimate  is  acceptable  for  all 
practical  purposes. 


Figure  31.  Limiting  Case,  Hemispheric  Radome  TM  Polarization 


8.3.3  Radome  Boresiglu  Scanning.  Circularly  symmetric  radome  structures  induce 
symmetrical  phase  front  distortions  on  incident  plane  waves  upon  propagating 
through  the  radome  along  the  boresight  axis,  i.e.,  through  the  tip.  This  limiting  case 
condition  parallels  the  preceding  hemispheric  radome  case  ana  occurs  regardless  of 
the  actual  radome  design,  i.e.,  is  independent  of  the  number  of  radome  layers  and 
electrical  properties.  Radome  designs  used  for  generating  data  in  Figures  28  and  29 
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were  used  to  establish  a  worst  case  BSE  error  estimate  for  an  azimuth  scan  angle 
of  O'.  Multi-layer  and  single  layer  radome  designs  with  fineness  ratios  in  the  range 
of  2  to  3  were  considered.  Radome  design  frequency,  number  of  layers,  and 
electrical  properties  were  varied  for  both  the  TE  and  TM  polarization  cases.  Given 
a  specific  number  of  radome  layers,  corresponding  taper  functions  (including  constant 
thickness),  and  a  specified  fineness  ratio,  the  system  BSE  at  0'  remained  virtually 
constant  for  all  variations  of  radome  electrical  properties,  E  -Field  polarizations,  and 
refractive  effect  combinations.  The  maximum  limiting  case  BSE  ootained  was 
primarily  dependent  on  radome  geometrical  properties  and  varied  between  .01 
and  .02  mRads  for  all  test  cases  considered.  As  with  previous  limiting  case  BSE 
bounds,  this  error  limit  is  acceptable  and  finalizes  the  analysis  and  modeling 
validation  under  limiting  case  conditions. 

8.4  Hemispheric  Radome  Validation:  Displaced  Aperture 

Initial  validation  of  analysis  and  modeling  results  with  experimentally  measured 
data  is  accomplished  using  a  "system"  comprised  of  the  32-element  quadrant 
symmetric  aperture  and  the  uniform  half-wave  hemispheric  radome  of  Section  8.3.2 
with  R  =  3.4X,  W  =  6.8A,  and  e,  =  5.3.  The  aperture  center  of  rotation/pivot 
point  is  displaced  one  wavelength  behind  the  radome  center  of  curvature, 
intentionally  introducing  a  small  measurable  amount  of  system  BSE.  Model  results 
are  validated  by  comparison  with  published  experimental  and  predicted  BSE  data. 
The  data  plotted  in  Figures  32  and  33  represents  measured  (*)  and  numerically 
predicted  (o)  results  for  ’system"  BSE  [2). 


BSE  (mRads)  «g  BSE  (inRads) 


Scan  Angle  (Degs) 


Hemispheric  Radotne,  Displaced  Aperture  TB  Polarisation 
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Figure  33.  Hemispheric  Radome,  Displaced  Aperture  TM  Polarization 
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The  dashed  line  represents  model  results  and  clearly  indicates  a  good  first-order 
approximation  to  plotted  curves.  Because  of  restrictions  placed  on  the  reflected 
E-Field  contributions  during  modei  development,  i.e.,  reflected  E-Field  equations  of 
Chapter  IV  are  valid  only  for  radomes  with  fineness  ratios  greater  than  two.  the  first- 
order  approximation  of  the  mode!  does  not  possess  the  higher-order  variations 
present  on  the  measured  and  numerically  predicted  data  (by  definition  all 
hemispheric  radomes  have  a  fineness  ratio  of  one).  Numerically  predicted  estimates 
of  the  published  data  were  obtained  via  a  computationally  intense  surface  integration 
technique.  Current  model  results,  obtained  from  a  geometric  optics  technique  which 
ranks  much  lower  in  computational  intensity,  compares  very  well  with  the  surface 
integration  results  In  fact,  data  presented  in  Figures  32  and  33  indicates  the  current 
technique  provides  a  superior  estimate  of  first-order  BSE  characteristics  (considering 
the  average  slope  of  the  data  presented). 

Consistent  with  previous  modeling  efforts,  the  correlation  between  measured  and 
predicted  data  is  better  for  the  TM  polarization  case  than  it  is  for  the  TE 
polarization  case.  These  differences  are  generally  "associated  with  the  planar  slab 
approximation  for  TE  and  TM  modes  in  the  highly  curved  radome  tip 
region"  [2,  21).  Locally  planar  approximations  are  less  valid  in  this  region  resulting 
m  poorer  transmission  coefficient  estimates.  Although  there  is  no  "highly"  curved  tip 
region  for  the  hemispheric  radome  under  consideration,  the  radius  of  curvature  in  the 
perpendicular  (TE)  plane  of  incidence  is  smaller  than  the  radius  of  curvature  in  the 
parallel  (TM)  plane  of  incidence  for  a  given  azimuth  scan  angle.  Considering  the 
terminated  hemisphere  geometry  and  displaced  aperture  location,  the  radius  of 
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curvature  decreases  for  TE  polarization  and  increases  for  TM  polarization  as  the 
azimuth  scan  angle  is  increased  This  condition  accounts  for  poorer  BSE  correlation 
at  higher  scan  angles  for  the  TE  polarization  test  case. 


8.5  Production  System  Validation 

8.5.1  Introduction.  A  production  radome.  radar,  and  monopulse  processing 
"system"  is  considered  for  final  analysis  and  modeling  validation.  The  production 
"system"  is  used  for  model  validation  for  several  reasons.  First,  measured  electrical 
and  mechanical  data  is  available  on  each  component  of  the  "system”.  Given 
measured  component  data,  intermediate  model  outputs  may  be  used  to  validate 
individual  component  responses  prior  to  validating  the  overall  model  response. 
Second,  acceptance  test/quality  control  data  from  several  production  units  is  available. 
Validated  analysis  and  modeling  results  are  obtained  by  comparison  with  average 
data  from  multiple  test  cases,  eliminating  the  possibility  of  validating  against  an 
'isolated"  test  case  which  "happens"  to  match  model  results.  Third,  acceptance  test 
results  on  production  units  are  well  documented  and  accomplished  under  tightly 
controlled  test  conditions.  Therefore,  model  inputs  may  be  established  such  that 
modeled  conditions  accurately  match  conditions  used  in  taking  measured  data. 
Lastly,  measured  system  BSE  data  is  available  on  the  overall  integrated  "system".  A 
final  acceptance  test  accomplished  by  radome  production  facilities  involves 
characterizing  the  system's  far-field  radiation  patterns  and  BSE  response.  A 
complete  radome.  radar,  and  monopulse  processing  system  is  utilized  for  these  tests. 
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Final  analysis  and  modeling  validation  is  accomplished  by  comparing  model  results 
with  the  final  system  BSE  measurements. 

The  following  sections  describe  the  validation  process  used.  As  with  previous 
validation  cases,  the  radar  aperture  and  monopulse  processing  portions  of  the  model 
are  completely  characterized  and  validated  prior  to  including  a  radome  in  the  system. 
All  component  specifications  and  measured  data  used  in  the  following  validation 
process  were  provided  by  a  radome  production  facility.  The  "system"  used  is 
comprised  of  a  mechanically  scanned  1368-element  aperture  corporately  fed  to 
produce  dual-plane  monopulse  processing  and  a  single  layer  tapered  radome. 

8.5.2  Modeled  Radar  Aperture.  A  mechanically  scanned  1368-element  planar 
array  is  used  for  validating  the  radar  and  monopulse  processing  portion  of  the  model. 
The  aperture  is  approximately  28-wavelengths  in  diameter  with  array  elements 
oriented  in  a  triangular  lattice  grid  pattern.  The  elements  are  slotted  waveguides 
possessing  negligible  cross-polarized  responses.  A  corporate  RF  feed  structure  is 
used  to  simultaneously  produce  a  single  sum  pattern  output  and  independent  azimuth 
(x-z  plane)  and  elevation  (y-z  plane)  difference  pattern  outputs.  The  measured  (*) 
data  plotted  in  Figure  34  represents  the  normalized  radial  amplitude  taper  applied 
across  the  aperture.  By  comparing  the  raised  cosine  and  cosine2  amplitude  tapers 
with  the  measured  taper,  it  is  evident  that  the  measured  taper  approximately  equals 
the  cosine  taper  for  radial  distances  less  than  one-third.  Beyond  this  point,  the 
measured  taper  is  approximately  the  average  value  of  the  cosine  and  cosine2  tapers 


These  conditions  define  the  amplitude  taper  used  for  model  validation  and  are 
summarized  by  Eq  (149)  where  x  is  the  normalized  radial  distance  from  the  aperture 
center  [22],  Equation  (149)  is  plotted  in  Figure  35  for  comparison  with  the 
measured  taper  using  a  pedestal  height  of  h  =  .185  .  Taper  values  calculated  in 
accordance  with  Eq  (149)  closely  approximate  the  measured  taper  and  are  acceptable 
for  modeling  purposes. 


Figure  34.  Amplitude  Taper,  Measured  vs.  Cosme/Cosme2 
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Figure  35.  Amplitude  Taper,  Measured  vs.  Calculated 
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The  modeled  amplitude  taper  of  Eq  (149)  is  used  in  conjunction  with  given  array 
element  locations  to  characterize  and  validate  the  radar  far-field  pattern 
characteristics  and  monopulse  processing  sensitivity.  Figures  36  and  37  show  the 
co-polarized  normalized  sum  and  difference  pattern  responses  of  the  modeled  radar 
system.  A  HPBW  of  2.488*  and  FSLL  of  -30.07  dB  were  calculated  from  model 
results  of  Figure  36.  These  radiation  parameters  compare  extremely  well  with 
specified/measured  HPBW  and  FSLL  values  of  2.5*  and  -30  dB,  respectively, 
validating  the  modeled  radar  geometry  and  amplitude  taper  function. 


133 


Normalized  Difference  Pattern  (v/v)  Normalized  Sum  Pattern  (dB) 


Figure  36.  Production  System,  Normalized  Sum  Pattern 


Figure  37.  Production  System,  Normalized  Difference  Pattern 
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The  co-polarized  normalized  difference  pattern  of  Figure  37  is  used  to 
characterize  the  radar  system’s  monopulse  sensitivity.  From  the  slope  of  the 
normalized  difference  pattern  at  an  azimuth  angle  of  0\  a  calculated  monopulse 
sensitivity  of  K  -  31.0  (v/v)/Rad  is  obtained.  This  model  value  identically  matches 
the  measured  sensitivity  value  provided  by  the  radome  production  facility.  The  near 
perfect  pattern  characteristics  and  monopulse  sensitivity  value  clearly  validate  the 
model’s  radar  response  for  the  production  system. 


Azimuth  Angle  (Rads) 


Figure  38.  Production  Aperture,  Uniform/Radial  Taper  Comparison 

The  modeled  production  system  radar  is  used  to  validate  the  model  response  for 
various  amplitude  tapers  Far-field  patterns  tsing  the  1368-element  28-wavelength 
diameter  aperture  were  generated  for  various  amplitude  tapers  provided  by  the 
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model.  Specifically,  uniform  and  radial  (approximate  cosine)  amplitude  tapers  were 
characterized.  Figure  38  is  a  plot  of  calculated  sum  patterns  comparing  a  uniform 
amplitude  taper  response  with  three  different  radial  tapers.  Table  1  provides 
comparison  data  for  calculated  HPBW  and  FSLL  values  from  data  of  Figure  38  and 
published  theoretical  data  for  circular  apertures  [17).  Calculated  HPBW  and  FSLL 
values  compare  very  well  with  corresponding  theoretical  values,  validating  mode! 
results  for  the  tapers  considered. 

Table  1.  Taper  Comparison,  Calculated, Theoretical  HPBWs  and  FSLLs 


TAPER 

Uniform 

Radial 

Radial2 

3 

Radial 

HPBW  (Dag) 

2.08/2.02 

2.44/2  52 

2.75/2.91 

3.15/3.34 

FSLL(-dB) 

17.53/17.60 

24.19/24.50 

30.48  /  30.50 

35.95/35.80 

( Calculated  /  Theoretical ) 


85.3  Modeled  Production  Radome  The  production  radome  considered  for 
validation  purposes  is  a  single  layer  tapered  radome.  The  specific  design  selected 
satisfies  the  circularly  symmetric  requirement  and  falls  within  the  moderate-to-large 
sized  class  of  radomes  (exceeds  the  5  to  10-wavelength  dimensional  requirement) 
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The  radome  is  produced  on  an  inner  contour  mandrel  by  building  up  a  layer  of 
dielectric  material  exceeding  the  desired  Outer  Mold  Line  (OML)  thickness 
specifications.  Once  properly  cured,  the  radome  OML  is  ground  to  specifications 
such  that  the  desired  radome  taper  is  achieved.  Radome  production  specifications 
and  acceptance  test  data  packages  were  obtained  on  three  identical  radome  units 
from  an  uniaentified  production  facility.  Acceptance  test  data  packages  for  the  three 
radomes  were  randomly  selected  from  approximately  40  available  production  units. 
Measured  "system"  BSE  data  used  in  the  following  validation  process  is  "average" 
data,  representing  averaged  BSE  values  from  the  three  selected  production  units. 

The  first  step  in  validating  model  results  with  the  production  radome  was  to 
determine/define  the  reference  ogive  surface  parameters  H  and  W.  Mandrel  data 
supplied  with  the  radome  specifications  was  plotted  and  compared  with  several  ogive 
surfaces  generated  by  varying  ogive  parameters  Rand  W.  Comparisons  revealed  that 
the  production  radome  inner  surface  could  be  very  closely  fitted  to  a  truncated 
tangent  ogive  surface,  i.e.,  a  non-tangent  ogive  surface  generated  by  displacing  the 
tangent  ogive  base  plane  toward  the  ogive  tip.  By  using  a  tangent  ogive  reference 
surface  with  a  compensated  aperture  gimbal  point  location,  an  accurate  inner  radome 
surface  equation  was  generated  using  values  of  R  =  234.4727  and  W  =  36.136X, 
resulting  in  an  overall  radome  length  of  L  =  90.267  The  specified  aperture  gimbal 
point  location  of  z%  =  15.933  A  was  re-established  as  ze  =  28.991a  for  the  modeled 
radome.  The  specified  gimbal  arm  length  of  3.983a  remained  unchanged  for 
modeling  purposes. 
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The  OML  grind  specifications  were  used  in  conjunction  with  established  ogive 
reference  surface  parameters  to  determine  the  layer  taper  function  t(0).  A  list  of 
radome  station  locations  (distances  from  the  radome  tip)  and  corresponding  grind 
thicknesses  were  provided  with  the  radome  production  data.  The  thicknesses 
provided  were  determined  on  the  original  radome  development  effort  by  empirical 
methods  using  countless  "build  and  test"  iterations.  As  such,  no  analytical  procedures 
were  applied  to  determine/establish  a  "function”  to  accurately  describe  the  radome 
taper.  Since  the  grind  thicknesses  were  measured  along  the  surface  normal  direction 
at  given  radome  station  locations,  the  procedure  developed  in  Section  3.4.2  for 
generating  "ideal'  taper  functions,  Eqs  (50)  thru  (55),  was  applied  to  convert  normal 
thickness  grind  specifications  into  corresponding  taper  function  values.  Specifically, 
grind  thickness  values  were  used  in  place  of  the  d,  values  calculated  by  Eq  (52). 

The  normalized  taper  function  profile  obtained  from  the  procedure  using  the 
grind  specifications  is  provided  in  Figure  39  as  "actual"  taper  values.  The  lack  of 
"smoothness"  at  scan  angles  around  5'  is  attributable  to  constant  normal  thickness 
grind  specifications  in  the  radome  tip  region  The  actual  radome  taper  as  a  function 
of  surface  normal  direction  begins  at  about  9.5k  from  the  radome  tip.  The  "ideal' 
Gaussian  taper  function  of  Eq  (56),  as  established  using  identical  ogive  surface 
parameters  and  dielectric  constant  for  a  half-wave  wall  design,  is  also  plotted  in 
Figure  39.  The  shape  of  the  two  curves  is  nearly  identical  with  a  difference  of  less 
than  .05k  in  normalized  taper  value  over  the  30'  scan  range.  In  fact,  a  shift  in  scan 
angles  of  approximately  5*  for  either  of  the  two  curves  results  in  a  near  perfect 
match.  Given  that  the  original  radome  taper  design  was  empirically  "tuned”  for 
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optimum  transmission  and  BSE  performance,  the  near  perfect  fit  of  "actual"  and 
"ideal'  curves  in  Figure  39  represents  an  important  step  in  validating  the  "ideal'  taper 
function  concept  developed  m  Section  3.4.2. 


Figure  39.  Production  Radonte,  Actual  vs.  "Ideal”  Taper 


8.5.4  Final  BSE  Validation.  The  actual  and  "ideal'  taper  functions  are  both  used 
m  the  final  BSE  validation  procedure  Acceptance  test  BSE  data  was  provided  at 
two  distinct  test  frequencies.  The  measured  (»)  data  points  in  Figures  40  and  41 
represent  averaged  BSE  values  at  frequencies  of  FD  and  Fih  respectively.  These 
i  frequencies  represent  the  radome  design  frequency  (F0)  and  the  upper  edge  limit  of 

I  an  8.5%  bandwidth  (Fsl).  The  data  was  obtained  using  linearly  polarized  array 
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elements  (et  =  x)  aiid  a  variably  polarized  reference  E-Field  (e,  =  x).  A  nominal 
dielectric  constant  value  of  s,  =  4.8  with  a  loss-tangent  of  .014  was  assumed  for  the 
tapered  layer.  The  aperture  was  scanned  in  the  monopulse  rzimuth  plane  (x-z  plane) 
and  corresponding  in-plane  BSE  estimates  were  computed.  The  solid  and  dashed 
curves  represent  mode!  results  using  the  actual  and  "ideal'  tapers,  respectively. 
Modeled  production  "system"  results  required  approximately  20  minutes  of  computer 
processing  time  to  calculate  and  outn”?  31  discrete  scan  angle  and  BSE  estimates, 
approximately  40  seconds  pel  0,/4>,  combination.  This  equates  to  a  processing  time 
of  approximately  30  mSec/element/scan  angle  for  the  1368-element  aperture  being 
modeled.  Considering  the  goal  to  limit  computational  intensity,  this  is  an  acceptable 
processing  rate  for  developmental  system  design  efforts. 

Comparison  of  measured  BSE  data  with  modeled  results  using  the  actual  taper 
reveals  a  BSE  prediction  erior  of  approximately  .5  mRads  or  less  for  both 
frequencies  considered,  extremely  good  prediction  results  considering  the  geometric 
optics  propagation  technique  being  implemented.  For  a  moderate  decrease  in  the 
nominal  dielectric  constant  value  (approximately  6%)  a  near  perfect  match  between 
actual  taper  and  measu.ed  BSE  data  could  be  achieved.  Radome  production 
tolerances  and  variations  in  dielectric  material  electrical  properties  could  easily 
account  for  a  .5  mRad  prediction  error.  Given  that  these  results  were  obtained  for 
!)  an  arbitrarily  selected  radome  design  satisfying  dimensional  and  symmetry 
requirements,  2)  multiple  randomly  selected  production  radome  units,  and  3)  multiple 
frequency  BSE  measurements,  they  clearly  validate  the  overall  analysis  and  modeling 
approach  developed  ,n  Chapters  HI  thru  VII  for  the  production  "system"  considered. 
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"Ideal'  taper  model  results  plotted  in  Figures  40  and  41  compare  favorably  with 
actual  taper  and  measured  BSE  data.  At  FD  "ideal'  BSE  predictions  are  consistently 
lower  than  actual  mper  results.  Assuming  predicted  BSE  results  are  co  isistently 
higher  than  measured  BSE  values,  as  indicated  in  figure  40  by  comparison  of  actual 
taper  and  measured  BSE  values,  the  lower  "ideal'  BSE  predictions  suggest  that 
implementation  of  the  “ideal'  taper  function  in  an  actual  radome  structure  may 
improve  measured  BSE  performance.  Data  presented  in  Figures  39  thru  41  relating 
actual  taper  and  "ideal  taper  characteristics  strongly  suggest  that  the  "ideal  taper 
function  development  process  produces  a  candidate  taper  function  which  1)  closely 
approximates  empirically  "tuned"  tapers,  potentially  resulting  in  fewer  "build-and-test” 
iterations  during  radome  design/development,  and  2)  yields  modeled  BSE  predictions 
which  are  consistent  with  both  measured  and  empirically  tuned  results. 

8.6  Refractive  Effects  on  BSE  Prediction 

Refractive  propagation  effects  on  BSE  prediction  are  characterized  using  the 
validated  production  "system"  model.  Refractive  effects  are  incorporated  into  the 
model  in  accordance  with  procedures  developed  in  Section  5.5  and  are  characterized 
in  both  principal  and  diagonal  scan  planes  at  several  test  frequencies. 

The  azimuth  monopuhe  processing  plane  (x-z  plane)  is  used  for  principal  scan 
plane  characterization.  Figures  42  and  43  are  modeled  results  obtained  for 
frequencies  of  Fn  and  F,„  respectively.  A  marginal  reduction  in  predicted  BSE  is 
indicated  over  the  30‘  scan  range,  corresponding  to  a  marginal  improvement  in  BSE 
prediction  when  compared  with  data  in  Figures  40  and  41. 
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The  marginal  improvement  in  BSE  prediction  is  achieved  at  the  expense  of 
computational  efficiency.  Additional  convergence  requirements  and  overhead 
calculations  are  required  when  overall  refractive  effects  are  included  in  model 
calculations.  Refractive  effect  calculations  increase  computer  processing  time  by 
approximately  80%,  from  20  to  36  minutes  for  31  discrete  data  point  calculations. 
A  severe  processing  penalty  is  incurred  if  second-order  refractive  effects  are  included 
;n  BSE  calculations  when  considering  the  marginal  "improvement"  in  predicted  BSE 
values. 

Diagonal  scan  plane  characterization  is  included  for  completeness.  Field  pattern 
anomalies  in  the  diagonal  scan  plane  generally  result  in  poorer  measured  "system" 
BSE  performance  tor  given  "system"  designs.  Cross-polarized  field  components  and 
asymmetrical  phase  front  distortions  are  introduced  along  diagonal  scan  planes, 
producing  "system"  BSE  in  both  the  azimuth  and  elevation  monopulse  processing 
planes.  Data  presented  in  Figures  44  thru  46  represents  modeled  BSE  results  for  a 
diagonally  scanned  aperture  and  is  "typical"  of  data  obtained  for  both  monopulse 
processing  planes.  The  aperture  was  scanned  from  the  origin  in  a  constant 
4>,  =  tr/4  Rad  scan  plane.  Data  presented  is  for  three  distinct  frequencies 
representing  the  design  frequency  FD  and  lower/upper  band  edges  of  a  6.5% 
bandwidth,  FL  and  F„  respectively.  The  figures  clearly  indicated  marginal  difference 
between  refractive  and  non-refractive  test  cases.  As  with  the  production  "system" 
refractive  results,  a  marginal  improvement  in  BSE  prediction  is  obtained  for  a 
substantiai  sacrifice  in  computational  efficiency. 
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IX.  Conclusions  and  Recommendations 


9.1  Conclusions 

A  unique  analysis  and  modeling  approach  was  developed  to  determine  radome 
depolarization  and  phase  front  distortion  effects  on  monopulse  receiver  tracking 
performance.  Previously  developed  ray-trace  receive  propagation  techniques  were 
modified  to  include  ray  refractive  effects  in  analyzing  moderate-  to  large-sized 
radar/radome  systems.  A  novel  multi-layer  tapered  radome  model  was  developed  to 
analyze  arbitrarily  shaped  radomes  using  a  tangent  ogive  reference  surface.  The 
radome  model  development  included  a  unique  'ideal'  taper  function  concept  based 
on  maintaining  constant  electrical  thickness  for  arbitrary  aperture  scan  angles. 
Reflected  E-Field  contributions  were  incorporated  into  the  analysis  using  calculated 
reflection  points,  equations  developed  for  an  arbitrary  reflecting  surface  were  used 
to  determine  specular  (primary)  reflection  point  coordinates.  Analysis  results  were 
implemented  in  a  Fortran  computer  model  which  predicts  the  "front-to-back" 
boresight  error  (BSF.)  performance  for  an  integrated  "system."  Model  results  were 
validated  using  empirical  limiting  case  data,  published  experimental  data,  and 
production  "system"  acceptance  test  data  Predicted  BSE  values  from  the  model  were 
consistently  within  ±  1  mRad  of  published  surface  integration  results  for  hemispheric 
radome  data  and  within  ±  .5  mRad  of  measured  production  system  acceptance  data, 
"excellent"  results  by  previously  established  standards. 
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9.2  Recommendations 


Four  areas  are  recommended  for  extending  the  current  work  into  future  analysis 
and  modeling  efforts  including,  1)  using  the  developed  multi-iayer  tapered  radome 
model  and  calculated  reflection  results  with  a  surface  integration  propagation 
technique,  2)  a  polarization  sensitivity  study  to  characterize  BSE  performance  under 
varying  polarization  conditions,  3)  validation  of  the  "ideal'  taper  function  concept  with 
additional  measured  data,  and  4)  using  the  current  analysis  and  modeling  procedure 
to  predict  image  lobe  level  and  location.  Each  of  these  recommendations  are 
addressed  in  the  following  paragraphs. 

Hemispheric  data  presented  in  Figures  22  and  23  indicates  a  significant 
improvement  in  BSE  estimation  may  be  realized  if  current  analysis  and  modeling 
procedures  are  implemented  using  surface  integration  propagation  techniques. 
However,  full  implementation  of  such  a  technique  will  increase  required  processing 
time,  possibly  to  an  unacceptable  level  for  developmental  system  analysis.  A  "hybrid" 
propagation  technique  should  be  considered  to  reduce  the  processing  burden.  From 
Figures  22  and  23  it  is  evident  that  the  current  ray-tracing  development  accurately 
predicts  first-order  BSE  characteristics.  Reflected  E-Field  components  generally 
contribute  to  the  higher-order  "variations"  present  on  predicted  and  measured  BSE 
curves,  "variations"  which  are  most  accurately  predicted  by  surface  integration 
techniques.  Future  studies  should  consider  using  a  "hybrid"  propagation  technique 
to  combine  first-order  ray-tracing  and  higher-order  surface  integration  results  into  a 
composite  BSE  estimate,  i.e.,  superimpose  surface  integration  variations  on  ray-trace 
receive  data.  Excellent  BSE  prediction  results  are  obtained  if  this  procedure  is 
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applied  to  data  in  Figures  22  and  23,  indicating  that  a  "hybrid"  propagation  technique 
may  provide  improved  BSE  estimates  at  processing  speeds  ranging  between  ray- 
tracing  and  surface  integration  processing  rates.  A  theoretical  basis  for  "hybrid" 
analysis  and  modeling  should  be  established  and  validated  if  appropriate. 

Polarization  "purity"  was  maintained  throughout  analysis  and  model  development 
of  the  current  research  effort,  i.e.,  E-Field  polarization  responses  were  accurately 
accounted  for  along  propagation  paths.  As  a  result,  current  results  are  well  suited 
for  polarization  sensitivity  analysis  for  varying  polarization  conditions.  Arbitrarily 
polarized  reference  E-Fields,  arbitrarily  polarized  aperture  elements,  independent 
co/cross-polarized  processing  channels,  and  accurate  radome  characterization  are 
essential  elements  for  accurately  characterizing  radome  depolarization  and  BSE 
performance  under  varying  polarization  conditions.  Previous  polarization  studies 
indicated  that  modulation  of  the  reference  E-Field  polarization  (vary  the  axial  ratio 
or  tilt  angle)  generally  results  in  improved/degraded  BSE  performance  depending  on 
the  specific  "system"  design.  Results  of  such  analyses  would  provide  the  radar  and 
electronic  countermeasure  technical  communities  valuable  insight  into 
enhanced/degraded  "system"  performance.  A  polarization  sensitivity  study  would  also 
provide  information  on  the  level/amount  of  radome  depolarization,  as  opposed  to  the 
"effects'  of  depolanzation.  The  "effects"  of  radome  depolarization  on  various  system 
parameters  have  been  extensively  analyzed  with  few  research  efforts  attempting  to 
characterize  the  depolarization  itself,  i.e.,  establish  the  level/amount  of  polarization 
"difference"  between  incident  and  transmitted  E-Fields.  The  analysis  and  modeling 
technique  developed  under  this  research  effort  provides  intermediate  results  for 
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predicting/establishing  radome  depolarizat'on  levels.  Improved  range  measurement 
techniques  are  required  to  provide  reliable  measured  depolarization  data.  When 
available,  a  comparative  analysis  between  predicted  and  measured  data  should  be 
accomplished, 

The  "ideal'  taper  function  concept  developed  in  Section  3.4.2  and  implemented 
in  Section  8.5.3  required  additional  validation  against  measured  BSE  data. 
Specifically,  an  empirically  "tuned”  production  radome  design  should  be  modified  and 
built  to  include  an  "ideal'  taper.  Measured  data  from  the  two  radome  designs  should 
then  be  analyzed  and  compared  to  determine  the  validity  of  the  "ideal'  taper  function 
concept.  Data  presented  in  Figures  39  thru  41  suggests  the  "ideal'  taper  function 
approach  may  reduce  the  number  of  developmental  "built-and-test"  iterations 
required  to  obtain  "near"  optimum  radome  performance,  potentially  reducing  radome 
development  time  and  cost. 

As  depicted  in  Figure  4,  reflected  rays  generally  form  an  image  lobe(s)  anytime 
a  significant  amount  of  energy  is  reflected  from  the  radome’s  inner  surface.  The 
relative  level  and  location  of  the  image  lobe(s)  provides  an  indication  of  "how  well" 
the  radome  is  performing.  High  image  lobe  levels  indicated  poor  transmission 
qualities  and  image  lobe  location  provides  an  estimate  of  where  taper  improvement 
is  required.  Accurate  prediction  of  image  lobe  level  and  location  would  greatly  aid 
in  radome  development  efforts.  Considering  the  ability  to  accurately  calculate 
specular  reflection  points  ( Chapter  IV)  and  far-field  radiation  patterns,  the  current 
modeling  technique  is  ideally  suited  for  future  image  lobe  analysis.  The  current 
analysis  and  modeling  technique  should  be  extended  to  allow  the  reference  E-Field 
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location  to  vary  beyond  the  aperture  field-of-view,  allowing  far-field  radiation  patterns 
to  oe  generated  within  ±  180'  of  the  aperture  scan  direction.  Modeled  results  would 
then  be  compared  with  measured  image  lobe  data  to  validate  the  image  lobe  analysis 
procedure. 


151 


Appendix  A:  Special  Case  Reflection  Points 


This  appendix  is  provided  to  augment  materia!  prese  ited  m  Chapter  IV,  Radome 
Reflection  Points.  It  provides  the  detailed  derivation  procedures  used  in  calculating 
special  case  critical  points.  Eqs  (A.l)  thru  (A  8)  are  similar  to  expressions  found  in 
Chapter  IV  and  are  repeated  here  for  completeness  and  clarity  reasons.  Eqs  (A.l) 
and  (A.2)  represent  the  scaling  process  used  for  generating  dimensionless  distance 
and  velocity  quantities  vvhere  R  is  the  radome  generating  radius  and  vn  is  the  velocity 
of  propagation  in  the  medium. 


Defining :  x(t)  = 


,  y(*1  « 


( V 

y  — 

-  I  R 

R 


zfel 


(A.1) 


and  d  -  —  ,  b  =  —  ,  c  =  —  .  t  =  ~ 
v'  v_  v_  R 


From  Boundary  Conditions : 
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For  the  tangent  ogive  surface  under  consideration,  Eq  (A.3)  provides  the 
expression  for  the  surface  equations  and  scaled  gradient  components  at  the  critical 
point.  The  sca'ed  surface  expression  of  Eq  (A.3)  and  the  three  equations  in 
Eq  (A.4),  in  conjunction  with  the  change  of  variable  expressions  of  Eq  (A.5), 
represent  a  system  of  four  equations  and  four  unknowns,  the  unknowns  to  be  solved 
for  are  a,  p,  T„  and  X. 


Scaled  Surface  Equation  :< 

=  o  -  /T?  -  \lx2  ~  f  *K  =  0 

R 

where  K  -  —  -  '1 
2  R 


(A.3) 


From  Velocity  Conditions  • 

(■ 5  -  *4  =  ol’t) 
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a  =  since  cosp 
S  =  sino  sin (5 
c  -  cos  a 


where  0  <  P  <  2r. 

and  0<  a  <  — 
2 


As  presented  m  Chapter  IV,  Eqs  (A.6)  and  (A.7)  are  used  for  deriving  critica! 
point  solutions  for  Case  I  where  it  is  assumed  that  both  x:  c  and  y  c  are  non-zero.  Two 
subcases  emerged  with  Eq  (A.8)  representing  the  Tc  solution  for  Case  I.A,  as  derived 
m  Chapter  IV.  The  second  sub-case,  Case  LB,  is  further  sub-divided  into  a  series  of 
cases,  all  of  which  cause  the  T<  expression  of  Eq  (A.3)  to  be  invalid,  i.e.,  cause 
division  by  zero 

For  Case  I,  xc*0  and  yc*0 

A  _  (a  -  '  ><  _  [b  ~  *  y]  _  (r  -  K)'! 1  -  A  (  } 

y< 


°-K  _ B  -  K 

where  =  of 


-  X,  and  y (r()  =  bic  -  y. 


(A.7) 


-  f.sina^cosp  -£„sinp )  =  x,fsina cosp  -k})  -  _vr(suiacosp 
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For  Case  I A :  sin  a{k^  cos  p  -  Jtxsinp)  *■  0 


(A-8) 


-  *r(sina  sinp  -  -  yr(sina  cosP  -  kz) 

sin  a  (;ty  cos  p  -  Jfcxsin  p ) 

For  Case  I.B.l,  it  is  assumed  that  sin(a)  =  0.  Eqs  (A.9)  thru  (A.13)  summarize 
the  results  of  making  this  assumption.  The  goal  is  to  establish  a  set  of  conditions, 
which  when  satisfied,  result  in  sin(a)  =  0  and  provide  an  alternate  solution  to  Case 
I.A.  Eqs  (A.10)  thru  (A.12)  form  such  a  set  of  conditions,  which  when  satisfied, 
result  in  the  critical  point  coordinates  as  given  by  Eq  (A.9). 

Case  /.Bis  Suppose  sin  a  =  0 
Then  a  =  sina  cos P  -  0  —  d  =  0  —  xc  =  xt  *  0 

(A.9) 

and  6  =  sinasinP  =  0  —  6  =  0  —  yc  =  yr  *  0 
d  -  0  &  £  =  0  ~  c  =  1  ~  Zc  ~  ic  *  z. 


From  Equation  A. 6:  "Condition  #  l’ 
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From  Eons  A3  and  A  A  :•  "Condition#  2" 
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From  Eqn  A.'i  with  xc  =  xr  and  yc  =  yr:  “Condtion  #3" 
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(A.13) 


To  verify  the  condition  expressed  in  Eq  (A.l  1),  it  is  first  necessary  to  establish  a 
maximum  sailed  critical  time.  This  maximum  time  is  based  on  the  longest  distance 
bewveen  any  given  antenna  element  and  the  tip  of  the  radome.  Assuming  a 
maximum  antenna  scan  angle  of  ir/2  Rads,  a  worst  case  condition  used  only  for 
calculation  purposes,  the  maximum  distance  between  an  element  and  the  radome  tip 
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will  be  given  by  =  L  +  W/2.  Using  this  relationship,  Eq  (A.13)  shows  the 
development  process  and  final  expression  for  the  maximum  scaled  critical  time. 

For  Case  I.B.2  it  is  assumed  that  fcycos/3  =  kx sin/3,  thereby  causing  the  numerator 
of  Eq  (A.8)  to  go  to  zero.  From  Eq  (A.4)  it  can  be  shown  that  for  X  to  be  non-zero, 

A 

kz  cannot  equal  one.  Since  A  is  a  unit  vector,  this  implies  that  kx  and  ky  may  not 
simultaneously  equal  zero  and  is  therefore  an  invalid  solution  to  the  assumed 
condition.  Further,  it  can  be  shown  from  Eq  (A.4)  that  neither  the  condition  that 
kx  =  cos/3  =  0  nor  the  condition  that  ky  =  sin/3  =  0  are  possible  solutions  since  both 
would  result  in  X  equalling  zero.  Therefore,  a  solution  exists  only  if  all  four  terms  of 
the  assumed  equation  are  non-zero.  This  being  the  case,  p  may  be  explicitly  solved 
for  as  shown  in  Eq  (A.  14)  and  three  special  cases  are  considered  as  a  result  of  the 
arctangent  function. 

For  Case  I.B  2:  k}  cos  p  =  kx  sin  p 

/*  s  (A.  14) 

kx  *  0  ,  ky  *  0  ,  cosp  *  0  ,  sinP  *0  —  P=tan'‘|  — J 


For  Case  I.B.l.i :  k„  -  k„ 


„  a  -3n  r 

—  P  =  —  or  -  —  a  -  b 

4  4 
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For  Case  I.B.2.1,  kt  and  are  assumed  to  be  equal  and  Eq  (A.6)  is  used  to  show 
that  x,  must  equal  y,  as  shown  in  Eq  (A.15).  For  Case  I.B.2.ii,  fc,  is  assumed  to  be 
equal  to  ~ky  and  Eq  (A.6)  is  again  used  to  show  that  this  assumption  results  in  x, 
being  equal  to  -yr  as  indicated  in  Eq  (A.  16).  For  both  of  these  cases,  X  is  solved  for 
and  equated  to  the  right  hand  expression  of  Eq  (A.6).  The  expression  formed,  which 
when  combined  with  the  surface  constraint  expression  of  Eq  (A.3),  forms  the  system 
of  two  equations  and  two  unknowns  (T(  and  a)  shown  in  Eq  (A.17).  As  for  Case  I.A, 
the  system  generated  is  composed  of  transcendental  non-linear  expressions  requiring 
numerical  solution  techniques 


For  Case  I.B.l.ii :  ks  =  -ky 


(A.16) 


Case  LB  2  :  and  Case  /  li.2  it  Equations : 

12 


(cos a  tc  *z,f  =  • 


(cosh  -k.) 


(i)— —  -k, :  ♦  (cosi  -y.f 

v/2  J 


(A.17) 


xr>zf)2  =  l-j^|(i)^r>i!-f 

i  i  v2  i 


Take:  (+)  for  p  =  t—  ,  (-)  /or  p  =  x™ 
4  4 
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The  final  sub-case,  Case  I.B.2.iii,  is  for  arbitrary  values  of  and  ky  such  that 
kt  f  +/-  ky  and  p  is  as  given  in  Eq  (A.  14).  For  this  value  of /3,  the  left-hand  side  of 
the  final  expression  ip  Eq  (A.7)  identically  equals  zero.  As  such,  the  right-hand  side 
equals  zero  and  may  be  used  to  explicitly  solve  for  a  as  given  m  Eq  (A.18).  With  a 
and  p  both  calculated,  only  Te  remains  unknown  and  may  be  determined  using 
Eqs  (A.2)  and  (A.3).  Provided  a  7C  value  satisfying  Eq  (A.3)  exists  and  satisfies  the 
inequality  specified  in  Eq  (A.13),  a  critical  point  exists  and  its  coordinates  may  be 
calculated  via  Eqs  (A.2)  and  (A.5). 


For  Case  I.B  2.iii :  kx  *  ky  ,  xr  *  C  and  yr  *  0 


•  -if  *rky-yrkx 
a  =  sin  ^ 

JcrsinP  ->fcosp  j 


for  0  <  a  <  cos‘!(jtj) 


(A.18) 


For  Case  II:  xc  -  0  and  yc  *  0 

*c=°  “  8i| 


-x. 


H  *  £ 


=  0 


(A.19) 


d  -  kx  =  XfJ 


=  0  -  a=kx  -  xc  =  kxte  *  xr  =  0 


For  the  second  major  case,  Case  11.  it  is  assumed  that  =  0  and  yc  0. 
For  xt  -  0  the  partial  derivative  of  the  reflecting  surface  g(x.v.z)  with  respect  to  x 
identically  equals  zero.  As  a  result,  only  equality  of  the  last  two  expressions  of 
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Eq  (A.6)  remains  valid  and  the  velocity  expression  relating  a  and  A,  is  used  to 
establish  the  relationship  shown  in  Eq  (A.19).  Two  sub-cases  are  considered  as  a 
result  of  this  relationship,  Case  II A  for  A',  equal  zero  and  Case  II.B  for  A,  not  equal 
zero. 


For  Case  II A.  I  :•  kx  =  0  and  sina  =  0 
A  =  0  —  x  =  0  "Conditional" 


sina  =0  — 


b  =  0  -  yc  =  y, 


(A.20) 


(a  =  Arc  ,  k  -  0,±1,±2,.:.;. 
c  =  cos  a  =  ±  1  but  c  >  k  >  0  —  c  =  1 


From  the  Velocity  Equations 


b-ky  =  kg 


-  =  k-~-  -  ksign(yr} 


(A.21) 


1ky< 0  and  y,>0  I 

or  |  "Condition  #2" 

k  >  0  and  yr  <  0  I 


For  a  =  smacos/S  =  A,  =  0,  Case  Il.A  is  further  subdivided  into  two  special 
cases.  Case  ll.A  1  provides  soiuuons  based  on  the  assumption  that  either  sine  =  0 
or  sina  -  cos/3  =  0  and  Case  II  A  2  provides  solutions  based  on  cos/3  =  0  only. 
Eqs  (A.20)  thru  (A.22)  summarize  the  results  for  Case  II.A.l.  The  three  conditions 
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For  X  =  sign[yr)k. 
c-k:  =  Xgt 


~  z^Nj 


Lh 

2 


(A.22) 


Into  the  Scaled  Surface  Equation 


I  *k 


N  2 


-  ~  |j, |  *  K  =  0  "Condition  #3" 


are  derived  by  considering  both  the  velocity  equations  and  surface  constraints  as 
shown.  Given  all  three  conditions  are  satisfied,  a  critical  point  exists  with 
coordinates  x(  =  0,  y c  =  y „  and  ze  as  given  by  Eq  (A.22). 

For  Case  II.A.2  results  are  not  as  easily  obtained  as  outlined  in  Eqs  (A.23) 
thru  (A.25).  As  with  the  previous  case,  the  condition  that  x,  =  0  must  again  be 
satisfied.  Under  the  assumption  that  cos£  =  0  and  given  the  range  of  /3  values 
specified  in  Eq  (A.5),  only  two  /3  values  are  possible,  it  12  and  3ir/2,  as  indicated  in 
Eq  (A.23).  Solving  for  X  from  the  velocity  equations  allows  for  Tc  to  be  explicitly 
solved  for  as  a  function  of  a  as  indicated  in  Eq  (A.24). 


For  Case  II A. 2  :  kx  =  0  and  cos  p  -  0 
kx  =  0  ■»  xr  =  0  "  Condition  tf  1" 

(A.23) 

cos  p  =  0  -  sin  p  =  (±)  1  -  p  =  ±  -2 

—  b  -  sinasmP  =  (±)sina 
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V- 


.  I 

.V 


Introducing  the  function  f(a),  the  t  c  expression  is  substituted  into  the  scaled 
surface  equation  resulting  in  the  final  expression  shown  in  Eq  (A.25),  a  single 
expression  of  one  variable  a.  Provided  an  a  exists  which  satisfies  this  expression,  t  c 
can  be  calculated  from  Eq  (A.24)  and  the  critical  point  coordinates  determined 
via  Eq  (A.2). 

From  Velocity  Equations 

5-ky  =  Xgj  >♦  X  =  »*«(&)[*,**“  *] 


c-k^xsA  -re*  - - 


(k.-cosa)  ,  -  .a 

i-i - i -  =  cosar  +  z, 

12  \2  '  ‘  'I 


Xy(r)sinccf  +(X;-cosa)' 


1  f  (fc-ccsa)  .  ) 

-  0  <  tc  =  — - v  -  ■  - - -  2, 

cosa  I  ,/(fcyTsina)J +(X.-cosal2  | 


V. 

•\  'A 


From  Scaled  Surface  Equation : 

Ik.  -  cosai 

Letting  .  /(a)  =  ~  : - 

<J{ky  r  sine)2  +  -  cosa)2 

then  i  -  *  [/(a)  -  and  i~  ~/2(a) 

cosa 1  ‘ 

-  J 1  -/2( a)  -  |(i)'ana[/(a)  -  Z,]  +  yr |  +  K  =  0 
For  0  <  a  <  cos ~l{kz) 


162 


For  Case  II.B  k,  is  assumed  to  be  non-zero.  As  a  result,  the  final  expression  of 
Eq  (A.  19)  is  used  to  explicitly  solve  for  Tt,  leading  to  test  Condition  #1  shown  tn 
Eq  (A.26).  Also,  a  =  kx  implies  that  /3  can  be  explicitly  solved  for  and  the  range  of 


possible  a  values  restricted  as  shown  in  Eq  (A  27) 


For  Case  II.B :  xc  =  0  ,  yc  *  0  and  kx  *  0 


♦  *,  =  0  -  U  =  -J- 


where  0  s  (.  <  •*  — -  >  0  "  Condition  #1" 


(A.26) 


Since  d  =  kt  -  sin  a  cos  p  *  0 
then  p  =  cos'1 


K 

k  sin  a  , 

sin  a 

s  1 


(A.27) 


for  0  <  a  <  cos  ‘(^,)  ~  sin'1(|itJj)  <  a  <  cos'  '(*«) 


From  Scaled  Surface  Equation : 


(A.28) 

K  =  0 


The  derived  expressions  for  $  and  Tc  are  next  substituted  into  the  scaled  surface 
equation  resulting  in  an  expression  involving  only  o  ar  snown  m  Eq  (A.28).  Provided 
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an  a  value  exists  which  satisfies  this  expression,  two  additional  conditions  must  be 
satisfied  before  declaring  the  existence  of  a  valid  critical  point.  After  calculating  the 
critical  point  coordinates  using  the  previously  derived  values  of  a,  0,  and  7C,  a  check 
must  be  made  to  ensure  the  velocity  equations  are  satisfied.  In  doing  so,  X  is  first 
calculated  and  checked  to  ensure  it  is  negative.  The  right  hand  equality  of  Eq  (A.6) 
is  then  used  to  ensure  the  velocity  equations  are  satisfied  for  the  calculated  X  value. 
These  checks,  identified  as  Conditions  #2  and  #3  in  Eq  (A.29),  if  satisfied  ensure  the 
critical  point  coordinates  calculated  are  valid. 


Provided  a  ,  P  and  ic  exist: 

X  =  signlyl''j(yky  -  sin  a  sin  p)  <  0  “Condition  #2“ 


then  (*  -  ky)sign{  yc)  =  — - 


(A.29) 


-2  =  _ (C0Sg  -  kz] 


ky  -  sin  a  sin  Pj2  +  (cos  a  -  k^ 


"Condition  #3" 


For  the  last  major  case  to  be  considered,  Case  III,  it  is  assumed  that  xc  ¥=  0 
and  y ,  =  0  and  the  derivation  procedure  parallels  the  process  previous!}  used  for 
Case  II  Since  y  t  =  0,  the  pa.tial  derivatives  of  g(x,y,z)  wit!,  icspect  to  y  van'shes. 
As  such,  only  equality  of  the  first  and  last  expressions  of  Eq  VA.6J  is  valid  and  the 
velocity  expression  involving  T>  and  ky  is  used  to  establish  the  relanundup  sho*.<.  m 


Two  subcases  are  considered  as  a  result  of  this  relationship,  Case  III.A  for  fcy 
equal  to  zero  and  Case  III.B  for  /cy  not  equal  to  zero.  For  B  =  sinasin/3  =  ky  =  0, 
Case  III.A  is  further  subdivided  into  two  special  cases.  Case  III.A.1  provides 
solutions  based  cn  sina  =  0  or  both  sina  -  0  and  cos/3  =  0,  and  Case  III.A.2 
provides  solutions  based  on  cos/3  =  0  only.  Eqs  (A.31)  thru  (A.33)  summarize  the 
results  for  Case  III.A.1.  The  three  conditions  shown  are  derived  by  considering  both 
the  velocity  equations  and  surface  constraint  equation  as  shown.  Given  all  three 
conditions  are  satisfied,  a  critical  point  exists  at  x„  =  x„  yc  =  0,  and  zc  as  given  by 
Eq  (A.33). 


For  Case  111 :  xc  *  0  and  yc  =  0 


y£=o  -  g. 


=  0 


b-k}  -  Xgf 


0  -  5  =  k}  ”  yc  =  kytc  *  yr  =  0 


(A.30) 


For  Case  III  A. I :  ky  =  0  and  sina  =  0 

ky  =  0  —  yr  =  0  "Condition#  V 

(  a  =  0  —  xc  =  xr  (A  31) 

sina  =  0  —  < 

[  a  =  kti  ,  k  =  0,±l,i2,.;.. 

.-.  c  =  cosa  =  ±  1  but  c  >  kz>  0  —  c  =  1 
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From  The  Velocity  Equations: 


a-kx  =  Xg ,  -  kx  =  -  ksign^x,} 


\<0  ~ 


j  k,  <  0  and  x>0 


\k  > 0  and  x<0 


’  Condition  #2" 


For  X  =  sign[xr}kx 


c-k:  =  Xg:  -  zc  =  , 


/nto  the  Scaled  Surface  Equation 

1  +  k,  i  , 

,  - — -  -  i,  *  K  -  0  "Condition  #3" 

\  9  I  M 


For  Case  III  A  i,  a  detailed  development  procedure  which  parallels  the  process 
used  for  Case  II  A.2  is  required  to  obtain  critical  point  solutions.  Eqs  (A.34) 
thru  (A.36)  summarize  result?  of  the  development  process.  As  with  the  previous 
case,  the  condition  that  y,  =  0  must  again  be  satisfied.  Under  the  assumption  that 
sin/3  -  0  and  given  the  range  of  p  values  specified  in  Eq  (A.5),  only  two  /3  values  are 
possible,  0  and  rr,  as  indicated  in  Eq  (A  34).  Solving  for  X  from  the  velocity 
equations  allows  for  T,  to  be  explicitly  solved  tor  as  a  function  of  a  as  indicated  in 
Eq  (A.35).  Introducing  the  function  f(a),  the  Tc  expression  is  substituted  into  the 
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scaled  surface  equation  resulting  in  the  final  expression  shown  in  Eq  (A.36),  a  single 
expression  of  one  variable  a.  Provided  an  a  exists  satisfying  this  expression,  7.  can 
be  calculated  from  Eq  (A.35  j  and  the  critical  point  coordinates  determined  via 
Eq  (A.2). 

For  Case  II!  A  2:  ky  -  0  and  sin  P  =  0 

k=  0  -  y  =0  "Condition  H  V 
y  '  (A.34) 

sinp  -  0  ~  rasp  =  (±)1  ~  p  =  0  or  n 

—  a  -  sin  a  cos  P  =  (±)  sin  a 


From  Velocity  Equations : 


-  X  =  svg/tjfyjffcj  r  since 


c-k.  =  Xg, 


.1  f*.-COSC  /  7  -|2 

*  1-.  V  ,  .-.A*  t  c  'I 


(A.35) 


0  .  -  ^  [  _ ( V 


^(xjsinuf  r(kz- cosa 

(k.  -cosa) 

|  y/tX^rsina  i2  +  (£.  -  cosa  f  j 


-z, <  r. 


For  Case  Ili.B  f.,  is  assumed  to  be  non  zero.  As  a  resj'..  the  final  expression  of 
Eq  (A.30)  is  used  to  explicitly  solw  for  7t,  icadmg  to  test  Condition  #\  shown  in 
Eq  (A.37)  Also,  o  =  xv  implies  that  B  tan  be  explicitly  solved  for  and  the  range  of 
possible  a  values  restricted  as  shown  in  Eq  (A.38). 

! 
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From  Scaled  Surface  Equation : 

(k-  cosa) 

letting.  f(a)  =  - - - - - 

y  (it,  r  sin  a  j2  +  (it,  -cosa  )2 

then  f  =  — L-[/(a) -£,]  and  zc  =/2(a) 
cosa 1  J 

-  Vl  -/2(«)  -|(i)tana[/(a)-i]  +  i,|  +  K  =  0 

(A.36) 

For  0  <  a  <  cos "'(fc,) 

For  Cose  ///.B :  i,  *  0  ,  >c  -  0  and  ky  *  0 

y,  =  °  -  r>-f' 

r 

(A.37) 

m  —  y 

M’ftere  0  <  f  <  f  —  — -  >  0  " Condition  ft  1  ’ 

C  «!XU*)  r 

s 

5/nce :  o  -  k  =  sinasinp  '0 

(  Jt  |  Jt 

then  p  =  sin ' 1  —  —  !  —  — —  s  1 

\  since  )  sina 

(A.38) 

for  0  <  e  <  cos '  '■(<:,)  —  sin''(|trj)  <  a  <  cos''(jt.) 

The  derived  expressions  for  £  and  T,  arc  next  substituted  into  the  scaled  surface 

|  equation  resulting  in  an  expression  involving  only  a  as  sr..'wn  ir.  Eq  (A  39' 

Provided 

1  an  a  value  exists  which  satisfies  this  expression,  two  additional  conditions  must  be 

16ft 

I 

l _ 

satisfied  before  declaring  the  existence  of  a  valid  critical  point.  Aftei  calcuia'ing  the 
critical  point  coordinates  using  the  previously  derived  values  of  a,  (3,  and  7„  a  ihcck 
must  be  made  to  ensure  the  velocity  equations  are  satisfied.  In  doing  so,  X  is  first 
calculated  and  checked  to  ensure  it  is  negative.  The  left-hand  and  right-hand 
expressions  of  Eq  (A.6)  are  equated  and  used  to  ensure  the  velocity  equations  are 
satisfied  for  the  calculated  X  value.  These  checks,  identified  as  Conditions  #2 
and  #3  in  Eq  (A.4),  if  satisfied  ensure  the  critical  point  coordinates  calculated  are 
valid. 


i  l  -r 

1  '  ~J~cosa  *  ' 


From  Scaled  Surface  Equation : 

f 

-sine  cos 
k 


-y  f  (  f  i 

— -  sine  cost  sin'‘|  — — 
ky  (  f sina  / 


(A.39) 


-  K  =  0 


Provided  a  ,  p  and  ic  exist : 

X  =  sign^xc'j^kI  -  sin  a  cos  P )  <  0  " Condition  #2" 


then  (a  -  X,)s!gnjrf)  = 
(cose  - k.f 


g/rn 


(A.40) 


iX,  -  stnocosp  )*  v  ^eosa  -  X,}* 


" Condition  #3" 
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Appendix  B:  Aperture  Mechanical  Scanning 


Mechanical  scanning  of  a  planar  phased-array  is  considered  in  this  appendix. 
Mechanical  scanning  is  employed  as  a  "simple"  method  for  varying  the  main  beam 
pointing  direction  while  maintaining  a  relatively  symmetrical  radiation  pattern; 
typically  not  the  case  for  electronic  beam  steering  via  phase  tapering  across  the 
elements.  As  such,  mechanical  steering  may  potentially  reduce/isolate  some 
monopulse  boresight  errors  being  considered  during  analysis.  Given  the  desire  to 
scan  the  antenna  main  beam  to  an  arbitrary  location,  the  task  at  hand  is  to  determine 
the  scanned  location  of  all  elements  within  the  array  relative  to  the  spherical 
coordinate  origin,  the  point  about  which  the  antenna  array  pivots  The  goal  is  to 
derive  simplified  mathematical  expressions  foi  the  rotated/scanned  element 
coordinates  in  spherical  form,  negating  the  need  to  perform  a  more  computationally 
intense  rectanguiai  coordinate  rotation  followed  by  a  rectangular-to-spherical 
conversion,  a  process  which  would  be  required  for  each  antenna  element  at  each 
desired  scan  angle. 

xT  -  x  (cos^cosCo  +  sin2  ^  -  y  cos<f>0sin4>0(l  -  cos  9^ 

yr  -  -x  cos 4i„ sin t+>0 { 1  -  cosG^)  *  y  (sin2<J>0cos60  -  cos2 <J>0)  (B.l) 

zt  -  -x  cos4>osin0o  -  y  sin4>0sin00 


Previous  mechanical  scanning  derivations  typically  generated  and  utilized 
rectangular  coordinate  rotation  equations  to  arrive  at  expressions  for  the  rotated 
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element  locations.  The  rotated  rectangular  coordinates  given  by  Eq  (B.l)  are  an 
example  of  the  type  of  expressions  usually  derived,  where  6^  is  the  antenna  scan 
angle  from  the  z-axis  in  the  <pa  direction  as  measured  from  the  positive  x-axis  [12]. 

From  Eq  (B.l)  it  is  readily  apparent  that  use  of  these  equations  to  determine  the 
rotated  element  locations  in  spherical  coordinates  is  computationally  intense, 
requiring  four  trigonometric  function  calculations,  two  squaring  operations,  ;n  excess 
of  fifteen  multiplications,  additions,  and/or  subtractions,  as  well  as,  invocation  of  a 
rectangular-to-spherical  conversion  routine.  The  computational  intensity  is  further 
compounded  when  considering  two  additional  factors.  First,  the  calculations  may  be 
required  for  antenna  arrays  which  typically  are  composed  of  100  or  more  antenna 
elements.  Second,  scanning  the  antenna  array  in  one  plane  may  itself  require 
calculations  to  be  carried  out  for  a  large  number  of  discrete  antenna  scan  angles, 
possibly  requiring  100  or  more  to  obtain  accurate,  reliable,  and  consistent  output 
data.  As  such,  minimizing  the  computational  burden  while  increasing  efficiency  is  the 
key  motivating  force  behind  the  following  discussion. 

Detailed  information  is  provided  on  the  derivation  process  used  in  developing  the 
final  equations  for  determining  the  rotated  antenna  element  locations  in  spherical 
coordinates.  In  developing  ihe  equations,  the  apetture  is  initially  assumed  to  be 
located  >n  the  x-y  plane.  The  aperture  is  rotated 'scanned  about  the  origin  to  an  angle 

from  the  z-axis  in  a  direction  specified  by  $>.,  from  the  positive  x-axis  (this 
convention  is  consistent  with  the  convention  used  in  establishing  ihe  expressions 
at  Eq  (B.l;.  For  any  given  amentia  element,  initially  located  at  ihe  point  (r,,ir/2.v.) 
m  spnsncai  coordinates,  <■  final  clement  location  is  required  ar.d  will  hr  specified  by 


the  coordinates  (rr0o0,)  after  rutauon/scanmr.g.  For  the  specified  coordinate  system, 
the  following  variable  ranges  are  defined.  rt  >  0,  0  <  80  <  nQ  0  <  6,,  8,  <  and 
•jr  <  <p„  <p„  <0.  <  x.  The  restriction  on  the  maximum  scan  angle  of  0„  - -  ir/2  results 
from  considering  only  .^an  angles  such  that  rays  emanating  normal  to  -.he  aperture 
surface  intersect  the  rademe  surface  under  consideiat-on. 


y 


Figure  B  1  Case  !  Scanning  Geometry,  c.  a>^n 


The  development  process  pioceeds  by  considenr.g  two  special  cases  which  are 
do  inguishab'e  by  the  vaiuc  of  relative  to  0O  for  a  given  ..menna  e.ement  Case  1 
applies  to  all  antenna  elements  where  0„  <  0,.  <  cs„  +  7  and  Case  II  applies  'o  ail 
elements  where  0„  -  rr  <  <p,  <  0O.  The  geometry  established  for  considering  each 
case  is  as  illustrated  in  Figures  B.  1  and  B  3  for  Case  I  and  Case  U.  respectively,  rot 
both  of  these  figures,  the  following  definitions  arc  made. 


p  =  axis  perpendicular  to  the  2-axis  located  <p0  from  the  x-axis 

q  =  axis  perpendicular  to  the  2-axis  and  p-axis  about  which  rotation 
by  8  occurs 

p'  =  axis  parallel  to  p  along  which  the  tip  of  Fp  follows  as  0O  vanes 

Fp  =  vector  "projection'  of  r  m  the  x-y  plane  for  a  specified  0,,  value 

D  =  distance  from  the  origin  to  p'  as  measured  along  the  9-axis 
(Note  that  D  =  rK„ia)  which  occurs  for  0„  =  it  12) 

E  =  distance  from  the  9-uxis  to  the  tip  of  F  as  measured  along  the 
p'-axis 

E'  =  distance  from  the  9- axis  to  the  tip  of  F'  as  measured  along  the 
p'-axis 

a  =  angle  measured  from  9- axis  to  Fp 
Y  -  angle  measured  from  9-axis  tc  F 

The  distance  E'  is  established  by  observing  that  variation  in  8„  from  0  to  it  12 
causes  r  to  rotate  about  the  9-axis  in  a  clockwise  direction  as  viewed  in  the  ->-9 
direction,  in  doing  so,  the  tip  of  F  traces  out  a  quaner-cncle  ip  the  z-p'  plane  as 
shown  in  Figure  S.2,  The  quarter-circle  generated  will  lie  either  in  the  'Icwer-haif 
of  the  z-p'  plane,  the  case  shown  in  Figure  8.2  tor  -  nil  <  <pt  <  p0  7/2.  or  m 
the  " upper-half'  of  the  z-p '  plane  which  occurs  when  <&,  -  ir  <  0.  <  e>„  -  7  '2  or 
<pc  +  7t  <  9'  <  1-  it! 2.  in  either  case,  the  tip  of ' ’protection "  rp  lies  on  the  p  ‘  -axis 
at  a  distance  of  E‘  -  £ccs(8„)  from  the  9- axis  The  "upper-half Howet -half 
relationship  of  the  quarter-circle  will  be  taken  into  consideration  later  when 
developing  an  expression  for  9,  With  all  geomeiuw  relationships  established  and  the 
quantities  rt,  0,.  sp0,  and  6„  given,  the  derivation  prowess  for  computing  <>,  for  Case  I 
proceeds  as  follows 


" -fto*'  •  4.1 
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Figure  B.2  Geometry  For  Establishing  E,cc-(90)  Quantity 


For  Case  I  Eqs  {B.2;  thru  (B.4i  were  developed  by  considering  the  geometry 
establishes  m  Figure  3.1.  With  the  quantities  7,  D.  and  E  determined  an  expression 
for  a  is  calculated  by  considering  tne  geometry  presented  in  Figure  B  2  1  *u-  a 

expression,  given  by  Eq  (3.2 ),  ,s  usee  in  »onjunctier.  with  the  georr.etrv  .stabbshed 
in  Figure  £  .  to  obtam  tne  rotate;  scanned  expression  a=  goer  >  Eo  ■  B  0* 


E  =  jr  |  sin(y)  =  -|rj  cos (<}>,-<)>,) 


a  =  tan 


E'\  Ecos(^ 

—  =  tan  1  - —  =  tan  1 

lOj  D 


(B.4) 


-cosf^-rjg cos (6^  | 

sin(4>e-<}>(1) 


.  .  it  ,  n  (  -cosf4>« -<)>„)  cos  (0^)  ] 

4>  -  <&,  +  o:  =  4>n  +  —  +  tan  i -  V 

'  °  2  0  :  I  | 


(B.6) 


For  <t>„  <  <j>,  <  <(>»»  = 


The  process  used  in  deriving  the  c*>r  expression  for  Case  I  is  next  applied  for  \ 

“  j 

Case  II,  <p0  -  it  <  <p<  <  t*a,  with  geometrical  relationships  as  estaohshed  in  Figure  B3  | 

Eqs  (B.7)  thru  (B.lOj  summarize  the  results  of  applying  the  procedure  over  this  | 

region  with  the  final  <p,  expression  given  bv  Ej  iB.ll,  T.-.e  onservation  is  made  that 
Eqs  ( B.6)  and  iB  11  ..ry  only  ir>  tne  +  -  sign  preceding  the  tr/2  quantity,  and  the 

1 

m:r  ,s  sig  or  -the-  -*■  nunerr.tor  cenommator  ot  tne  arctangent  argument  j 

Applying  *h'  ..t<r  gen  tropertv  giver  n  (B  121  to  Eu  B  1 it  can  be  shown 
"i  Eqs^B.6)..  .  .  c:ua!  'de-icai  expressions  "herenre  'egardless  rf 

the  rela'r  Mish  b  be-«"  j  a,  -er  ettu-tior  .an  oe  asec  to  calculate  tne  unai 
0  f c  me  antenna  tes*  m.!'  .  tisK-c-u 


B. 


=  4>„  -  -  +  a 


* 

- -  tan 


I  cos (<j)t C 
I  -sin(4>,-' 


cosfG^) 


(B.11) 


For  h-n  <  <t>e<  <t>a 


|  =  tan-1  f— )  +  tr  (B.12) 

1-Bj  UfiJ 


\Fp\  =  D1  +  E2cos2(6J 

-  \j  I?  I  sin2 (<t>,  -  4>^)  +  jrj  cos2  (<]>,-(}>.,)  cos2  (0„)  (B-13) 

=  |fj  ^  1  -  cns2(4>(-<tg  sin2  (0^) 


The  next  task  is  to  calculate  0„  the  rotated/scanned  0  value  of  the  antenna 
element  being  considered.  The  first  step  in  determining  0,  is  to  calculate  the 
magnitude  of  Fp.  Referring  to  Figure  B  2,  the  magnitude  of  rp  is  obtained  by 
appivmg  Pythagorean’s  Theorem  to  the  triangle  shown,  resulting  in  the  final 
expression  as  given  by  Eq  (B.13).  This  magnitude  expression  is  valid  for  each  of  the 
two  cases  previously  considered,  regardless  of  the  relationship  between  <f>,  and  0C. 

Using  this  magnitude  expression,  the  plane  containing  Fp  and  7  is  considered  to 
calculate  A0„  the  change  in  the  antenna  element’s  0  component  resulting  from 
rotation  by  0O  about  the  <,  axis.  Considering  this  plane,  the  expression  for  AG,  is 
obtained  as  given  by  Eq  ^B.  14)  This  equation  is  based  on  the  fact  that  the 


1~? 


magnitude  of  r  rtmains  unchanged  hroughout  t  n  ,  "ci/stunninu  "rocess  t- 
|?p|  =  |  r  | ,  since  the  aperture  pivot'  •  jt 'he  or  gin 


cos(A8^i  =  if- 


Figure  B.4  Regions  For  D'-ermin.:  r  Aduiion/Suntraction  . .  _vc 

As  mentioned  previously,  the  rt  lation  of  r  an  >ut  tne  q-axi  as  80 1  am  rom 
0  to  jt/2  produces  a  quarter-circle  when  lies  esthe-  .n  ne  "upper  ’tall"  or  Von/"  ml!" 
of  the  z-p '  plane.  This  "upper-halfTltmir-^alf  re,a  10;  .mp  determines  wnernc  - 
value  of  A3(  given  by  Eq  (B  14}  is  to  he  ..Mrattec  *rr  m  or  aoucc  ne  mitn.  8 
value.  Since  all  antenna  elements  are  in  the  it  .n-  r.nor  i  auon,  ttieir  inn  u 
8,  component  is  equal  to  it!2.  Theretore.  an  elements  routed  .aiue  is  simplv 


I 
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~'2  ±  A0„.  T  e  question  is,  "When  should  A0t  be  added  or  subtracted?".  To  answer 
nr-  question  -on.  e>  ihe  diagram  shown  in  Figure  B.4.  Three  separate  cases  exist 
or  letc;-,  ed  ’he  an  enna  element’s  initial  $>c  value  relative  to  <p0. 

_  ji  -iuenr.a  cements  located  in  Region  A,  including  the  positive  p-axis,  the  tip 
*  7  trace--;  .  jt  a  quarter-circle  lying  :n  the  "lower-half'  of  the  z-p'  plane;  for  this 
egion  6,  =  -  1  *  A0..  For  antenna  elements  located  in  Region  B,  including  the 
-egative  p-axis,  tip  ,f  7  traces  out  a  quarter-circle  lying  in  the  "upper-half  of  the 
'  -lane,  u  his  region  0,  =  tr/2  -  A6t.  Lastly,  points  located  on  thea-axis  do  not 
experience  an-  -furge  in  0  value,  in  this  case  0,  =  0t  =  ir/2.  Coilsidering  these 
;  Lnal  expression  for  0,  is  obtained  and  may  be  summarized  as  given  by 


6r  =  f  I’]  cos’1  {  v'  1  -  cos2 (4>f -<(>„) sin^e^)  } 

[  (-0  for  0  s  |  <t>€  -  <J)0  |  s  — 

lake  \ 

j  (-)  M  |  s  |  -  <(><,  |  s  it 


(B.15) 
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where  -  ji  <  d>  -  4>  sn 

T  e  c 


Appendix  C:  Aperture  Electronic  Scanning 


Electronic  scanning  of  a  planar  phased  array  is  considered  in  this  appendix.  A 
phased  array  aperture  is  an  antenna  whose  mam  beam  maximum  direction  or  pattern 
shape  is  primarily  controlled  by  the  relative  phase  of  the  element  excitation  currents 
on  the  array  [13].  Assuming  the  elements  and  array  feed  network  form  a  passive 
reciprocal  structure,  the  phased  array  will  exhibit  identical  far-field  patterns  both 
transmission  and  receptions  modes  [14].  The  following  development  is  based  on  a 
receive  formulation,  consisteru  with  rnunc  pulse  receiver  tracking  requirements. 


Figure  C.l  Aperture  Electrical  Scanning  Geometry 

In  devi  >pmg  expressions  for  the  phase  weight  0„(m)  ot  an  arbitrary  array 
element  located  at  (x„,y„,z.),  far-field  conditions  are  assumed  As  such,  the  E-Field 
incident  on  the  array  is  considered  to  be  a  Uniform  Plane  Wave  (UPW)  As 
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established  in  Figure  Cl,  the  UPW  surface  is  given  by  f(x,y,z)  =  0  with  the  plane 
wave  itself  being  specified  by  E,.  Additional  vector  quantities  identified  in  Figure  C.l 
are  defined  as  follows: 

E0  =  Incident  E-Field  at  gimbal  point  (0,0, zj 

E„  =  Incident  E-Field  at  m'1  element  location  (x*,y0,z„) 

D  =  Vector  from  gimbal  point  to  reference  plane 
D„  =  Vector  from  element  location  to  reference  plane 
Tt  =  Vector  from  coordinate  origin  to  gimbal  point 
r„  -  Vector  from  coordinate  origin  to  m*  array  element 
Letting  0,  and  <p,  represent  the  spherical  coordinate  angles  which  define  the 

A 

desired  aperture  scan  direction,  unit  direction  vector  k  is  formed  as  indicated  m 
Eq  (C.l).  Since  k  represents  the  direction  from  which  the  aperture  response  is 
desired  to  be  a  maximum,  a  UPW  mciden*  fiom  this  direction  is  considered.  This 
condition  dictates  that  Vf(x,y,z)  =  k  for  all  (x^.z,)  on  the  reference  surface.  The 
resulting  expression  for  the  reference  surface  is  as  given  by  Eq  (C.2). 


Given  0  and  4> 


k  =  ksi  *  k  y  *  ktz  =  sinOjCos^i  *•  sin O^sin^y  +  cos0,z 


(C.1) 


F'r  V/(x,y,z)  =  k 
/(*,}’, z)  kxx  *  ky y  +  kt z  *  C 


(C.2) 


The  unknown  C  in  Eq  (C.2)  is  solved  for  by  noting  that  the  vector  sum  of  r, 
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and  D  must  lie  on  the  reference  surface.  This  constraint  allows  for  the  value  of  C 
to  be  solved  for  as  specified  in  Eq  (C.3),  which  when  substituted  into  Eq  (C.2)  results 
in  the  final  expression  for  f(x,y,z)  as  given  by  Eq  (C.4).  This  equation  completely 
describes  the  reference  surface  for  given  scan  angles  and  specified  gimbal  point 
location. 

Given  r  +  5  satisfies  f(x,y,z)  =  0 
f[Dkx , Dky ,  z,  -  Dkt)]  =  D(k]  +  *$+£)  +  ztkt  ♦  C  «  0 
-  C  =  -(O  +  zgkz) 

f(x,y,z)  =  kxx  *  kyy  <■  kt  (z  -  z,)  -  D  =  0 

Given  Fm  *  Dm  satisfies  f(x,y,z)  =  0 

/[(*«  +  D«K^  >  +  Dnky) ,  +  *>„*,)]  =  0 

-  (*mkx  *  ymky  +  z„kz)  *  Dm  {k\  *  k]  *  k\ )  -  (D  *  zfkz)  =  0 
-  (Dm-D)  --  lxnkx  +  ymkr(zn-Zt)k£\ 

For  a  pattern  max, mum  to  occur  at  given  scan  direction  k ,  responses  from  all 
elements  within  the  array  must  coherently  sum.  Since  a  far-fielu  assumption  is  made, 
all  rays  emanating  from  the  UPW  reference  surface  which  are  incident  on  the 
aperture  elements  are  approximately  parallel.  Therefore,  an  arbitrary  array  element 


(C.3) 


(C.4) 


(C.5) 
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may  be  chosen  for  developing  the  final  expression  for  <#>w(m).  Given  arbitrary 
element  location  (x„,y„„z„)  as  indicated  in  Figure  Cl,  the  vector  sum  of  r„  and  D„ 
is  constrained  to  lie  on  the  reference  surface  and  must  satisfy  Eq  (C.4).  Enforcing 
this  constraint  results  in  the  final  expression  relating  D„  and  D  as  given  by  Eq  (C.5). 
The  difference  between  D„  and  D  corresponds  to  phase  delay  variances  experience 
along  the  two  propagation  paths. 

Given  £  =  £;  e1*' 

Eo  -  -  jr  evl*.»-V4»l  (C.6) 

£  =  £^  f 

For  Coherent  Summation 

<t>K(o)  -=  K0Dm  -  *,  -  «J >»  (C.7) 

-  4>»  =  K,(Da  -  D )  *  4>„(o)  ,  K0(Dn  -  D ) 


Given  the  incident  E-Field  E„  expressions  for  the  E-Fields  incident  at  the  gimbal 
point  and  element  location  may  be  expressed  as  given  in  Eq  (C.6).  Free-space 
propagation  constant  A'c  is  introduced  to  convert  distance  variations  into  appropriate 
radian  phase  variations.  For  coherent  summation  to  occur,  the  final  phase  terms  of 
the  two  incident  E-Fields  in  Eq  (C.6j  must  be  equivalent.  Eq  (C.7)  shows  the  result 
of  equating  the  two  phase  terms.  In  simplifying  the  final  expression,  the  phase 
reference  for  all  aperture  elements  was  arbiirar.iy  chosen  to  be  the  aperture  center 
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at  (0,0, zt),  eliminating  0,(0)  from  the  final  expression.  Results  of  Eq  (C.5)  and 
definitions  established  in  Eq  (C.l)  are  substituted  into  Eq  (C.7)  to  arrive  at  the  final 
phase  weight  expression  given  by  Eq  (C.8).  Given  arbitrary  element  locations  and 
a  desired  scan  direction,  this  expression  is  used  to  determine  phase  weights  for  all 
elements  within  the  phased  array  aperture. 


4>»  =  -*„{*„*, +  y „  V  (z»  '  z*)ki} 

=  sine^cos*,  +  >„sm0j)  -  cos0s(z„  -  z,)} 

A0 


(C.8) 


184 


Bibliography 

1.  Sherman,  Samuel  M.,  Monopuke  Principles  and  Techniques.  Massachusetts: 
Artech  House,  1984. 

2.  Siwiak,  K.,  T.B.  Dowling  and  L.R.  Lewis,  "Boresight  Errors  Induced  by  Missile 
Radomes,"  IEEE  Trans.  Antennas  Propagat.,  AP-27,  Nov  ’79,  pp.  832-841. 

3.  Huddleston,  G.K.,  H.L  Bassett,  and  J.M.  Newton  "Parametric  Investigation 
of  Radome  Analysis  Results:  Salient  Results",  Vo!  I  of  IV,  Technical  Report, 
AFOSR-77-3469,  Bolling  AFB,  1981. 

4.  Burks,  D.G.,  E.R.  Graf,  and  M.D.  Fahey,  "A  High-Fiequency  Analysis  of 
Radome-lnduced  Radar  Pointing  Error,"  IEEE  Trans.  Antennas  Propagat., 
AP-30,  Sep  ’82,  pp  947-955. 

5.  Kozakoff,  D.J.  "Improved  Radome  Analysis  Methodology",  Conference  Notes, 
MICOM  Workshop  on  EM  Windows  and  Domes,  Redstone  Arsenal,  1981. 

6.  Klemer,  D.P.  "Effects  of  Refraction  by  an  Ogive  Radome  on  Radome 
Boresight  Error",  Technical  Report  #675,  Lincoln  Laboratory,  Massachusetts 
Institute  of  Technology,  1984. 

7.  Huddleston,  G.K.  and  H.L.  Bassett,  "Radomes,"  in  Antenna  Engineering 
Handbook.  R.  Johnson  and  H.  Jasik,  Eds.,  New  York:  McGraw-Hill,  1984. 

8.  Balanis,  C. A.,  Advanced  Engineering  Electromagnetics.  New  York:  John  Wiley 
&  Sons,  1989. 

9.  Ruck,  G.T.,  D.E.  Barrick,  W.D.  Stuart,  arid  C.K.  Krichbaum,  Radar  Cross 
Section  Handbook,  Volume  2.  New  York:  Plenum  Press,  1970. 

10.  Richmond,  J.H.  "Efficient  Recursive  Solutions  for  Plane  and  Cylindrical 
Multilayers",  Technical  Report  #1968-1,  Contract  #NOW  65-0329-d,  The 
Ohio  State  University  Research  Foundation,  1968. 

11.  Henderson,  L.W.  "Introduction  to  PMM",  Technical  Report,  WRDC-TR- 
89-9007,  Wright-Patterson  AFB,  1989. 

12.  McGrath,  D.T.,  "Analysis  of  Radome  Induced  Cross  Polarization”.  Technical 
Manual,  WL-TM-92-700-XPN,  Wright-Patterson  AFB,  1992. 

13.  Stutzman,  W.L,  G.A  Thiele,  Antenna  Theory  and  Design.  New  York:  John 
Wiley  &  Sons,  1981. 

14.  Stark,  L,  "Microwave  Theory  of  Phased-Array  Antennas,  A  Review,"  Proc.  of 
the  IEEE,  Vol.  62,  Dec  1974. 


185 


15.  Balanis,  C.A.,  Antenna  Theory:  Analysis  and  Design.  New  York:  John  Wiley 
&  Sons,  1982. 

16.  Leonov  A.I.  and  K.l.  Fomichev  (Translated  by  William  F.  Barton),  Monopulse 
Radar.  Massachusetts:  Artech  House,  1986. 

17.  Barton,  D.K.  and  H.R.  Ward  "Handbook  of  Radar  Measurement’,  New 
Jersey:  Prentice  Hall,  Inc.,  1969. 

18.  Young,  D.M.  and  R.T.  Gregory,  A  Survey  of  Numerical  Mathematics,  Vol  I 
of  II,  New  York:  Dover  Publications,  Inc. 

19.  Gelfand,  I.M.  and  S.V.  Fomin,  Calculus  of  Variations,  New  Jersey:  Prentice- 
Hall,  1963. 

20.  Munk,  B.A.,  G.A.  Burrell,  and  T.W.  Kornbau,  "General  Theory  of  Periodic 
Surfaces  in  a  Stratified  Dielectric  Medium,"  Technical  Report, 
AFAL-TR-77-219,  The  Ohio  State  University,  1977. 

21.  Tricoles,  G.,  "Radiation  Patterns  and  Boresight  Error  of  a  Microwave  Antenna 
Enclosed  in  an  Axially  Symmetric  Dielectric  Shell,"  Journal  Article,  Optical 
Society  of  America,  Vol.  54,  No.  9,  pp.  1094-1101,  Sep.  1964. 

22.  Tang,  R.  and  R.W.  Burns,  "Phased  Arrays,"  in  Antenna  Engineering  Handbook, 
R.  Johnson  and  H.  Jasik,  Eds.,  New  York:  McGraw-Hill,  1984. 

23.  Hirsch,  Herbert  L.  and  Douglas  C.  Grove,  Practical  Simulation  of 
Radar  Antennas  and  Radomes.  Massachusetts:  Artech  House,  1978. 

24.  Moeschhn,  L.  and  F.  Ehrensvard,  "Computer  Analysis  of  Radome 
Influence  on  Antenna  Radiation  Patterns  by  Surface  Integration 
Technique."  Proc.  of  The  Twelfth  Symposium  on  Electromagnetic 
Windows,  Jun  ’74,  pp.  8-11. 

25.  Orta,  R.,  R.  Tascone,  and  R.  Zich,  "Performance  Degradation  of 
Dielectric  Radome  Covered  Antennas,"  IEEE  Trans.  Antennas 
Propagation,  AP-36,  Dec  1988,  pp.  1707-1713. 

26.  Rhodes,  D.R.,  Introduction  to  Monopulse.  New  York:  McGraw-Hill, 

1959. 


186 


REPORT  DOCUMENTATION  PAGE 


CM8  Vo  0704  0738 


1  AGENCY  USE  ONLY  Leave  ounx) 


3  REPORT  TYPE  AND  OATES  COVERED 
Dissertation 


4  TITLE  AND  SUBTITLE 

RADOME  DEPOLARIZATION  EFFECTS  ON  MONOPULSE 
RECEIVER  TRACKING  PERFORMANCE 


|  6  AUTHOR(S) 

!  Michael  A.  Temple,  Capt,  USAF 


1  7  PERFORMING  ORGANIZATION  NAMt(S)  AND  ADORESS(ES) 


j  Air  Force  Institute  of  Technology,  WPAFB,  OH  45433-6583 


9  SPONSORING.  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFTT/DS/EN  G/93-03 


tO  SPONSORING  MONITORING 
AGENCY  REPORT  NUMBER 


David  A.  Hime 

Wright  Laboratory,  AAWD 

WPAFB,  OH  45433-6534 


12a  DISTRIBUTION /AVAILABILITY  STATEMENT  I  12b  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  unlimited 


13.  ABSTPACT  (Maximum  200  words) 

Borcsight  Error  (BSE),  defined  as  the  angular  deviation  between  the  true  position  and  the  apparent 
position  of  a  target  as  indicated  by  a  radar,  is  an  important  figure  of  merit  for  a  tracking  radar.  A 
significant  contributor  to  system  BSE  is  the  protective  radome.  This  research  effort  employed  a  GO 
technique  to  investigate  the  effects  of  a  radome  on  BSE,  expanding  previous  ray-trace  receive  techniques 
to  include:  1)  a  uniquely  defined/developcd  mathematical  description  for  each  surface  within  arbitrary 
multi-layer  tapered  radomes,  2)  an  "ideal*  taper  function  concept  for  obtaining  optimum  BSE  prediction 
performance,  3)  a  generalized  technique  for  calculating  specular  reflection  points  within  the  radome,  and 
4)  the  total  refractive  effects  along  ray  propagation  paths.  Computer  model  results  were  compared  with 
limiting  case  data  (BSE  =  0°),  published  experimental  data,  and  production  system  acceptance  test  data. 
For  all  limiting  cases,  ’system"  modeling  error  was  less  than  .06  mRad.  "Excellent*  (BSE  within  ±1  mRad) 
results  were  obtained  using  a  hemispheric  radome  with  a  displaced  aperture  gunbal  point,  predicted  BSE 
values  were  within  ±  1  mRad  of  published  experimental  data.  Likewise,  BSE  predictions  for  the  production 
system  were  within  ±0.5  mRad  of  measured  data  over  a  30*  scan  range.  Ray  refractive  effects  on  BSE 
prediction  were  characterized  using  the  validated  model. 
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